Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC31.4% 103.6 / 15.9 51intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base A (34/1);  new base C (0/16);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 2.37e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG  >  NC_000913/2001784‑2002053
                                                                                                                                       |                                                                                                                                       
aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                    >  1:484169/1‑140 (MQ=255)
                agagagTTATATGGTTTACCGGTATGATTATCCTTTTCGGCGGGTTTTTAAACCGGGCGATTAGGGGGTGAAAGGGGGGTTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                    <  2:1602405/135‑1 (MQ=255)
                  tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc                                                                                                                  >  1:578321/1‑140 (MQ=255)
                  tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc                                                                                                                  >  1:456397/1‑140 (MQ=255)
                    ttttATAGGAGTTATGGGGAGGTTTATCCTTTTCTGGGGGTTTTTTAATGGGACGATGAGGGGTGGAAAGGAGGGGTTATTGGGGGT‑TAAAAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                <  2:303424/138‑1 (MQ=255)
                    cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                >  2:813999/1‑140 (MQ=255)
                       tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                             >  1:1240824/1‑140 (MQ=255)
                              gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  2:1102833/1‑140 (MQ=255)
                              gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  2:1025839/1‑140 (MQ=255)
                              gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  1:171565/1‑140 (MQ=255)
                              gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  1:1505678/1‑140 (MQ=255)
                              gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGCGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  1:303424/1‑140 (MQ=255)
                               ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGAGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                     >  1:480065/1‑140 (MQ=255)
                               ttACCGGAATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTTTTGGGGGG‑TTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                     <  2:456397/140‑1 (MQ=255)
                                 aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAcctgtct                                                                                                   >  1:1602405/1‑134 (MQ=255)
                                 aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                   >  2:426137/1‑140 (MQ=255)
                                   cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                 >  1:808038/1‑140 (MQ=255)
                                                    tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                >  1:584696/1‑140 (MQ=255)
                                                     ctcgaGGTTTTTTAACGGGACGTTTAGTGGGAAAAATGGGGGGTTATGTGGGGG‑TTAAAGATAAATTCCAGGAAGAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                               <  2:1458363/136‑1 (MQ=255)
                                                        cAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                            >  1:743884/1‑140 (MQ=255)
                                                            ttttttAATCGGACGATTGGGGGGTAAAAGGGGGGTTTTTTGGGGGG‑TTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                        <  2:584696/139‑1 (MQ=255)
                                                            gTTTTAAACGGGACGTTTGGTGGGTGAAAGGAGGGGTTATTGGGGGG‑TTACAGGAAAATTCCGGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                        <  2:578321/140‑1 (MQ=255)
                                                             tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                       >  1:372586/1‑140 (MQ=255)
                                                                ttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                    >  2:167144/1‑140 (MQ=255)
                                                                      ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                              >  2:669106/1‑140 (MQ=255)
                                                                       gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                             >  1:345804/1‑140 (MQ=255)
                                                                       gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                             >  1:644918/1‑140 (MQ=255)
                                                                       gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                             >  1:300602/1‑140 (MQ=255)
                                                                        aCGATTAGTGGGTAAAAGGGGGGTTTTTTGGGGGG‑TTACAGTAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                            <  2:300602/140‑1 (MQ=37)
                                                                               gggggTGAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                     <  1:919718/138‑1 (MQ=38)
                                                                                      aaaTGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                              >  1:1337500/1‑140 (MQ=255)
                                                                                         ggAGGGTTTATTGGGGGG‑TAAAAGAAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                           <  2:1627315/139‑1 (MQ=37)
                                                                                          gAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                          >  1:1627315/1‑140 (MQ=255)
                                                                                                 tATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                   >  1:1096526/1‑140 (MQ=255)
                                                                                                 tATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                   >  1:1229206/1‑140 (MQ=255)
                                                                                                  ttttGGGGG‑TTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                  <  2:1443433/139‑1 (MQ=255)
                                                                                                   tttGGGGGTTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                  <  1:147637/140‑1 (MQ=255)
                                                                                                     gggggg‑TTACAGGAAAATCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                               <  2:1337500/139‑1 (MQ=37)
                                                                                                       gggg‑TTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                             <  1:1025839/140‑1 (MQ=37)
                                                                                                         ggtttACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                            <  1:1344542/137‑1 (MQ=255)
                                                                                                             tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                        >  2:1612207/1‑140 (MQ=255)
                                                                                                              aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                       >  1:125132/1‑140 (MQ=255)
                                                                                                              aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                       >  1:964723/1‑140 (MQ=255)
                                                                                                                    aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                 <  1:669106/140‑1 (MQ=255)
                                                                                                                    aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                 >  1:390768/1‑140 (MQ=255)
                                                                                                                    aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                 <  2:1137179/140‑1 (MQ=255)
                                                                                                                    aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                 <  2:372586/140‑1 (MQ=255)
                                                                                                                    aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                 <  2:743884/140‑1 (MQ=255)
                                                                                                                     aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                >  1:1230307/1‑140 (MQ=255)
                                                                                                                       ttCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg              >  1:375302/1‑140 (MQ=255)
                                                                                                                         ccAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc            <  1:426137/140‑1 (MQ=255)
                                                                                                                          cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc           >  1:1595602/1‑140 (MQ=255)
                                                                                                                            ggCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg         >  2:1152273/1‑140 (MQ=255)
                                                                                                                                 aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag    >  2:649188/1‑140 (MQ=255)
                                                                                                                                   aaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg  <  2:480065/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG  >  NC_000913/2001784‑2002053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑