Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA52.4% 17.6 / 35.3 21L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (10/0);  new base A (0/11);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 2.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05

AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357745‑358028
                                                                                                                                                  |                                                                                                                                            
aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                         >  1:442555/1‑151 (MQ=255)
                                  cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag                                                                                                         >  1:665246/1‑149 (MQ=255)
                                     tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGCGaaaa                                                                                                    >  1:159741/1‑149 (MQ=255)
                                        cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGaaaaa                                                                                                   >  1:48590/1‑146 (MQ=255)
                                         aGCGGCTACCGTGATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCAATGTGCAAAAATGCCTTCGGGGGGTCTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                       <  1:486819/143‑1 (MQ=21)
                                               tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGgaaaaaacttaaaac                                                                                         >  2:161111/1‑139 (MQ=25)
                                                 ccGTGTTTCATTCCCCCCAACACCCGCGCTTTTCCCCACCCCAAGGTGAAAAAATCCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                         <  1:127859/149‑1 (MQ=255)
                                                 ccGTGATTCATTCCCCCCACCACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                        <  1:495704/150‑1 (MQ=37)
                                                   gtattccttcccccccAACACCCGCGCTTTTACCCAACCCCATTTGCAAAAAGCCCTCGGGGGGGGTTGCTGGCCGTTTTTTTTCCCGCCCCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                       <  2:405949/135‑1 (MQ=255)
                                                   gTGATTCATTCCCGACACCACCCGCCTATTTACCCACCCCCTTGGGCAAAAAGTCTTGCGGGGCGGCTGTCTGCCATTTTTTTTCCCACCAAATAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                              <  2:170574/142‑1 (MQ=255)
                                                     ttttCTTTCCCCCCACCACCCGCCCTTTTCCCCAACGCCTTGGGAAAAAATCCCTCCGCGGCGGCTGTCTCCGAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                                    <  2:227747/148‑1 (MQ=255)
                                                     gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAaaaaa                                                                                                    >  1:65965/1‑133 (MQ=255)
                                                                                    ccacccccATGTGAAAAAATGCCTCGGCGGCGGCTGTCGCCAAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                    <  1:469587/145‑1 (MQ=255)
                                                                                               tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCCCCCTTTTGTGTGGGAATCGCTGCTGGtggaaatgt                                                    >  1:192424/1‑132 (MQ=255)
                                                                                                                         cTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAAACCGCCCTGTTGTGTTGGAAACGCTGCTGGGGAAAATATGGGAAGGGGAGGGGGGGGGATCAGcccc                       >  1:120951/1‑142 (MQ=255)
                                                                                                                           gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAAAAAAAAATTTAAAACCCCCACTTTTTGTGGGAAATGCTGCCGGGGACAAAGTGGGAAAGGGAAGGGGGGGGaacaaaacccacaaac                 >  1:646819/1‑133 (MQ=255)
                                                                                                                             cAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  1:452029/143‑1 (MQ=255)
                                                                                                                               ttttttttcCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa          <  2:766793/145‑1 (MQ=255)
                                                                                                                                  tttttccccgccaCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:297687/143‑1 (MQ=255)
                                                                                                                                          cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGAATAGCTGCCGGGGAAAAAGTGGGAAAGGGAGGGCGGGGGAACAAATCCCAAAACCCCCTAATTGAggt   >  1:362344/1‑144 (MQ=255)
                                                                                                                                            ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  2:777677/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357745‑358028

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