Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA42.3% 34.2 / 11.1 26L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (12/3);  new base A (0/11);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357742‑358035
                                                                                                                                                      |                                                                                                                                                     
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                      <  2:1029326/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCTGCCGCAGTCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                     <  1:967669/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                     >  1:908025/1‑150 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                 >  2:431597/1‑151 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                 <  2:1029039/151‑1 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                              >  1:956022/1‑151 (MQ=255)
       ctAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                  <  1:793542/146‑1 (MQ=255)
        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                             >  1:785286/1‑151 (MQ=255)
                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                     >  2:790741/1‑150 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                           >  1:134257/1‑151 (MQ=255)
                                  gtccTGTGTCCGCGGCTCCCGTGATTCTTTCCCCCAAACCACCGCTATTTTACCCAACGCGATGGGCAAAAAAGCGTTGGGGGCGGTTGTCTGCAACTTATTTTCCGGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                                   <  2:117488/148‑1 (MQ=255)
                                    ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                        >  2:733225/1‑144 (MQ=255)
                                      tgtgTCAGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                               <  2:261310/149‑1 (MQ=255)
                                        tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGGAGTGGCggggggag                                                                                                                    >  1:68381/1‑140 (MQ=255)
                                        tgacggCGGCTCCGGTGTTTCTTTCCCGCCAACACCCGCATTTTTCCCCACCCCCATGTGAAAAATTGCTTCCGGGGCGGCTGTCTGCCGGTTTTTTTCCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                          <  1:9844/133‑1 (MQ=255)
                                               ggCTCCCGTATTTCATTCCCGCAAACACCCGCGCATTTACCCACCGCCTTGTGAAAAAATGCCTCCGCGGCGGTTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                                       <  2:999392/148‑1 (MQ=255)
                                                             aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAAGAGAAAAAATTTAAAGccgtgactcttt                                                                                       >  1:788053/1‑141 (MQ=25)
                                                                                  tttACCCAACGCAATGAAAAAAAACCCCTCCGGGGCGTCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                  <  1:255967/148‑1 (MQ=255)
                                                                                   ttttCCAACCGCTATGTAAAAAAACGCTTTCGGGAGGGCTTTTTGCAAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                    <  2:839742/143‑1 (MQ=255)
                                                                                       ccACCGCCATGGAAAAAATTGCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                              <  2:372900/147‑1 (MQ=255)
                                                                                                  gtaaaaaaaTCCCTCGGCGGCGCCTTTTTGCAATCTTTTTGCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                                    <  1:751062/143‑1 (MQ=255)
                                                                                                    taaaaaaTGGCTTCCGGGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGCGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                                <  2:664413/148‑1 (MQ=255)
                                                                                                            gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAGCCCCCCATCTTTTGTGGGGAACCCTGCTGGGGAACAATGGGGGAAggaagggggggg                                         >  1:255091/1‑139 (MQ=255)
                                                                                                             ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGATTTAAAACCCTCATTCTTTGTGGCTATCGCTCGCTGTAAAGAATTGGGGAAGGGagggcgggcg                                        >  1:719707/1‑140 (MQ=255)
                                                                                                                     gggggcttccTGCCATTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                                <  2:24668/140‑1 (MQ=255)
                                                                                                                        ggcTGTCTGCCATCTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                             <  2:956036/149‑1 (MQ=255)
                                                                                                                                   aGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAAATTTAAACCCTCCTGTTTTTTGGGAATAGCTGCTGGGGACAAAATGGAGAGGGGAGGGCGGCGAAAAAGCCCCCGAACCCCGCCaa                 >  2:26525/1‑150 (MQ=255)
                                                                                                                                    gtggtttgcccgcagAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGGTTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                <  1:27151/141‑1 (MQ=255)
                                                                                                                                          tttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGGGAAAAAACCTTATGCCCCTCATGTTGTGTTGGAAACGCTGCTGGTAACAATGTGGGAAGGGGAGGGGGGCGGAACAGCTCCCCCACCGCGCCaaatttg            >  2:86141/1‑143 (MQ=255)
                                                                                                                                                  agaacaaAACCGCCGTCCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAGCGCTGCTGGTAACAATGTGGAAAGTGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCGTCaa  <  1:192266/148‑1 (MQ=255)
                                                                                                                                                    aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa   <  2:27866/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357742‑358035

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