Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT44.7% 36.8 / 16.1 38N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (16/5);  new base T (0/17);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT  >  NC_000913/357728‑358022
                                                                                                                                                    |                                                                                                                                                     
gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTAt                                                                                                                                                      <  1:601997/149‑1 (MQ=255)
      gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTccgc                                                                                                                                              >  1:1083580/1‑151 (MQ=255)
           ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGca                                                                                                                                          >  2:460062/1‑150 (MQ=255)
                   gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                 >  2:243355/1‑151 (MQ=255)
                    cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                <  1:1209719/151‑1 (MQ=255)
                     gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                               >  1:1014721/1‑151 (MQ=255)
                       aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                             >  2:892951/1‑151 (MQ=255)
                              cAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                         >  1:726604/1‑148 (MQ=255)
                                           ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                          >  1:841571/1‑150 (MQ=255)
                                             gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                         >  1:975387/1‑149 (MQ=255)
                                              cagccttgtgTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                       <  1:51395/144‑1 (MQ=21)
                                              cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                          >  1:84397/1‑147 (MQ=255)
                                                gcccTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCACCACCCGCGCATTTACCCACCGCCATGCCAAAAAATGCCTGCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                     <  2:410939/148‑1 (MQ=255)
                                                    tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgcggg                                                                                                   >  2:278737/1‑143 (MQ=255)
                                                     gtgtCAGCGGCACCCGTAATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCAATGTGAAAAATTGCCTTCGCGGCGGCTGTCTGCCA‑GTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aa                                                                                                <  1:571519/149‑1 (MQ=255)
                                                        tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaggaaa                                                                                                >  2:832191/1‑138 (MQ=255)
                                                             ggCTCCCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCAACCGCCTTGTGCAAAAATCCCTTGGGGGCGGCTGTCTCCCA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                   <  1:542443/139‑1 (MQ=255)
                                                             gcctccgttatttCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCATTGTGAAAAAATGCCTCCGGGGGGTTTGTCTCCCA‑GTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTaa                                                                                     <  2:231432/140‑1 (MQ=255)
                                                                  ccGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑TTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGccc                                                                                 <  1:118006/150‑1 (MQ=25)
                                                                             ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑AAAAAAA‑TTTaaaaccccccctttttg                                                                      >  1:29861/1‑135 (MQ=25)
                                                                               ccGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTCGGGGCGGCTGTCTGCCA‑GCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg                                                                    <  1:706877/150‑1 (MQ=255)
                                                                                 gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAGCTTTTAAAACCG‑CCTGTTGTGTGGcg                                                                 >  1:1096841/1‑150 (MQ=25)
                                                                                         ccGCCCTTTTACCCAACGCATTTTGAAAAAATGCTTCCGGGGCGGCTTTCTGCCA‑GTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                            <  1:1146738/148‑1 (MQ=255)
                                                                                           gcgcATTTATCCACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTTTTTCCCA‑GCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                        <  1:128241/150‑1 (MQ=255)
                                                                                               aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTAAACCCCCCTGTTGTGTGGGGATCGCTGCTgg                                                      >  1:879485/1‑148 (MQ=25)
                                                                                                ttactccacccctaTGTGAAAAAAGGCCTTCGCGGCGGTTGTCGGCCA‑GTTTTTTTCCGGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                       <  1:258390/133‑1 (MQ=255)
                                                                                                     ccACCGCCTTGTAAAAAAATGCCCCCGCGGCGTGTGTCGCCCA‑CCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                             <  2:408038/151‑1 (MQ=255)
                                                                                                       acctccttgtgAAAAAAGGCCTCCGCGGGGGCTGTCTCCCG‑GCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                            <  1:122420/143‑1 (MQ=255)
                                                                                                               gtgCAAAAATTCTTCCGCGGGGGTTGTCTGCCA‑GCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                     <  2:192610/149‑1 (MQ=255)
                                                                                                                tgaaaaatTCCCTTCGCGGCGGCTGTCGCCCA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                   <  1:1200253/147‑1 (MQ=255)
                                                                                                                   aaaaaTGCCTTCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                               <  1:775079/151‑1 (MQ=255)
                                                                                                                   aaaaaTCCCTTCGGGGGGGTTGTCTGCCA‑GTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                               <  1:125454/151‑1 (MQ=255)
                                                                                                                   aaaaaTCCCTTCGCGGCGGTTGTCGCCCA‑TTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                 <  2:1166190/149‑1 (MQ=255)
                                                                                                                    aaaaAGCCCTCCCGGCGGCCTTCCGCCACGCTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                   <  2:660123/147‑1 (MQ=255)
                                                                                                                          ccTCCGGGGCGGCTGTCTGCCA‑GTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                        <  1:857688/151‑1 (MQ=255)
                                                                                                                              cgcgGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAAT‑TTTAATCCCG‑CATTTTTTGGGGGCATAGCTGCTGGGTAAAATGGGGGAAGGGGAgggggggggaacaa                     >  1:943409/1‑137 (MQ=255)
                                                                                                                                             ccA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg          <  1:601815/146‑1 (MQ=255)
                                                                                                                                                  cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAA‑GAAAAGA‑GTTTAAACCGCCCCGTTGTGGGGCAATAGCCGCGTGGGACAAAATGGGAAGGGGAGAGGGGCGGATCAGCTCCACAAACgcgcccatt  >  1:551112/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT  >  NC_000913/357728‑358022

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