Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,965,196 G→T 54.5% C54* (TGC→TGA fleQ ← transcriptional regulator FleQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,965,1960GT54.5% 2.2 / 42.7 33C54* (TGC→TGAfleQtranscriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base T (9/9);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTG  >  NC_002947/4965051‑4965343
                                                                                                                                                 |                                                                                                                                                   
aaaGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAgcactcac                                                                                                                                                  <  1:55813/149‑1 (MQ=255)
 aaGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCTGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGcaca                                                                                                                                                 <  3:262788/149‑1 (MQ=255)
          gcggcgAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTcgc                                                                                                                                        <  2:12476/149‑1 (MQ=255)
                    tCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAATCCGGATTGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCTGCTCGAGGa                                                                                                                              <  7:103479/149‑1 (MQ=255)
                    tCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGa                                                                                                                              <  2:74290/149‑1 (MQ=255)
                     cTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGa                                                                                                                              >  1:26262/1‑148 (MQ=255)
                                  ggcggAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCg                                                                                                                >  5:239864/1‑149 (MQ=255)
                                       aagaaATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCctgc                                                                                                           >  1:122252/1‑149 (MQ=255)
                                         gaaATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCctgctg                                                                                                         <  2:122252/149‑1 (MQ=255)
                                         gaaATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCct                                                                                                             <  6:31890/145‑1 (MQ=255)
                                         gaaATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCct                                                                                                             >  5:31890/1‑145 (MQ=255)
                                                   cTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCAtgg                                                                                               <  2:26262/149‑1 (MQ=255)
                                                   cTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCGTGTGCCAGCCAtgac                                                                                              >  4:176403/1‑147 (MQ=255)
                                                          cAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGCCATGACTGGcag                                                                                        <  6:239864/149‑2 (MQ=255)
                                                                  gaaggaaCTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCCCGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAgg                                                                                <  2:154590/149‑1 (MQ=255)
                                                                  gaaggaaCTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCCGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAgg                                                                                <  1:220259/149‑1 (MQ=255)
                                                                      gaaCTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGttt                                                                             >  3:185473/1‑148 (MQ=255)
                                                                               ccAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAgcacgca                                                                                                                                                   >  6:165054/1‑69 (MQ=255)
                                                                               ccAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAgcacgca                                                                                                                                                   <  5:165054/69‑1 (MQ=255)
                                                                               ccAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCCGTCATGGCTGGCACCCGCGAGGTTGTCTTCGCCgc                                                                   >  8:169701/1‑148 (MQ=255)
                                                                                 aGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTAACAGAACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAgg                                                                 <  7:169701/149‑1 (MQ=255)
                                                                                         aCTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAg                                                         >  5:92427/1‑149 (MQ=255)
                                                                                         aCTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAg                                                         >  3:101008/1‑149 (MQ=255)
                                                                                               ttAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGcaccac                                                   >  7:27014/1‑149 (MQ=255)
                                                                                                  agcagcCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGCGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCg                                                 >  2:88531/1‑148 (MQ=255)
                                                                                                   gcagcCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAg                                                                  <  5:110313/130‑1 (MQ=255)
                                                                                                   gcagcCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAg                                                                  >  6:110313/1‑130 (MQ=255)
                                                                                                      gcCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAgg                                            <  6:92427/149‑1 (MQ=255)
                                                                                                        ccAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGca                                                       >  5:45677/1‑136 (MQ=255)
                                                                                                        ccAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGca                                                       <  6:45677/136‑1 (MQ=255)
                                                                                                          aaGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTcgc                                        >  4:237288/1‑149 (MQ=255)
                                                                                                                         aGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGGGCGCGCGCTgtc                        <  4:101008/149‑1 (MQ=255)
                                                                                                                           cATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtcg                       >  2:94512/1‑149 (MQ=255)
                                                                                                                           cATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTgtcg                       <  1:94512/149‑1 (MQ=255)
                                                                                                                                    tGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAgg                                            <  2:126816/119‑1 (MQ=255)
                                                                                                                                    tGCCGATGAGCACTCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAgg                                            >  1:126816/1‑119 (MQ=255)
                                                                                                                                               acgcacAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCtt   >  3:207511/1‑149 (MQ=255)
                                                                                                                                                cgcacAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTg  <  4:185473/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAAGCACCCGGCGGCGAAGGTCTTCCGGCAGCTCGGCGGAAGAAATTTCACCTAAAAGCAAAACCGGAAGGAACTCATCCCAGGTAGCCACTGTCTTAAGCAGCCCAAGCAGATTGCCCGGAGCATTCACGGTGCCGATGAGCACGCACAGCACTTCGCGGCTCGAGGACAATGGCTCGACCGCCTGCTGCCAGTCATGGCTGGCACACGCGAGGTTTTCTTCGCCGAGGAAGTTCAGCACCACCGCAAGGTCGCGGCGGCGCGCGCTGTCGTCATCGATCAGCAGAATCTTG  >  NC_002947/4965051‑4965343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: