Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2940TC44.4% 39.6 / 12.6 36intergenic (+100/+129)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (3/17);  new base C (16/0);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.35e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  NC_000913/3423147‑3423443
                                                                                                                                                   |                                                                                                                                                       
ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCGAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                      <  1:137717/151‑1 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGCCCGCCCCCCtaccgtc                                                                                                                                                   >  2:422607/1‑145 (MQ=255)
     aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGTCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCg                                                                                                                                                 <  2:488277/151‑1 (MQ=255)
       cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCggt                                                                                                                                               >  1:490717/1‑149 (MQ=37)
           gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                           <  1:411370/151‑1 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCACCCCCCCGGCCGGTGGTTTGCAACAtt                                                                                                                                 >  1:565765/1‑148 (MQ=255)
                         ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTTGCAACGTTCaaa                                                                                                                             >  1:708242/1‑150 (MQ=255)
                               cGCAGGTAATCCATTAATTGAATTTTAGTTCGAAAAGAAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                       <  2:176541/150‑1 (MQ=255)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTTTCCCACACCCCCCAGGCCCGTGATTCGCTACGCTATAATCTGCTGCCcgg                                                                                                                 >  1:544845/1‑150 (MQ=255)
                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTTCTTTGCCAACTTTAAAACCCCCCCCCGCGGGttt                                                                                                          >  2:375797/1‑151 (MQ=255)
                                             tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCTCCCGGCCGTTTCTTCGCAAAGTTTAAAACCCCCCCccgcggggtt                                                                                                          >  1:126513/1‑144 (MQ=255)
                                              aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGCTGGTTCTCCAACATTCAAATCTCCCCCcgggggtgt                                                                                                           >  2:493710/1‑144 (MQ=255)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGGCGTTGCTTCCCCACGTTCAAATCCCCCCCCCGCGGGttttt                                                                                                        >  1:60524/1‑144 (MQ=18)
                                                       ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg                                                                                              <  1:568577/150‑1 (MQ=255)
                                                          gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAAAGTTCAAAACCCCCCCCCGCGGGGTTTTTCTACTCggg                                                                                              >  1:426826/1‑147 (MQ=17)
                                                            tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCACCCCCCCGGGCGGTGCTTCCCAACGGTCAAAATCGCTCCccgccgggttttccctc                                                                                                   >  1:285645/1‑127 (MQ=255)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGTTGCTTCGCCACGTTCCAATCCGCCCCCCGCGGGTTTTTTCTTCTCcgg                                                                                              >  2:624078/1‑144 (MQ=255)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc                                                                                  <  1:159967/151‑1 (MQ=255)
                                                                   gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCg                                                                                 <  1:667363/151‑1 (MQ=255)
                                                                    cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCCCAACGTTCAAAACCCCCCCCCGGGGGTTTTTTCTACTcccgggggggttcac                                                                                   >  1:161424/1‑134 (MQ=9)
                                                                     aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac                                                                               <  2:181257/151‑1 (MQ=255)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                             <  1:10450/151‑1 (MQ=255)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                             <  1:234260/151‑1 (MQ=255)
                                                                         aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                           <  1:432542/151‑1 (MQ=255)
                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGGACTGCCCCCCCCCCCCCGGCCCCGTTCTTTCCCACGCTTCAATCCGCTCCCTGCGGGTTTTTTCTTCTCCGGGGAGCGGTGACCGACAAaaaacaaa                                                                     >  1:247696/1‑146 (MQ=255)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa                                                                  <  1:533593/150‑1 (MQ=255)
                                                                                             cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCGGTTCTTTGCAACATTCAAATCCCCCTCCCGCGGGTTTTTTCCACTCCGCAAAGTGGTCACCGACAAAAAACAGATAaaaaaaaaaac                                                        >  1:333243/1‑147 (MQ=9)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGGCGGTGCTTTGCCACGTTCCAATCCGCTCCCCGCGGGGTTTTTCTCCTCGGGGGGGGGTTTCCCGCCACACAACAGATAAAAAAAAAAGCGCAggc                                                  >  2:291938/1‑149 (MQ=9)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTCCTGCCCCCCCCCCCTCCGGCCCTTGCTTCGCAAACGTCAAATCCCCTCCCGGCGGATTTTTTCCACTCGGGAGAGGGTTCACCGACAAACAACAAATAAAAAAAAAAGCCCAGTc                                                  >  1:238718/1‑151 (MQ=255)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                              <  2:67518/151‑1 (MQ=39)
                                                                                                        tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                              <  2:589558/149‑1 (MQ=38)
                                                                                                         ttgatGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCCCCCCCCGGGGCGGTTCTTCGCCAACTTCAAATACCCCCCCCGCGGGGTTGTTCTACTCCGGGGGGGGGTTACCGACCAAAAACAGATAAAAAAAAAAGGCCAGCCTTcctcc                                           >  2:534130/1‑148 (MQ=255)
                                                                                                                         tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                            <  2:515605/150‑1 (MQ=38)
                                                                                                                          aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                           <  2:489646/150‑1 (MQ=38)
                                                                                                                           tGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                         >  2:255381/1‑151 (MQ=38)
                                                                                                                                           cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc         <  2:66243/151‑1 (MQ=40)
                                                                                                                                                  cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc  >  2:672817/1‑151 (MQ=40)
                                                                                                                                                   |                                                                                                                                                       
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  NC_000913/3423147‑3423443

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