Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA42.3% 34.2 / 11.1 26L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (12/3);  new base A (0/11);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357742‑358035
                                                                                                                                                      |                                                                                                                                                     
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                      <  2:1029326/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCTGCCGCAGTCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                     <  1:967669/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                     >  1:908025/1‑150 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                 >  2:431597/1‑151 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                 <  2:1029039/151‑1 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                              >  1:956022/1‑151 (MQ=255)
       ctAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                  <  1:793542/146‑1 (MQ=255)
        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                             >  1:785286/1‑151 (MQ=255)
                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                     >  2:790741/1‑150 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                           >  1:134257/1‑151 (MQ=255)
                                  gtccTGTGTCCGCGGCTCCCGTGATTCTTTCCCCCAAACCACCGCTATTTTACCCAA‑CGCGATGGGCAAAAAAGCGTT‑GGGGGCGGTTGTCTGCAACTTATTTTCCGGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑Aag                                                                                                                   <  2:117488/148‑1 (MQ=255)
                                    ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                        >  2:733225/1‑144 (MQ=255)
                                      tgtgTCAGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCAC‑CGCCTTGTGAAAAAATGCCTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAgag‑a                                                                                                               <  2:261310/149‑1 (MQ=255)
                                        tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGGAGTGGCggggggag                                                                                                                    >  1:68381/1‑140 (MQ=255)
                                        tgacggCGGCTCCGGTGTTTCTTTCCCGCCAACACCCGCATTTTTCCCCAC‑CCCCATGTGAAAAATTGCTTC‑CGGGGCGGCTGTCTGCCGGTTTTTTTCCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                          <  1:9844/133‑1 (MQ=255)
                                               ggCTCCCGTATTTCATTCCCGCAAACACCCGCGCATTTACCCAC‑CGCCTTGTGAAAAAATGCCTC‑CGCGGCGGTTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCtt                                                                                                       <  2:999392/148‑1 (MQ=255)
                                                             aTTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAAGAG‑AAAAAATTTAAAGccgtgactcttt                                                                                       >  1:788053/1‑141 (MQ=25)
                                                                                  tttACCCAA‑CGCAATGAAAAAAAACCCCTC‑CGGGGCGTCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGt                                                                  <  1:255967/148‑1 (MQ=255)
                                                                                   ttttCCAACCGCTATGTAAAAAAACGCTTT‑CGGGAGGGCTTTTTGCAAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGct‑gct‑g                                                                    <  2:839742/143‑1 (MQ=255)
                                                                                       ccAC‑CGCCATGGAAAAAATTGCCTC‑CGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACa                                                              <  2:372900/147‑1 (MQ=255)
                                                                                                  gtaaaaaaaTCCCTC‑GGCGGCGCCTTTTTGCAATCTTTTTGCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAg                                                    <  1:751062/143‑1 (MQ=255)
                                                                                                    taaaaaaTGGCTTCCGGGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGCGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGaa                                                <  2:664413/148‑1 (MQ=255)
                                                                                                            gCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑AAGAG‑AAAAAAATTTTAGCCCCCCATCTTTTGTGGGGAACCCT‑GCTGGGGAACAATGGGGGAAggaagggggggg                                         >  1:255091/1‑139 (MQ=255)
                                                                                                             ccTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑AAGAG‑AAAAAGATTTAAAACCCTCATTCTTTGTGGCTATCGCTCGCT‑GTAAAGAATTGGGGAAGGGagggcgggcg                                        >  1:719707/1‑140 (MQ=255)
                                                                                                                     gggggcttccTGCCATTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                                <  2:24668/140‑1 (MQ=255)
                                                                                                                        ggcTGTCTGCCATCTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                             <  2:956036/149‑1 (MQ=255)
                                                                                                                                   aGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑GAGAGAAAAAAAATTTAAACCC‑TCCTGTTTTTTGGGAATAGCT‑GCT‑GGGGACAAAATGGAGAGGGGAGGGCGGCGAAAAAGCCCCCGAACCCCGCCaa                 >  2:26525/1‑150 (MQ=255)
                                                                                                                                    gtggtttgcccgcagAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGGTTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                <  1:27151/141‑1 (MQ=255)
                                                                                                                                          tttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑GAGGG‑AAAAAACCTTATGCCCCTCATGTTGTGTTGGAAACGCT‑GCT‑GGTAACAATGTGGGAAGGGGAGGGGGGCGGAACAGCTCCCCCACCGCGCCaaatttg            >  2:86141/1‑143 (MQ=255)
                                                                                                                                                  agaacaaAACCGCCGTCCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAGCGCT‑GCT‑GGTAACAATGTGGAAAGTGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCGTCaa  <  1:192266/148‑1 (MQ=255)
                                                                                                                                                    aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa   <  2:27866/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357742‑358035

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