Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913734,6360CT38.5% 11.4 / 18.9 27intergenic (+534/‑1014)ybfB/ybfCuncharacterized protein YbfB/uncharacterized protein YbfC
Reads supporting (aligned to +/- strand):  ref base C (5/11);  new base T (5/5);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACA  >  NC_000913/734418‑734865
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     
gAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGcgccg                                                                                                                                                                                                       >  1:195015/1‑251 (MQ=34)
       gCTATCTTTACGCCCCGCTGGGCCGCGGGGTGGCAAAACGGGTGTGGCGCCGGGACCGGGCCCTGACGGGCTGGATGTCGCTGTCCCGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACTACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCAGCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTg                                                                                                                                                                                                <  2:150213/251‑1 (MQ=11)
               tACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGACACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCGGCACCCTGCCGGATAc                                                                                                                                                                                        >  2:129492/1‑251 (MQ=16)
                 cGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATAccc                                                                                                                                                                                      >  1:153892/1‑251 (MQ=33)
                  gccccGCTGGGCCGCGGGGTGGCAAAACGGGTGTGGCGACGTGAACGGTACCTGAGGGGCTGGAGGTCGCGGTCAGGGAAGCTGCAAGTGACCTGGTACGGCTGGGACGGCGACTGCTTGACCACGGTACAGAACGACAGAACCCGCGCCCGGACGATTGATCGGCCGGGGAGCTGCACGCCACTCAGCAGGGTTGAAACCGCTCCCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATAccc                                                                                                                                                                                      <  1:172371/248‑1 (MQ=11)
                                                        gcGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGgtgt                                                                                                                                               >  2:102257/2‑251 (MQ=9)
                                                            gAACGGGACCTGCCGGGCTGGAGGTCGCTGTCACTGAAACCGCAAGGGACCTGGTACGGCTGGGACGGTGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTTGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTTTTccc                                                                                                                                           <  1:184641/251‑1 (MQ=11)
                                                             aaCGGGCCCGGACGGGCTGGATGGCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccg                                                                                                                                          <  2:68130/251‑1 (MQ=9)
                                                                          gggCTGGATTTCGCTGCCAGGGAACCCCCAAGTGACCGGTTCCGGCTGGGACGGCGCCCGCCTGCCCACGATACAGAACGACAGAACCCGCACCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCAGCCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGc                                                                                                                             <  1:76282/251‑1 (MQ=11)
                                                                                  tGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTc                                                                                                                     >  1:84706/1‑251 (MQ=14)
                                                                                         gTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGCCGGCGACCGCCTGACCAAGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGACGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGAGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGAAGCCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGTGGAGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGc                                                                                                              >  1:106925/1‑251 (MQ=11)
                                                                                               aaccacgaaaGTGACCGGTTACGGCTGGGACGGCGCCCGCCTGACCAAAATACAGAACGACAGACCCCGCATCCAGACGATTTATCAGCGGGGGAGCTTCACCCCCCTCACCAAGGTGGAACCCGCCCCCGGTGAGCTGGCGAAAACGCCGCGCCGCAGCCTGGCGGATGCCCTCCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaaa                                                                                                        <  2:152059/243‑1 (MQ=11)
                                                                                                         ggagccgggtccggcggGGCCGGCGACCGCCGGACCACGAAAAGGAACGAAAGACCCCCCAGCCAGACGATTTATCAGCGGGGGAGCTTCACGCCACTCAGCAGGGTTGAAACTGCCCCGGGTGAGCGGGCGAAACCGCAGCGCCGCAGCCGGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCGGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCt                                                                                              <  2:70391/240‑1 (MQ=2)
                                                                                                           gcccTGGTCCGGCGGGGAGGGCGCCCGCCTGCCCACGATACAGAACGACAGAACCCGCACCCAGACGATTTATCAGCCGGGGAGCTCCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                                                             <  1:78713/248‑1 (MQ=11)
                                                                                                           gaccggggagggcggggACGGCGCCCGCCTGCCCACGATACAGAACGCCAGACCCCGAATCCATACAATTTATCAGCCGGGGAGCTTCACCCCACTCATCAGGTTTGAACCTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                                                            <  1:229176/237‑1 (MQ=11)
                                                                                                             ccgggTACGGCTGGGACGCCGCCCGCCGGACCACGATACAGACCGACAGAACCCGCATCCAGCCGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCGGTGTGTTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                                                          <  1:158861/248‑1 (MQ=11)
                                                                                                              ggggTACGGCGGGGACGGCCACCGCCTGACCACAATACAAACCGACAGACCCCGCACCCAGACGATTAACCGGCGGGGGGGCTTCACGCCACTCATCGGGGTTGAACCGGCCCCCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGAGGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTg                                                                                         <  2:242687/249‑1 (MQ=11)
                                                                                                               tgttACGGCGGGGCCGGGGCCCGCCTGCCCACAATACAGAACAACAGACCCCGCGTCCATACGTTTTATCACCCGGGAAGTTTCACGCCACTCATCAGGGTTGAAACCGCCTCCGTTGAGCTGGGGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACTGTGTCAGGGTGGTGTTCCCTCCGGTGCTTGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                                                          <  1:158859/246‑1 (MQ=11)
                                                                                                                      gcggggAGGGCGCCCCCCTGACCACGATACAGAACGAAAGAACCCGCTTCAAGACGTTTTACCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCTGGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtg                                                                                 <  1:222657/248‑1 (MQ=11)
                                                                                                                                     gCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAACGGGTGAGTGAGGAAAgccgc                                                                  >  2:205561/1‑251 (MQ=14)
                                                                                                                                      ccTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgccgcc                                                                 >  1:16703/1‑251 (MQ=31)
                                                                                                                                            cACGATACAAAACAACAGAACCCGCATCAAGACGATTTATCAGCGGGGGAGCTCCCCGCACCTCATCAGGTTTGAACCCGCCACCGGTGAGCTGGCGAAACCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                                                           <  1:60159/251‑1 (MQ=11)
                                                                                                                                                      aaCGACAGAACCCGCATACAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGATGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg                                                 >  2:130200/1‑251 (MQ=9)
                                                                                                                                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc                                            >  2:179131/1‑251 (MQ=11)
                                                                                                                                                                          gccgttttACCAGCGGGGGAGCTCCCCGCCACTCATCGGGTTTGAACCCGCCACCGGTGAGCTGGCAAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGCCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGa                             <  1:70927/246‑1 (MQ=11)
                                                                                                                                                                                         ggggAGCTCCACGCCACTCATCAGGGTGGAACCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATgg              <  1:269281/251‑1 (MQ=21)
                                                                                                                                                                                                     cccACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGAAACCGCAGCGCCGCAGCCGGGCGGATCCCCTTCAGCAGCCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTaca  <  1:61970/250‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     
GAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACA  >  NC_000913/734418‑734865

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