Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2780GC43.8% 34.3 / 11.7 32intergenic (‑1008/‑704)yohF/yohPputative oxidoreductase YohF/putative membrane protein YohP
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base C (5/9);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCC  >  NC_000913/2228029‑2228523
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                
cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGc                                                                                                                                                                                                                                                                 <  1:580701/251‑1 (MQ=255)
 gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGcc                                                                                                                                                                                                                                                                <  1:560598/251‑1 (MQ=255)
                 aaTATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAAGGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAg                                                                                                                                                                                                                                              <  2:623138/250‑1 (MQ=255)
                                acCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTTGCGGGGTTTCTCGCGGTTCGCCCGCCGTGCTTTTAACCGACCCCCGCGGGCAGGAACCTCCGCCCGAGTCTTTTCCCCGataaccgaaagccaac                                                                                                                                                                                                                              >  1:359677/1‑238 (MQ=255)
                                    tCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGTTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                                                                                                                                         <  1:570700/249‑1 (MQ=255)
                                       acagaagaaGTAATCCCTTTTTTGTATAGTTGGCGATAGGGCGCTATAATGGGGTTCTGTTCTTTTGCTATTTGCTGCAAGACCGTGTTTAACGCCGCCGTGGTTTGCCACTCCCAGTCCCGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGGGCTTTAACCGTACCCAGCCGGGCAGTAACCGCCGCCCGATTTTTTCCCCAGATATCGAGAAGCCATCccgcc                                                                                                                                                                                                                      <  2:429291/246‑1 (MQ=255)
                                                                ataGCTGGCGATGGGTCGCGATAATGGTTTTCTGCTCTTTTTTTATTTGCGGCAAGACCGTGTTAACCGCCGCCTGGGTTTGCCACCCCAAGTACAGGGGGGCTACGCGGCCAGCCAGAAATGGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGAACCCGGGCGGGAAGTAACCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                                                                                                                         <  1:411284/250‑1 (MQ=255)
                                                                 taGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCCCCCCCCCGGCGTTTAACGTACCCAAGCGGGCCGGAACCTCCGCCCGGTTCTTTTCCCCGAAAACGAGATTGCCGCCCGCCGCTAAACCaaaaacacccttgggggaac                                                                                                                                                                                        >  1:437293/1‑235 (MQ=255)
                                                                      ggCGACCGGTCGCTATAATGGTTTTCTGCTCTTTTTCTATTTGCTGAAGGACCGTGTTTAACGCCGCCTGGGTTTGCCCCTCCAAGTACGGGGGGGTTCCGCTGCCAGCCAGCAATGGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGTTTTAACCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACggg                                                                                                                                                                                 <  1:467327/251‑1 (MQ=255)
                                                                            cggtggtgttagtggtttcttGCTTTTTTTCTATTTGCTGCAGGACGGTGTTTAACGCCGCCTGGGTTTCCCCCCCCAAGTACGGGGGGGTTACGCTGCCAGCCAGAAATGGGGGGGTTGGCTCGGGTCCCGCCCCCCGTGCTTTACCCGTCCCCAGGCGGGAAGTAACCCCCGCCCGTTTCTTCCCCCAGATATCGATATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTc                                                                                                                                                                           <  1:682271/231‑1 (MQ=255)
                                                                                ggggtttagggtgttctctttttttttTATTTGCGGCAAGACGTGGTTAACCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGCTTACGCGCCCAGCCAGAAAGGGGCGGGTTTGCTCGCGTCCCGCCGCCCGGTTTTTAACCGTCCCCAGGCGGGCAAAACCCTCCCCCCGATTCTTTCCCCAGATACCGAGATGCCAGACCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGTGCCGGTCGTAc                                                                                                                                                                       <  1:167505/230‑1 (MQ=255)
                                                                                                  ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCGCCCGGGCTTTTACCCTAACCCAGCGGGCCGTAACCTCCCCCCGATTCTTTTCCCCGATATCGGGATGCCAGCCCCCCCTTAAGCCAAAAATAACGTTTGTGTACCAccggccgggcggaacctcggcccggtcgttc                                                                                                                                                      >  1:90245/1‑227 (MQ=255)
                                                                                                            aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTCCCCCCCCCCGGCTTTAAACCCACCCAGGCGGGGCGGAAACTCCCCCCCCTTCTTTTCCCCGATATCCGGAAGCCCACCCCCCCCTTAAACCAAAAGAACCGTGGGGGACCACGGGCCGGGCCTACCTGCGGGCGGGTCGtttcaagcaaaa                                                                                                                                             >  1:465123/1‑239 (MQ=255)
                                                                                                                 ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTTCGCCCCCCCTGCTTTTACCCTACCCAGGCGGGCAGTAAACTCCCCCCCGGTTTTTTCCCCCATATCGAGATGCCCGCCCCCCGTTAAGCCAAAATAACGGTTGTGTACCCCGGGCCCGGCGTACCCTCGGCCGGGTCGTTTCGAGCAAAACGGccccc                                                                                                                                     >  1:140211/1‑250 (MQ=255)
                                                                                                                               tttcccgcgccgtgtttttGCCACTCCGAGTACTGGCGGCCCGCGCTGCCCGCCAGCAGGTGGCGTGTTTCCCCGCGTTCCGCCGCCCGTTTTTTACCCGTCCCCGGGGGGGCAGAACCCCCCCCCCGTTTCTTCCCCAAGTGAGCGAGCTCCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                                                                                                                         <  2:388356/234‑1 (MQ=255)
                                                                                                                                  aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCCCGCGGTCCGCCGCCCCTGCTTTTACCGTACCCCAGCGGGGAGGAAACCCCCCCCCGTTTTTTTCCCCAATATCGGGATGCCCGCCCCCCGTTAAACCCAAAAAACCCGTTGGGGACCCCCGGGCCGGCGTGCCTGCGGCCGGGGCGGTCCGAGCAAAACGCCCCCTTTAGCCCTTCTGccac                                                                                                                     >  1:549064/1‑249 (MQ=255)
                                                                                                                                         ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAAATGGCGGGTTTGCTCGCGGTTCGCCGCCCGGGCTTTTACCCTCACCAGGCGGGGAGGAAACCCCCCCCGGTTCTTTTCCCCGATAACGAGATGCCCGCCCCCCGGTTAACCAAAAAGACCGGTGGTGGACCCAAGGGCGGTCGGACCTCCGGCCGGGTCGTTCAGAGGAAAAGGCCCCCCTAAACCCTCTGCCGACAttttt                                                                                                              >  1:342712/1‑251 (MQ=255)
                                                                                                                                            gggttggccccccccAGTACGGGCGGGCTACGCGCCCAGCCAGCAATTGGCGGTTTTGCTCGCGTCCCGCCCCCCGTTCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAAATATCGAAATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCaga                                                                                                                                                   <  1:424360/200‑1 (MQ=255)
                                                                                                                                                 tccctcccccggtacgggggggCTACGCTGCCAGCCAGCAAGTGGCGGGTTTGCGCGCGTCCCGCCGCCCGTGTTTAACCCGTCCCCAGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGAGGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCt                                                                                                       <  1:10942/231‑1 (MQ=255)
                                                                                                                                                  cccactccattgccggggtggCTACGCGGCCGCCCGGCAGTTGGTGGTTTTGCTCGCGTGCCCCCGCCCGTGCTTCAACCTTCCCCGGGGGGGCAGTACCCCCCGCCCGATTCGTTCCCCAAACTTCGAGTTCCTCCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCTCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTcc                                                                                                     <  1:159590/232‑1 (MQ=255)
                                                                                                                                                   ccccccccAGTCCAGGCGGGCTCCCCGCCCCGCCAGCAATGGGTGGTTTTCCTCCCGTCCCGCCGCCCGTGTTTTACCCGTACCCAGGGGGGCAGTACCCCCCCCCCTTTCTTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCa                                                                                                    <  1:129016/246‑1 (MQ=255)
                                                                                                                                                                                     agcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAAACGTACCCAGGCGGGCAGTAACCCCCGCCCGCTTCTTTCCCCAGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTCACCTTGGTGTAACACGGGCCGGGCGGACCTGCGGCCGGATCGTTCAGAGCAAAAGGCCCCCTTAAGCCTTCTGCCGACATTTTTTTCCGCTCCCTATTCGCCGAACAGTCGGTCTGCGGACCATCCtt                                                                  >  1:130654/1‑250 (MQ=255)
                                                                                                                                                                                           gggCGGGTTTCCCCGCGTTCCCCCGCCGGTGCTTTACCCGCCCCCAGGGGGGCATTACCCCCCCCCCGATTCTTTCCCAAGATATGGAGTAGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGc                                                            <  1:200232/250‑1 (MQ=255)
                                                                                                                                                                                                ggTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGACCGGTCAGAGCAAACGGCCCCATTAAGCCTTCCGCCGACATTTTTTGCCGCTCCATAGCCGCCGAACAGTCGATCGGCGGACCATCCTGAGTGGCAACTg                                                       >  2:716151/1‑251 (MQ=255)
                                                                                                                                                                                                        cgcgTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGcc                                               >  1:308876/1‑251 (MQ=255)
                                                                                                                                                                                                        cgcgTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGAACATCCTGAGTGGCAACTGCCTTCGcc                                               >  1:480201/1‑251 (MQ=255)
                                                                                                                                                                                                           gTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTgg                                            >  1:622384/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccccccGTGCTTTCCCCGCACCCAGGCGGGTGGTACCCCCCCCCGATTTCTTCCCCCAGATATCGGGATCCAGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCa                                      <  1:16476/248‑1 (MQ=255)
                                                                                                                                                                                                                                       tACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTTCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGt                    >  2:497417/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ccaccccgggcggcaagaccccaccgcccgcTTCTTCCCCCATATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTTCCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCt                      <  2:513643/228‑1 (MQ=255)
                                                                                                                                                                                                                                          cccAGGGGGTCAGTACCCCCCCCCCGATTCTTCCCCCAAATACCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTc                  <  1:74399/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCat      >  1:684056/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           ccTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATcc  >  1:479916/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                
CGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCC  >  NC_000913/2228029‑2228523

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