Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC50.0% 20.3 / 13.5 28K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (12/2);  new base C (0/14);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 5.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357735‑358025
                                                                                                                                                    |                                                                                                                                                    
tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTccgcc                                                                                                                                                 >  1:227066/1‑151 (MQ=255)
    ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑Gcaa                                                                                                                                            <  2:521094/151‑1 (MQ=255)
          aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAAc                                                                                                                                      >  1:701323/1‑151 (MQ=255)
                    cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTtctc                                                                                                                            >  2:579320/1‑151 (MQ=255)
                      tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTtctctc                                                                                                                          >  1:72053/1‑151 (MQ=255)
                       cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTcc                                                                                                                         >  1:680592/1‑151 (MQ=255)
                                    ccGCTGCCCGGTGTCGGCGTCTCCCGTGTTTCTTGCCCGCCAACACCCGCTCTTTTCCCCAC‑CGCAA‑TGGACAAAAAGGCCTCCGCGGCGCTTGTCGGCGA‑TCTTTTTTCCCCCC‑AAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            <  2:5075/151‑1 (MQ=255)
                                       cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTCTTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGAgggg                                                                                                            >  2:431060/1‑148 (MQ=255)
                                         gcctttggtCAGCGGCTACGGGTTTTTATTCCCGCAAACACCCGC‑CCTTTCCCCAACCGCAG‑TGTGAAAAAATGCCCCGGGGGGGGCTGTCGGCGA‑GATTTTTGCCCGCC‑GAAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                              <  2:317348/137‑1 (MQ=255)
                                          aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  2:288079/1‑145 (MQ=255)
                                                      ggCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGggaaaggaaaaaaatt                                                                                             >  1:658782/1‑142 (MQ=255)
                                                           ccGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCCAC‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCG‑GTTTTTTTCCCGCC‑GAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGc                                                                                       <  1:680620/148‑1 (MQ=25)
                                                              tGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattttaaacc                                                                                      >  1:607957/1‑134 (MQ=25)
                                                               tattcttccccccaaacaccCGCTCTTTTCCCCAA‑CGCCA‑TGTGCAAAAGCCCCTCCGGGGCGGCGGTCTGCCA‑GCTATTGTCCCGCG‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCa                                                                                 <  2:572864/136‑1 (MQ=255)
                                                                        ccGCCAACACCCGCCCATTTACCCAA‑CGCAA‑TGTAAAAAAAGGCCTCCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCC‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTtgtgtg                                                                        <  1:124819/150‑1 (MQ=255)
                                                                          gCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAAAAA‑GTTTTAGGCCGCCCTGTTGTgggggg                                                                     >  1:418748/1‑149 (MQ=25)
                                                                           ccAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAA‑ATTTaaagccccccagt                                                                              >  1:737141/1‑130 (MQ=25)
                                                                           caaccacccGCTCATTTCCCCAC‑CGCCA‑TGAAAAAAATTCCCTCCGCGGCGGCGGTCTACCA‑TTTTTTTCCCCGCA‑AAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCa                                                                     <  2:227020/143‑1 (MQ=255)
                                                                                         tttACCCAA‑CGCCA‑TGTGAAAAAATCCTTCCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCC‑GAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTaa                                                      <  1:486689/150‑1 (MQ=255)
                                                                                         ttaacccAG‑CGCCA‑TGTGAAAAAAGCCCTGCGGGGGGGTTTTC‑GCCAGGTTATTTTCCCGCC‑ACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTaa                                                      <  2:202815/145‑1 (MQ=255)
                                                                                          ttcatcAA‑CCCCATTGTGAAAAAATCCCTCCGGGGCGGCTCTC‑GCCC‑GGTTTTTTTCCCCCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                  <  2:696274/134‑1 (MQ=255)
                                                                                                aa‑CGCAA‑TGTGAAAAAATGCTTTCGCGGCGGCTTTCTCCCA‑GTTTTTTTCCCGCC‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtgtg                                               <  1:85798/150‑1 (MQ=255)
                                                                                                    cgcggttcacAAAAATCCCTGCGGGGGGTCTGTCGGCGT‑GTTTTTTTCCCGCC‑AAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGa                                             <  2:535933/141‑1 (MQ=255)
                                                                                                      cattgtgAAAAAAGCCTTTCGGGGGGGCTTTCTCCCA‑GTTTTTTTCCCGCC‑GAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg                                          <  1:841424/148‑1 (MQ=255)
                                                                                                              aaaaaTGCCTTCGCGGCGGTTGTCTCCCA‑GTTTTTTTCCCGCC‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAg                                     <  1:285316/148‑1 (MQ=255)
                                                                                                                    gCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAGCCTTTAATCCCG‑CCCGTTTTGTGGGCATAGCTGCTGGGTGAAAATG‑Gggaagggaaggggggg                               >  1:393016/1‑142 (MQ=255)
                                                                                                                            ggcggcTTTCCGCCT‑GTTTTTTCCCCGCC‑AAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                         <  2:418646/146‑1 (MQ=255)
                                                                                                                                       gCCA‑GTTTTTTTCCCGCC‑GCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc              <  2:225829/146‑1 (MQ=255)
                                                                                                                                         caattttttttcCCGCC‑GAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTCGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc            <  2:143612/140‑1 (MQ=255)
                                                                                                                                                ttttcCCGCC‑ACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg  <  1:579459/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357735‑358025

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