Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,972 T→G 17.1% intergenic (‑122/+35) fliL ← / ← fliK flagellar protein FliL/flagellar hook‑length control protein FliK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,9720TG17.1% 86.0 / 9.0 40intergenic (‑122/+35)fliL/fliKflagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base G (3/4);  total (18/22)
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.64e-01

CTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGC  >  NC_002947/4950835‑4951108
                                                                                                                                         |                                                                                                                                         
cTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                 <  7:103006/148‑1 (MQ=255)
  tCGCTCTTCGCCATGCCAATAATCCGTCGGCCATCGGGGGTTCAAAGGAGAAGGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGAGGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                              <  7:277595/149‑1 (MQ=255)
      tctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagc                                                                                                                            >  8:266944/1‑147 (MQ=255)
        ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagc                                                                                                                        >  1:353963/1‑149 (MQ=255)
        ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagc                                                                                                                        >  5:204560/1‑149 (MQ=255)
                    ataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAg                                                                                                           >  7:201454/1‑149 (MQ=255)
                    ataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAg                                                                                                           <  7:266944/149‑1 (MQ=255)
                       atCCGGCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTCCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGtt                                                                                                         <  5:305206/148‑1 (MQ=255)
                           gtcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAg                                                                                                    >  3:129751/1‑149 (MQ=255)
                                ccATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACa                                                                                                            >  3:323762/1‑136 (MQ=255)
                                ccATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACa                                                                                                            <  4:323762/136‑1 (MQ=255)
                                     gggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTaa                                                                                           <  3:236392/148‑1 (MQ=255)
                                                      tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGcc                                                                          >  7:161247/1‑148 (MQ=255)
                                                         cttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGAc                                                                       >  7:241773/1‑148 (MQ=255)
                                                            ggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACaa                                                                     <  8:201454/147‑1 (MQ=255)
                                                            ggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAac                                                                    <  2:64103/148‑1 (MQ=255)
                                                             gcttgAGCATGTGTTATGCCAGATTTTGCAGGGGTTGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAacc                                                                   <  7:295944/148‑1 (MQ=255)
                                                             gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAacc                                                                   <  4:129751/148‑1 (MQ=255)
                                                                gcagaaaGTGTTTTCCCATTTTTTCCCGGGTAGGGATTCCAGCCGGTTTTGGGGTGCGGGCTTCTTCGGGGCCGCCCCGGCTAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAaccacc                                                                <  2:353963/143‑1 (MQ=255)
                                                                  aGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctg                                                              <  8:241773/148‑1 (MQ=255)
                                                                   gCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCGGAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAac                                                                    <  1:254690/142‑1 (MQ=255)
                                                                   gCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAac                                                                    >  2:254690/1‑141 (MQ=255)
                                                                               ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                    <  5:139090/145‑1 (MQ=255)
                                                                               ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                    >  6:139090/1‑145 (MQ=255)
                                                                               ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                >  7:81326/1‑149 (MQ=255)
                                                                                 aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                              <  8:161247/149‑1 (MQ=255)
                                                                                  gATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCggg                                             <  6:204560/149‑1 (MQ=255)
                                                                                       tGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCt                                         <  8:81326/148‑1 (MQ=255)
                                                                                                        agcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGCAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCACAGGCCGGGCGCTATCGACCGCTTCGGcccc                       >  6:82456/1‑149 (MQ=255)
                                                                                                                             ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgc                                            <  4:134026/107‑1 (MQ=255)
                                                                                                                             ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgc                                            >  3:134026/1‑107 (MQ=255)
                                                                                                                             ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTg                 >  4:219244/1‑134 (MQ=255)
                                                                                                                             ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTg                 <  3:219244/134‑1 (MQ=255)
                                                                                                                             ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGAGCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc  >  2:191954/1‑149 (MQ=255)
                                                                                                                             ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAccc                   <  7:180444/132‑1 (MQ=255)
                                                                                                                             ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAccc                   >  8:180444/1‑132 (MQ=255)
                                                                                                                             ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCACAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAc                                   <  7:91095/116‑1 (MQ=255)
                                                                                                                             ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCACAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAc                                   >  8:91095/1‑116 (MQ=255)
                                                                                                                                        cTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTg                 >  6:292537/1‑123 (MQ=255)
                                                                                                                                        cTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTg                 <  5:292537/123‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGC‑CAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGC  >  NC_002947/4950835‑4951108

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: ‑