Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA52.4% 17.6 / 35.3 21L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (10/0);  new base A (0/11);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 2.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05

AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357745‑358028
                                                                                                                                                  |                                                                                                                                            
aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc                                                                                                                                         >  1:442555/1‑151 (MQ=255)
                                  cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag                                                                                                         >  1:665246/1‑149 (MQ=255)
                                     tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGCG‑aaaa                                                                                                    >  1:159741/1‑149 (MQ=255)
                                        cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGG‑aaaaa                                                                                                   >  1:48590/1‑146 (MQ=255)
                                         aGCGGCTACCGTGATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCAATGTGCAAAAATGCCTTCGGGGGGTCTGTCTGCCAGCTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                       <  1:486819/143‑1 (MQ=21)
                                               tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAG‑gaaaaaacttaaaac                                                                                         >  2:161111/1‑139 (MQ=25)
                                                 ccGTGTTTCATTCCCCCCAACACCCGCGCTTTTCCCCACCCCAAGGTGAAAAAATCCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                         <  1:127859/149‑1 (MQ=255)
                                                 ccGTGATTCATTCCCCCCACCACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGccc                                                                                        <  1:495704/150‑1 (MQ=37)
                                                   gtattccttcccccccAACACCCGCGCTTTTACCCAACCCCATTTGCAAAAA‑GCCCTCGGGGGGGGTTGCTGGCCGTTTTTTTTCCCGCCCCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg                                                                                       <  2:405949/135‑1 (MQ=255)
                                                   gTGATTCATTCCCGACACCACCCGCCTATTTACCCACCCCCTTGGGCAAAAAGTCTTGCGGGGCGGCTGTCTGCCATTTTTTTTC‑CCACCAAATAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCttt                                                                                              <  2:170574/142‑1 (MQ=255)
                                                     ttttCTTTCCCCCCACCACCCGCCCTTTTCCCCAACGCCTTGGGAAAAAATCCCTCCGCGGCGGCTGTCTCCGAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCa                                                                                    <  2:227747/148‑1 (MQ=255)
                                                     gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAa‑aaaa                                                                                                    >  1:65965/1‑133 (MQ=255)
                                                                                    ccacccccATGTGAAAAAATGCCTCGGCGGCGGCTGTCGCCAAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                    <  1:469587/145‑1 (MQ=255)
                                                                                               tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCC‑CCCTTTTGTGTGGGAATCGCTGCTGGtggaaatgt                                                    >  1:192424/1‑132 (MQ=255)
                                                                                                                         cTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAG‑AAAAAAATTTTAAACCGCCCTGTTGTGTTGGAAACGCTGCTGGGGAAAATAT‑GGGAAGGGGAGGGGGGGGGATCAGcccc                       >  1:120951/1‑142 (MQ=255)
                                                                                                                           gCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAA‑AAAAAAATTTAAAACCCCCACTTTTTGTGGGAAATGCTGCCGGGGACAAAGTGGGAAAGGGAAGGGGGGGGaacaaaacccacaaac                 >  1:646819/1‑133 (MQ=255)
                                                                                                                             cAGCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  1:452029/143‑1 (MQ=255)
                                                                                                                               ttttttttc‑CCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa          <  2:766793/145‑1 (MQ=255)
                                                                                                                                  tttttccccgccaCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:297687/143‑1 (MQ=255)
                                                                                                                                          cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAATTTTAAACCCCCCTGTTTTGTGGGAATAGCTGCCGGGGAAAAAGT‑GGGAAAGGGAGGGCGGGGGAACAAATCCCAAAACCCCCTAATTGAggt   >  1:362344/1‑144 (MQ=255)
                                                                                                                                            ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  2:777677/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357745‑358028

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