Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA41.9% 30.8 / 17.2 31L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (12/6);  new base A (0/13);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.42e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035
                                                                                                                                                        |                                                                                                                                                    
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                      >  2:490167/1‑151 (MQ=255)
 tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                      >  2:527592/1‑150 (MQ=255)
  gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                     <  2:749787/150‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                  <  1:1015633/150‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                 <  1:290396/151‑1 (MQ=255)
            gcgcATCAGGAATTTCAGGCCGCAGACCTGTGTCAGCGTCTCCCGTGTTCCATTTCGGCCAACAC‑CCGCGCATTTATCCACCGCCATGTGCAAAAATGCCTCCGGGGCGG‑CTGTCTGCCGGCTATTTTCCCGCCGCAACAAAAccgccg                                                                                                                                            <  1:69528/149‑1 (MQ=18)
             cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                         >  2:182004/1‑151 (MQ=255)
             cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                         >  2:992792/1‑151 (MQ=255)
               cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                       >  2:1020948/1‑151 (MQ=255)
                aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                                      >  2:913867/1‑151 (MQ=255)
                            cAGCCGCAGACTTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAC‑CCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGGGG‑CTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                            <  2:742618/149‑1 (MQ=21)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg                                                                                                                        >  2:538379/1‑151 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg                                                                                                                        >  1:981118/1‑151 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGCCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAACGCCGTTCTCTCCAGTAGTGGC‑gg                                                                                                                        >  1:565187/1‑151 (MQ=255)
                                      cttgggtCGGCGGCTCCGGTGATTCATTCCCCCCAACAC‑CCGCGCATTTACCAACCGCCATGTGAAAAAATCCTTCCGCGGCGG‑CTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                              <  2:261121/132‑1 (MQ=255)
                                      ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                      >  2:1040925/1‑144 (MQ=255)
                                             cggcgcCTCCCGTGATTCTTTCCCCCAAACAC‑CCGCGCTTTTACCAAACCCCTTGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATa                                                                                                         <  2:618777/146‑1 (MQ=255)
                                                   cgaCCGTGATTCATTCCCCCAAACAC‑CCGCGCTTTTACCAAACGCAGTGTGCAAAAATGCCTTCGCGGCGG‑CGGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTaa                                                                                                  <  2:124193/147‑1 (MQ=255)
                                                         gtgtttaTTTCCCCCCACAACCCGCCCTTTTCCCAACCGCAATTTGAAAAAATCCTTTCGGGGCGG‑CTTTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTc                                                                                           <  1:576754/145‑1 (MQ=255)
                                                               aTTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGGGG‑AAAAACTTTTAAACCCGCCAttt                                                                                       >  2:710721/1‑147 (MQ=25)
                                                                  tCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGG‑AAGGG‑AAAATGTTTTAATCCCgcctatttt                                                                                     >  2:889254/1‑140 (MQ=25)
                                                                   ccGCCACCAC‑CCGCGCTTTTACCAACCGCCTTGTGAAAAAATGCTTCCCGGGCGG‑CTTTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTtgtg                                                                                   <  2:883859/148‑1 (MQ=255)
                                                                   ccGCCAACAC‑CCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTCCGCGGCGG‑CTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                                <  2:1052872/151‑1 (MQ=255)
                                                                       ccacaacccGCCCTTTTCTCCACCCCCATGGGAAAAAATGCCTCCGCGGCGG‑TTGTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                               <  1:24456/142‑1 (MQ=255)
                                                                         accacccGCCCTTTTACCCAACGCCATGTGCAAAACTGCCTCCGCGGCGG‑CTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                            <  1:722332/145‑1 (MQ=255)
                                                                           aa‑cCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGG‑GAGAG‑AAAAAAATTTTAAGACGTCATGTTGTGTgggcgtct                                                                          >  1:302348/1‑143 (MQ=25)
                                                                                          ccACCGCCTTGTGAAAAAATCCCTCCGGGGCGG‑CTTTCTGCCACTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                            <  2:601488/149‑1 (MQ=255)
                                                                                              ccccATGTAAAAAAATGCCTGCGCGGCGG‑CTGTCTCCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                      <  2:1006391/149‑1 (MQ=255)
                                                                                              cGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTCCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                      <  1:494925/151‑1 (MQ=255)
                                                                                                   tgtgCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGG‑AAGAGAAAAAAACTTTAA‑CCCCCCCTGTTGTGTGGGGAACGCCGCTGGGTACaaaaagggaa                                                     >  2:369196/1‑139 (MQ=255)
                                                                                                         aaaaGGCCTTCGGGGGGGCCTGCCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                            <  1:745135/151‑1 (MQ=255)
                                                                                                                    ggggcgg‑cTGTCCCCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                                 <  1:866264/148‑1 (MQ=255)
                                                                                                                         gg‑cTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                             <  2:244550/149‑1 (MQ=255)
                                                                                                                             tttCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                          <  1:578633/147‑1 (MQ=255)
                                                                                                                                                   cacaaaaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCaa  <  2:499825/149‑1 (MQ=255)
                                                                                                                                                      aacaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAAAAAAATTAAACCCGTCATGTTGTGTTGCGATCGCTGCTGGGGAAAATGTGGGAAGGGGAGGGCGGGGGAGCAACTCCAGAAACACCCTAATTGAGTTCGcc      >  1:162129/1‑144 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035

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