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

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  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)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCACCCCCCCGGCCGGTGGTTT‑GCAACAtt                                                                                                                                 >  1:565765/1‑148 (MQ=255)
                         ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTT‑GCAACGTTCaaa                                                                                                                             >  1:708242/1‑150 (MQ=255)
                               cGCAGGTAATCCATTAATTGAATTTTAGTTCGAAAAGAAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCt                                                                                                                       <  2:176541/150‑1 (MQ=255)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTTTCCCACACCCCCCAGGCCCGTGATTC‑GCTACGCTATAATCTGCTGCCcgg                                                                                                                 >  1:544845/1‑150 (MQ=255)
                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTTCTTT‑GCCAACTTTAAAACCCCC‑CCCCGCGGGttt                                                                                                          >  2:375797/1‑151 (MQ=255)
                                             tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCTCCCGGCCGTTTCTTC‑GCAAAGTTTAAAACCCCC‑CCccgcggggtt                                                                                                          >  1:126513/1‑144 (MQ=255)
                                              aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGCTGGTTCTCCAACATTCAAATCTCCC‑CCcgggggtgt                                                                                                           >  2:493710/1‑144 (MQ=255)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGGCGTTGCTTC‑CCCACGTTCAAATCCCCC‑CCCCGCGGGttttt                                                                                                        >  1:60524/1‑144 (MQ=18)
                                                       ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCggg                                                                                              <  1:568577/150‑1 (MQ=255)
                                                          gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTC‑GCAAAGTTCAAAACCCCC‑CCCCGCGGGGTTTTTCTACTCggg                                                                                              >  1:426826/1‑147 (MQ=17)
                                                            tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCACCCCCCCGGGCGGTGCTTC‑CCAACGGTCAAAATCGCT‑CCccgccgggttttccctc                                                                                                   >  1:285645/1‑127 (MQ=255)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCCGCCGTTGCTTC‑GCCACGTTCCAATCCGCC‑CCCCGCGGGTTTTTTCTTCTCcgg                                                                                              >  2:624078/1‑144 (MQ=255)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc                                                                                  <  1:159967/151‑1 (MQ=255)
                                                                   gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCg                                                                                 <  1:667363/151‑1 (MQ=255)
                                                                    cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCCCCCCCCCGGCCCGTGCTTC‑CCAACGTTCAAAACCCCC‑CCCCGGGGGTTTTTTCTACTcccgggggggttcac                                                                                   >  1:161424/1‑134 (MQ=9)
                                                                     aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac                                                                               <  2:181257/151‑1 (MQ=255)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                             <  1:10450/151‑1 (MQ=255)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                             <  1:234260/151‑1 (MQ=255)
                                                                         aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                           <  1:432542/151‑1 (MQ=255)
                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGGACTGCCCCCCCCCCCCCGGCCCCGTTCTTT‑CCCACGCTTCAATCCGCT‑CCCTGCGGGTTTTTTCTTCTCCGGGGAGCGGTGACCGACAAaaaacaaa                                                                     >  1:247696/1‑146 (MQ=255)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa                                                                  <  1:533593/150‑1 (MQ=255)
                                                                                             cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCGGTTCTTT‑GCAACATTCAAATCCCCC‑TCCCGCGGGTTTTTTCCACTCCGCAAAGTGGTCACCGACAAAAAACAGATAaaaaaaaaaac                                                        >  1:333243/1‑147 (MQ=9)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGGCGGTGCTTT‑GCCACGTTCCAATCCGCT‑CCCCGCGGGGTTTTTCTCCTCGGGGGGGGGTTTCCCGCCACACAACAGATAAAAAAAAAAGCGCAggc                                                  >  2:291938/1‑149 (MQ=9)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTCCTGCCCCCCCCCCCTCCGGCCCTTGCTTC‑GCAAACGTCAAATCCCCT‑CCCGGCGGATTTTTTCCACTCGGGAGAGGGTTCACCGACAAACAACAAATAAAAAAAAAAGCCCAGTc                                                  >  1:238718/1‑151 (MQ=255)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                              <  2:67518/151‑1 (MQ=39)
                                                                                                        tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                              <  2:589558/149‑1 (MQ=38)
                                                                                                         ttgatGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCCCCCCCCGGGGCGGTTCTTC‑GCCAACTTCAAATACCCC‑CCCCGCGGGGTTGTTCTACTCCGGGGGGGGGTTACCGACCAAAAACAGATAAAAAAAAAAGGCCAGCCTTcctcc                                           >  2:534130/1‑148 (MQ=255)
                                                                                                                         tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                            <  2:515605/150‑1 (MQ=38)
                                                                                                                          aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                           <  2:489646/150‑1 (MQ=38)
                                                                                                                           tGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                         >  2:255381/1‑151 (MQ=38)
                                                                                                                                           cGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc         <  2:66243/151‑1 (MQ=40)
                                                                                                                                                  cTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc  >  2:672817/1‑151 (MQ=40)
                                                                                                                                                   |                                                                                                                                                       
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GCAACGTTCAAATCCGCT‑CCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  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: ‑