Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913734,6820AG52.0% 2.4 / 24.5 29intergenic (+580/‑968)ybfB/ybfCuncharacterized protein YbfB/uncharacterized protein YbfC
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base G (5/8);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGG  >  NC_000913/734433‑734919
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
tACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGACACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCGGCACCCTGCCGGATAc                                                                                                                                                                                                                                               >  2:129492/1‑251 (MQ=16)
  cGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATAccc                                                                                                                                                                                                                                             >  1:153892/1‑251 (MQ=33)
   gccccGCTGGGCCGCGGGGTGGCAAAACGGGTGTGGCGACGTGAACGGTACCTGAGGGGCTGGAGGTCGCGGTCAGGGAAGCTGCAAGTGACCTGGTACGGCTGGGACGGCGACTGCTTGACCACGGTACAGAACGACAGAACCCGCGCCCGGACGATTGATCGGCCGGGGAGCTGCACGCCACTCAGCAGGGTTGAAACCGCTCCCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATAccc                                                                                                                                                                                                                                             <  1:172371/248‑1 (MQ=11)
                                         gcGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGGGTGgtg‑t                                                                                                                                                                                                     >  2:102257/2‑251 (MQ=9)
                                             gAACGGGACCTGCCGGGCTGGAGGTCGCTGTCACTGAAACCGCAAGGGACCTGGTACGGCTGGGACGGTGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTTGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTT‑TTccc                                                                                                                                                                                                 <  1:184641/251‑1 (MQ=11)
                                              aaCGGGCCCGGACGGGCTGGATGGCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCccg                                                                                                                                                                                                <  2:68130/251‑1 (MQ=9)
                                                           gggCTGGATTTCGCTGCCAGGGAACCCCCAAGTGACCGGTTCCGGCTGGGACGGCGCCCGCCTGCCCACGATACAGAACGACAGAACCCGCACCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCAGCCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGc                                                                                                                                                                                   <  1:76282/251‑1 (MQ=11)
                                                                   tGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTc                                                                                                                                                                           >  1:84706/1‑251 (MQ=14)
                                                                          gTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGCCGGCGACCGCCTGACCAAGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGACGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGAGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGAAGCCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGTGGAG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGc                                                                                                                                                                    >  1:106925/1‑251 (MQ=11)
                                                                                aaccacgaaaGTGACCGGTTACGGCTGGGACGGCGCCCGCCTGACCAAAATACAGAACGACAGACCCCGCATCCAGACGATTTATCAGCGGGGGAGCTTCACCCCCCTCACCAAGGTGGAACCCGCCCCCGGTGAGCTGGCGAAAACGCCGCGCCGCAGCCTGGCGGATGCCCTCCAGCAGTCCGGCGGCGAAGACGGGGGCAGTGGGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaaa                                                                                                                                                              <  2:152059/243‑1 (MQ=11)
                                                                                          ggagccgggtccggcggGGCCGGCGACCGCCGGACCACGAAAAGGAACGAAAGACCCCCCAGCCAGACGATTTATCAGCGGGGGAGCTTCACGCCACTCAGCAGGGTTGAAACTGCCCCGGGTGAGCGGGCGAAACCGCAGCGCCGCAGCCGGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGGGGCGGTGGGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCt                                                                                                                                                    <  2:70391/240‑1 (MQ=2)
                                                                                            gcccTGGTCCGGCGGGGAGGGCGCCCGCCTGCCCACGATACAGAACGACAGAACCCGCACCCAGACGATTTATCAGCCGGGGAGCTCCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                                                                                                                   <  1:78713/248‑1 (MQ=11)
                                                                                            gaccggggagggcggggACGGCGCCCGCCTGCCCACGATACAGAACGCCAGACCCCGAATCCATACAATTTATCAGCCGGGGAGCTTCACCCCACTCATCAGGTTTGAACCTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                                                                                                                  <  1:229176/237‑1 (MQ=11)
                                                                                              ccgggTACGGCTGGGACGCCGCCCGCCGGACCACGATACAGACCGACAGAACCCGCATCCAGCCGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATCCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCGGTGTGTTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                                                                                                                <  1:158861/248‑1 (MQ=11)
                                                                                               ggggTACGGCGGGGACGGCCACCGCCTGACCACAATACAAACCGACAGACCCCGCACCCAGACGATTAACCGGCGGGGGGGCTTCACGCCACTCATCGGGGTTGAACCGGCCCCCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGAGGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTg                                                                                                                                               <  2:242687/249‑1 (MQ=11)
                                                                                                tgttACGGCGGGGCCGGGGCCCGCCTGCCCACAATACAGAACAACAGACCCCGCGTCCATACGTTTTATCACCCGGGAAGTTTCACGCCACTCATCAGGGTTGAAACCGCCTCCGTTGAGCTGGGGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACTGTGTCAGGGTGGTG‑TTCCCTCCGGTGCTTGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                                                                                                                <  1:158859/246‑1 (MQ=11)
                                                                                                       gcggggAGGGCGCCCCCCTGACCACGATACAGAACGAAAGAACCCGCTTCAAGACGTTTTACCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCTGGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtg                                                                                                                                       <  1:222657/248‑1 (MQ=11)
                                                                                                                      gCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAACGGGTGAGTGAGGAAAgccgc                                                                                                                        >  2:205561/1‑251 (MQ=14)
                                                                                                                       ccTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgccgcc                                                                                                                       >  1:16703/1‑251 (MQ=31)
                                                                                                                             cACGATACAAAACAACAGAACCCGCATCAAGACGATTTATCAGCGGGGGAGCTCCCCGCACCTCATCAGGTTTGAACCCGCCACCGGTGAGCTGGCGAAACCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                                                                                                                 <  1:60159/251‑1 (MQ=11)
                                                                                                                                       aaCGACAGAACCCGCATACAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGATGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg                                                                                                       >  2:130200/1‑251 (MQ=9)
                                                                                                                                            cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc                                                                                                  >  2:179131/1‑251 (MQ=11)
                                                                                                                                                           gccgttttACCAGCGGGGGAGCTCCCCGCCACTCATCGGGTTTGAACCCGCCACCGGTGAGCTGGCAAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGCCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGa                                                                                   <  1:70927/246‑1 (MQ=11)
                                                                                                                                                                          ggggAGCTCCACGCCACTCATCAGGGTGGAACCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATgg                                                                    <  1:269281/251‑1 (MQ=21)
                                                                                                                                                                                      cccACTCATCAGGGTTGAACCCGCCACCGGTGAGCTGGCGAAACCGCAGCGCCGCAGCCGGGCGGATCCCCTTCAGCAGCCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTaca                                                        <  1:61970/250‑1 (MQ=255)
                                                                                                                                                                                                                               aaaCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATcgc               >  1:158576/1‑251 (MQ=33)
                                                                                                                                                                                                                                   gcggccccgcAGCTTGGCGGATCCCCTTCGGCAGCCCGGGGGCGAAGCCGGTGGCTGT‑TGGTGGTCCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGcc           <  1:9341/245‑1 (MQ=11)
                                                                                                                                                                                                                                          cgcAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCCGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATg                                                                     >  1:261720/1‑186 (MQ=21)
                                                                                                                                                                                                                                            cAGCCTGGCGGCTACCCTTCAGCAGTCCGGCGGCGCAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAACGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGCGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGCAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTgg  >  2:172563/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
TACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGG  >  NC_000913/734433‑734919

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: ‑