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

AGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357746‑358033
                                                                                                                                                   |                                                                                                                                                 
aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAAcc                                                                                                                                             <  1:396226/151‑1 (MQ=255)
 gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAAccg                                                                                                                                            <  2:792969/151‑1 (MQ=255)
  cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                                           >  1:809644/1‑151 (MQ=255)
            cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                 >  1:165323/1‑151 (MQ=255)
                   ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                          >  2:342744/1‑151 (MQ=255)
                        gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                                    >  2:695955/1‑151 (MQ=255)
                                     gtCAGGGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                                               <  1:311373/143‑1 (MQ=255)
                                        aGCGGCTAACGTGATTCAAATCCGCCAACAATCGCGCATTGATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTC‑TGCCA‑GCTATTTTT‑CGGCCGCAACAAAAACGCCGTTCTGTCCAGTTGTGGC‑GGGGCAggacacaaagtg                                                                                                     >  2:722223/1‑141 (MQ=18)
                                           ggCTCCCGTTATTCATTCCCCCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCCTCCGCGGCGTCTGTC‑GCCCG‑GTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATa                                                                                                       <  2:775215/144‑1 (MQ=255)
                                                  gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑aaaaaaaatttaacccc                                                                                           >  1:68215/1‑135 (MQ=25)
                                                   tGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGG‑AAAAAACTTTatctc                                                                                             >  2:383828/1‑143 (MQ=25)
                                                      ttCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAATTCCCTCGGGGGCGGCTGTC‑TCCCA‑GTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAt                                                                                       <  1:784696/149‑1 (MQ=255)
                                                            cccccaaacaccCCCGCTTTTACCCACCCCAATGTCAAAAAAGCCCTCCGCGGCGTCTTTC‑TCCAA‑GTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                               <  2:217394/145‑1 (MQ=255)
                                                              ccgcaacagCCGCCTATTTTACCAACGCCATGTTAAAAAATGCCTTCGGGGGGGGTGTCTTGCCG‑GCTTTTTTT‑CCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                               <  2:876475/147‑1 (MQ=255)
                                                                     aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAACTTTTAACCCCCCCTGTTGTGTgggg                                                                            >  1:796708/1‑143 (MQ=25)
                                                                        cacgcATTTACCCACCGCCGTGAACAAAATTCCTCCCGCGGCGCCTGTC‑TCCGA‑GCTGTTTCC‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                   <  2:178621/149‑1 (MQ=255)
                                                                                   ccAACGCCATTTCAAAAATTGCTTCCGCGGCGGGTGCC‑TCCCA‑GTTATTTTC‑CCGCCCCAAAAAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                            <  1:275398/147‑1 (MQ=255)
                                                                                       ccccttgtgAAAAAAGCCCTCCGGGGGGGCTGTC‑TGCCA‑GCTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                    <  1:137734/146‑1 (MQ=255)
                                                                                            tgtcaaaaaaTCCCTGCGCGGCGCCTTTC‑TCCCG‑GTTTTTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                 <  2:332029/144‑1 (MQ=255)
                                                                                             gcaaaaaaaGGCCTTCGCGCCGGCTG‑C‑TGCCAGGCTTTTTTC‑CCGCGGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                     <  2:271228/140‑1 (MQ=255)
                                                                                                 aaaaaTGCTTGCGCGGCGGCTGTC‑GGCCA‑GTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                           <  1:661949/150‑1 (MQ=255)
                                                                                                 aaaaaTGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAAATTTAAACCC‑CCACGTTTTGTGGGGATCGCTGCTGGGGAAAAAGTGGGAAGGGGAgggg                                          >  1:895985/1‑148 (MQ=255)
                                                                                                  aaaaTGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAG‑AAAAAAAATTTAACCCCTCCTGTTGTGTGGGGAACGCCGCTGGGGAAAATATGGGAAGGGGAggggg                                         >  1:832003/1‑147 (MQ=255)
                                                                                                               ggcggcTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAG‑AAAAAAATTTAAACCCCCCATGTTGTGTGGGGAGAGCTGGTGGTGAAAATGTGGGAAGGGGAGAgcgggggaacaacc                              >  1:770323/1‑139 (MQ=255)
                                                                                                                 gcgcTGTC‑CGCCC‑TCTTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                 <  2:816091/143‑1 (MQ=255)
                                                                                                                            ccA‑GTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                 <  1:294247/150‑1 (MQ=255)
                                                                                                                                             ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc  <  1:45382/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTC‑TGCCA‑GCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357746‑358033

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: ‑