Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA42.9% 39.2 / 23.4 42L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/7);  new base A (0/18);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 2.26e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357745‑358028
                                                                                                                                                    |                                                                                                                                            
atgtgtcgggcATCAGGCAATTCCAGCCGCGGACATGTGTCAGCGGCTACCGTGATTCATTTCCTCCAACAACCGCGTCGTCATCCAACCCCCTGTCCCATAATGCCCT‑CTCGGCGGCT‑GTCTTCCAGCTATTTTT‑CCGCCGCAACaaac                                                                                                                                          <  2:97376/141‑2 (MQ=17)
aaGCGTAGCGCATCAGGCAATTCCAGCCGAAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCACCAAAAc                                                                                                                                         >  2:808749/1‑151 (MQ=255)
     tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                     >  2:8948/1‑150 (MQ=255)
              aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                            >  2:336118/1‑150 (MQ=255)
               ggaaaccccACCCGCAGCCCGGTGTCCGCGGCTCCGGTGATTAATTCCCGCCAGCACCCGCGTATTAACCCAGCGCCA‑GGGCAAAAAACCCCT‑CGGGGGGGGTGGTCTCCCGGCTTTTTTC‑CCCCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                           <  2:9119/143‑1 (MQ=255)
               ggCAATTCCAGCCGCAGCCTTGTGTCGGCGGCTACCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAATTGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCag                                                                                                                          <  1:927432/151‑1 (MQ=255)
                 cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                         >  1:422261/1‑150 (MQ=255)
                    ttCCAGCCCAATCCCTTTGCCAGCGGCTCCGGTGTTTCTTTTCCGCCAACACCCCCGCTTTTTTCCAACGCCATGTACAAAATTCCTTG‑CGGGGGGGCT‑GTCTGCCAGTTTTTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                      <  1:848376/150‑1 (MQ=255)
                    tcccACCCGCAGCCTGGTGTCAGCGGCTACGGTGATTCTTTCCCCCCACCAACCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTT‑CGGGGGGGCT‑GTCTGCCAGCTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                      <  1:944520/148‑1 (MQ=18)
                     tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGGGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    >  2:35156/1‑151 (MQ=255)
                        aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                 >  1:362075/1‑151 (MQ=255)
                          ccGCAGCCCTGTGTGAGGGGCTCGTGGTTTTCTTTCCCGCAACCACCCGCGCTTTAACCCACCCCCTTGTGC‑AAAATCGCTT‑CGGGGGCGGT‑GGCCGGCAGTTTTTTTTGCCGCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGc                                                                                                                   <  2:270161/147‑1 (MQ=255)
                          ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTATGGCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                       >  1:377940/1‑143 (MQ=255)
                           cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                 >  2:670049/1‑148 (MQ=255)
                                   tgtttaggCGGCTCCCGTGATCCTTTCCCGCCAACACCCGCGTTTTTACCCACCCCCTTGTGAAAAAATCCTTT‑CGGGGCGGCT‑GTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                          <  1:60502/140‑1 (MQ=255)
                                      gtaaGCGGCTCCTGTGTTTCATTCCCCCCAACACCCGGTTATTTAACAAACCCCATGTGAAAAATTGCTTC‑CGGGGCGGCT‑CTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt                                                                                                    <  1:927840/146‑1 (MQ=255)
                                            ggCTCCCGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCCACCCCCGTGTCAAAAAATCCCTT‑CGCGGCGGCT‑GTCTGCCGGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                                     <  1:1047150/142‑1 (MQ=255)
                                            ggCTACCGTGATTCTTTCCCCCAAACACCCGCCTTTTTACCAACCGCCATGTG‑AAAAATTCCTTCCGGGGCGGCT‑GTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTaa                                                                                           <  1:463497/151‑1 (MQ=255)
                                              cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaagt                                                                                               >  1:240837/1‑140 (MQ=255)
                                               tACCGTGATTCATTTCCGCCAGCAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGACGAG‑AGAAAG‑CGTTTAgac                                                                                        >  2:205751/1‑149 (MQ=21)
                                                 ccGTGATTCATTCCCCCCAACACCCGCGCATTTATCCACCCCCTTGTAAAAAAATGCCTC‑CGGGGCGGCT‑GTCTGCCAGTTATTTTC‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGc                                                                                         <  1:695044/148‑1 (MQ=255)
                                                             tCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                  >  2:808750/1‑113 (MQ=255)
                                                              ccccaaacaacCGCCCTTTAACCCACCGCCATGT‑CAAAAATGGCTTCCGCGCCGGC‑‑GTCTGCCAGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTtgtg                                                                            <  2:707340/143‑1 (MQ=255)
                                                              ccGCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                  <  1:808796/112‑1 (MQ=255)
                                                              ccGCCAACACCCGCCCTTTTACCAACCGCTATGTGAAAAAGTGCCTC‑CGCGGCGGCT‑GTCTGCCGGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTgg                                                                         <  2:62156/151‑1 (MQ=255)
                                                                    acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑GAAAAA‑CATTAAACCCCTcc                                                                                   >  2:230993/1‑134 (MQ=25)
                                                                                  acccAACACCGTGTACAAAATTGCCTC‑CGCCGCGTCT‑GTTTGCCATCTTTTTTC‑CCGCCGCAAAAAAACCGCCTTTCTCTCCAGTAGTGGCGGGGCAAGAT‑AAATAG‑CTTTGAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                        <  2:760515/146‑1 (MQ=255)
                                                                                     cacccccATGTAAAAAAATCCCTC‑GGGGGGGGCT‑GTCTGCCATTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                  <  2:72720/146‑1 (MQ=255)
                                                                                       aCGCCATTTACAAAAGTGCCTT‑CGCGGGGGGT‑GTCTGCGTGTTTTTTTC‑CCCCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                   <  2:1059963/148‑1 (MQ=255)
                                                                                        cGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAA‑CTTTAAACCC‑CCCTGTTTTGTGGGAAACGCtggtggggaaaa                                                      >  1:611827/1‑138 (MQ=25)
                                                                                           cATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAGAATTTAAACCC‑CCCTTTTGTGTGGGGATCGCTGCTGGGGACAAAGTgg                                                 >  1:353179/1‑146 (MQ=25)
                                                                                             tgggaaaaatTCCCTC‑CGGGGCGGCT‑GTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                           <  1:1020648/145‑1 (MQ=255)
                                                                                               tgCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAAAGAAAAAAC‑CTTTAAACCC‑CCATGGTGTGTGGCAATCGCTCCTggggaaaaagtgggaa                                              >  1:525807/1‑131 (MQ=255)
                                                                                                gCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑AAAAAG‑CATTAAAAAGGGCGTGCTGTTTGGCCATCGCCTCTGGGGACAATGGGGGAAgg                                            >  1:18427/1‑146 (MQ=255)
                                                                                                  aaaaaTCCCTC‑CGCGGCGGCT‑TTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                       <  1:336143/149‑1 (MQ=255)
                                                                                                                              ccAGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                     <  2:654314/141‑1 (MQ=255)
                                                                                                                                    aTTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAG‑AAAAAA‑AATTTATACCCCCATTTTTTGTGGGAATAGCTGGTGGGGACAATAGGGAAAGGGGAGTGGGGGGGATCAACCCCCGAAACccgcccatt        >  1:929622/1‑143 (MQ=255)
                                                                                                                                     ttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg     <  2:1013718/150‑1 (MQ=255)
                                                                                                                                     ttttc‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt        <  1:180035/147‑1 (MQ=255)
                                                                                                                                     ttttc‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt         <  1:912511/146‑1 (MQ=255)
                                                                                                                                        ttgccgcggCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt   <  2:715880/147‑1 (MQ=255)
                                                                                                                                           ccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  1:543420/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  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 < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑