Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA40.0% 46.1 / 19.1 41intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (21/3);  new base A (0/16);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.54e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001763‑2002059
                                                                                                                                                       |                                                                                                                                                    
acCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGa                                                                                                                                                     >  2:658576/1‑150 (MQ=255)
 ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGaaa                                                                                                                                                   <  1:302448/151‑1 (MQ=255)
            ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCccgcc                                                                                                                                        >  2:15791/1‑151 (MQ=255)
            ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCccgcc                                                                                                                                        >  1:145093/1‑151 (MQ=255)
                agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTg                                                                                                                                    >  2:346284/1‑151 (MQ=255)
                   gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                 >  2:401547/1‑151 (MQ=255)
                            aaaaaCGGCCCTGCTTTATATGAGTTAGCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAAGCGGACGATTAGTGGGTAAAAGGAGGGTTTATTGGGGGGTTAAAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCa                                                                                                                         <  2:477292/148‑1 (MQ=255)
                              aaacgccctggctttATATGATTTATCGGGATGTTTATCCGTTTCTGCAGGTTTTTAAATCTGACGTTTAGTGGGTAAAAGGAGGGGTTATTGGGGGGTTAAAGG‑AAATT‑CCC‑GGACAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                        <  1:516297/136‑1 (MQ=255)
                              aaaCAGCCCTGCGTTATATGAGTTACCGACATAATTATCCGTTTCTGGAGGTTTTTTAAGCGGACGATTGGTGGGTGAAAGGAGGGGTTTTTGGGGGGTTAAAGGTAAATC‑CCA‑GC‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                          <  1:166032/147‑1 (MQ=255)
                               aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                     >  2:188961/1‑151 (MQ=255)
                                         cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                           >  1:586158/1‑151 (MQ=255)
                                              tatGAGTTTCCGGATTGATTATCCGTTTCTGCGGGTTTTTAAATCGGACGATTAGTGGGTGAAAGGGGGGTTTATTGGGGGGTTACAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                        <  1:376965/149‑1 (MQ=255)
                                                         ggCATGTTTCTCCTTTTCCGGAGGTTTTTTAAGCGAACGATTAGTGGGTAAAAGGGGGGGTTTTTTGGGGGTTAC‑GTTAAATTCCCA‑GT‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                             <  1:73961/149‑1 (MQ=255)
                                                                       tttCTGCGGGGTTTTTAATCGGACAATTGGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGGAAAATC‑CCA‑GG‑CAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                 <  1:17880/147‑1 (MQ=255)
                                                                       tttCTGCAGGGTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTTTTTGGGGTTTACAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                               <  1:222736/149‑1 (MQ=255)
                                                                       tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                              >  2:533730/1‑150 (MQ=255)
                                                                       tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                              >  2:84982/1‑150 (MQ=255)
                                                                            gcggggTTTTTAAGCGGAGGATTAGTGGTTGAAATGGGGGTTTTTTGGGGGGTTAAAGGAAATTC‑CCA‑GA‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                        <  1:29070/148‑1 (MQ=255)
                                                                              aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       >  1:206263/1‑150 (MQ=255)
                                                                               ggttttttAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       <  2:262905/146‑1 (MQ=255)
                                                                                                 taagggggTGAAAGGAGGGGTTATTTGGGGGTAAAAGGAAAATT‑CCA‑GG‑CACAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                           <  1:363295/138‑1 (MQ=38)
                                                                                                  tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  1:460968/1‑151 (MQ=255)
                                                                                                   aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  1:60281/1‑150 (MQ=255)
                                                                                                    gTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                 >  2:370589/1‑150 (MQ=255)
                                                                                                        gTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                             >  1:686972/1‑150 (MQ=255)
                                                                                                          gAAATGGGGTTTTTTTTGGGGTTTACAGTTAAATT‑CCA‑GA‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                   <  1:710263/142‑1 (MQ=37)
                                                                                                             aTGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                        >  2:344171/1‑150 (MQ=255)
                                                                                                               gAGGGGTTTTTGGGGGGTTACAGTTAAATC‑CAA‑GA‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     <  1:401455/151‑1 (MQ=37)
                                                                                                                            ggggtttaaagaaaaTT‑CCCGGG‑CAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                               <  1:360697/133‑1 (MQ=25)
                                                                                                                            ggggTTACAGGTAAATC‑CCA‑GG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                        <  2:481829/151‑1 (MQ=255)
                                                                                                                             ggttaaaaGGAAAATC‑CCA‑GG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                       <  1:553801/144‑1 (MQ=25)
                                                                                                                              gtttACAGGAAAATT‑CCG‑GG‑CAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      <  2:320318/149‑1 (MQ=255)
                                                                                                                                ttACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGAGCTTTATGTTGCCGGAAGCGGCGCAAACGCCTTATCCGGCCTGCAAAAATGTGCAAATTGAATAAATTGCAATTCAAc                          >  1:520384/1‑145 (MQ=255)
                                                                                                                                   cAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                 >  2:129710/1‑151 (MQ=255)
                                                                                                                                    aGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                >  2:599846/1‑151 (MQ=255)
                                                                                                                                     ggaaaaTC‑CCA‑GG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg               <  2:111897/148‑1 (MQ=255)
                                                                                                                                     ggTAAATT‑CAA‑GC‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                <  1:658489/150‑1 (MQ=255)
                                                                                                                                      gTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg               >  1:580641/1‑150 (MQ=255)
                                                                                                                                        aaaTC‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa            <  2:180189/151‑1 (MQ=255)
                                                                                                                                         aaTT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa            >  2:259527/1‑150 (MQ=255)
                                                                                                                                                      aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  2:594506/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001763‑2002059

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑