Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA50.0% 22.9 / 22.5 40L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (16/4);  new base A (0/20);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.54e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357740‑358026
                                                                                                                                                        |                                                                                                                                          
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                            >  1:133534/1‑151 (MQ=255)
        cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                    >  1:375463/1‑151 (MQ=255)
           aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                 >  1:535071/1‑151 (MQ=255)
             cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                               >  1:349805/1‑151 (MQ=255)
              gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                              >  2:541577/1‑151 (MQ=255)
                       attcccAGCCGCAGCCTTGTGTCAGCGGCTACCGTGATTCTTTCCCCCAAACACCCGC‑GCTTTTATCCACCGCCTTGTGAAAAAATCCTTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  2:763287/146‑1 (MQ=18)
                         tcccAGCCGAAGCCCTGTGTCAGCGCCTACCGTAATTCATTCCCCCCAACACCCGC‑GCTTTTATCCAACCCCTTGTGAAAAAATGCCTTC‑GGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    <  2:163983/148‑1 (MQ=255)
                          tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                  >  1:45126/1‑151 (MQ=255)
                           ccAGCCGCAGACTTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGC‑GCTTTTACCCAACGCCATGTGAAAAAATGCCTCC‑GGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                    <  2:296652/148‑1 (MQ=255)
                               ccGGAGCCCTGTGTCTGGGGCTCCCGTGATTCATTCCCGCCAACACCCGC‑GCATTTACCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                               <  2:268957/149‑1 (MQ=18)
                               ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                             >  2:421789/1‑151 (MQ=255)
                                 gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                               >  1:715467/1‑147 (MQ=255)
                                    gccttgtgTTAGGGGCAACCGT‑ATTCCATTCCGCCAACTACCGCGGCTTTTACCCAACCCAATGTGAAAAAATCCTTTC‑GCGGCGGCTGTCGCCCAGCTATTTCACCCCCCAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  2:402837/144‑1 (MQ=255)
                                       cgtggccggCGACTCCCGTGATTCTTTCCCGCAAACACCCGC‑CCATTTACCCACCGCATTGTGCAAAAAGCCCTCG‑GCGGGGCTTGTCCGCCGGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                        <  1:42925/140‑1 (MQ=255)
                                                      ccttttttCGTTCCCGCCACCAACCGC‑GCATTTACCCAGCCCTATGTGCAAAAAT‑CCTTCGGCGGGGTCTTCCCGCTATCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAg                                                                                        <  2:188540/141‑1 (MQ=255)
                                                      ccGTTATTCATTCCCGCCAACACCCGC‑GATTTTATCCAACGCCATGTGAAAAAATGCCTTC‑GGGGCGGCTGTCTCCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCg                                                                                    <  2:316783/151‑1 (MQ=255)
                                                       cGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑G‑aaaaaaactttaaaccc                                                                                    >  1:526064/1‑136 (MQ=25)
                                                            ttCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGA‑A‑AAAAACCATTaacacc                                                                                     >  2:433699/1‑140 (MQ=25)
                                                                ttcccGCAACCACCCGC‑CCATTTACCCACCCCCATGGCAAAAAATGCCTCC‑GCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                           <  1:391667/147‑1 (MQ=255)
                                                                  tCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAAATTTTAACCCGCCAAGTtttgtg                                                                        >  2:650954/1‑151 (MQ=25)
                                                                      ccAACACCCGC‑GCTTTTACCCACCCCCATGTGAAAAAATCCCTTC‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                    <  1:65122/151‑1 (MQ=255)
                                                                        accacccGC‑GCTTTTCCCCACCGCTATGGGAAAAAATCCCTTC‑GGGGCGGCTTCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                  <  1:606031/146‑1 (MQ=255)
                                                                           accccc‑gcTTTTACCCAACGCCATTTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                               <  2:919631/146‑1 (MQ=255)
                                                                           aacCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCCGA‑G‑GAAAAGCTTTAAAAGCGGCCTGTTGTGTgggga                                                                    >  2:706203/1‑143 (MQ=25)
                                                                            acCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGG‑AAAAAGCTTAAACCCCgcacatttt                                                                            >  2:322970/1‑130 (MQ=25)
                                                                              cgc‑gCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAA‑G‑AAAAAAATTTTAACCCCCCACGTTTTGTGGGAATc                                                                  >  2:605593/1‑145 (MQ=25)
                                                                                              cGCCATTTGCAAAAATCCTTTC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                              <  2:28460/150‑1 (MQ=255)
                                                                                              cGCCATGTGAAAAAAGGCTTTC‑GCGGCGGCTGCCTGCCGACTTTTTCCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                              <  2:873466/150‑1 (MQ=255)
                                                                                                  aTGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑GAAAAAAATTTTAAACCCG‑CATGTTTTGGGGGAAACGCTGCCGGGGACaaaatgggaaagg                                         >  2:511680/1‑140 (MQ=255)
                                                                                                   tgtgAAAAAATGCCTGC‑GGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                          <  2:375506/149‑1 (MQ=255)
                                                                                                    gtgAAAAAATCCCTCG‑GGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                        <  1:942787/150‑1 (MQ=255)
                                                                                                     tgaaaaatTGCTTTC‑GCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                        <  2:749586/146‑1 (MQ=255)
                                                                                                     tgCAAAAATCCTTTC‑GCGGGGGCTGTCTGCCAGTTATTTTTCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                       <  1:779015/150‑1 (MQ=255)
                                                                                                        aaaaaTGCCTCC‑GCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                      <  2:535129/148‑1 (MQ=255)
                                                                                                                     gcggcgGCTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                         <  1:749937/149‑1 (MQ=255)
                                                                                                                                 tGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAACTTTTAAACCCCCATGTTTTGTGGGGATAGCTGCTGGTGAAAAAGTGGGAAGGGGAGGGGGGCGGATCAGCCCCAGAAACccgc           >  2:15777/1‑151 (MQ=255)
                                                                                                                                  cccAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  1:63488/150‑1 (MQ=255)
                                                                                                                                   ccattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                <  1:312226/137‑1 (MQ=255)
                                                                                                                                    gaggtttttccccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  1:419819/138‑1 (MQ=255)
                                                                                                                                          ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg   <  2:334764/150‑1 (MQ=255)
                                                                                                                                             ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCCAGG‑G‑AAAAAACATTAAACCCCCCCTGTTGTGTGGCGATAGCCGCTGGTGACAATTTGGGAAGGGGAGTGCGGGGGGTCAGCCCCACAACCCCGCCAAATGaat  >  1:502373/1‑146 (MQ=255)
                                                                                                                                                        |                                                                                                                                          
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357740‑358026

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