Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA54.5% 3.6 / 25.7 22L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (8/2);  new base A (0/12);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 1.41e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-02

GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357749‑358028
                                                                                                                                                 |                                                                                                                                           
gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                    >  2:749669/1‑150 (MQ=255)
     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                                                >  2:106096/1‑149 (MQ=255)
              atttcccccggCAGCCCTGTGCCGGCGGCTCCCGTTTTTCTTTCCCCCAACCACCCGCGATTTTACCCACCCCCTTGTACAAAATTGCTTCCGGGGCGGTTGTCTGCCAGTTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                     <  2:153865/141‑1 (MQ=255)
              attcccccccccAGACCTGTGTCGGCGGCTACGGTGATTCTTTCCCCCAACCACCCGCGCTTTTACCCACCGCAGTGTGAAAAAATGCTTTCGGGGCGGCTTTCTGCCAGTTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                     <  2:746321/140‑1 (MQ=255)
              attcccAGCCGCAGCCTTGTGTCGGGGCCCCCCGTGTTTCTTTCCCGCCAACACCCGCACTTTTCCCCAACCCCATGTGAAAAAGTCCCTCCGCGGCGGCTGTCTCCCAGCTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                     <  2:613064/147‑1 (MQ=255)
               attccccgccccgACCTGTGGCGGCGGCTTCCGGTATTTATTTCCCCCAAAAACCCCCGCTTTTCCCCACCGCAAGTTGAAAAATGCCTTCCGGGGCGCCTTTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                     <  1:591658/139‑1 (MQ=255)
                  cccgCCGCAGCCTTGTTCCAGCGCACACCGGGTTTCATTCCCCCCAACACCCGCCCTTTACCCCAACGCCATGTGAAAAATGGCTTTCGGGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGc                                                                                                                 <  1:596116/148‑1 (MQ=255)
                           gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagg                                                                                                        >  2:229730/1‑146 (MQ=255)
                            cccTGTGCCAGGGGCTACCGGGTTTCTTTCCCCCAAACAACCGCCCATTTACCCACCCCCGTGGGAAAAAAGGCCTTCGGGGGGGCTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                       <  2:244679/150‑1 (MQ=255)
                             ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGccaaaaa                                                                                                      >  2:430780/1‑145 (MQ=255)
                                       cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaacgagcaaaaacttt                                                                                            >  1:675878/1‑137 (MQ=255)
                                           ttccGTGATTTCTTTCCCCCCAACACCCGCGCTTTTACCAAACGCTATGTGCAAAATTGCCTCCGGGCCGTTTTCCGCCGGGTTTTTTTCCCCCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                         <  1:141499/149‑1 (MQ=255)
                                             ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAaaaaaaaattttaaacc                                                                                       >  2:217178/1‑137 (MQ=25)
                                                 gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAGATTTTAGccccccc                                                                                   >  2:35545/1‑146 (MQ=25)
                                                      taTTCCCGCCAACACCGGCGATTTAACCCACCCCCATGTGCAAAATTGCCTTCGGGGGGGCGGTCTCCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                              <  1:747483/150‑1 (MQ=255)
                                                                                    ccacccccttgtcAAAAATTCCTTCCGGGGCGGTTGTCTGCCACTTTTTTTCCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                   <  2:744250/139‑1 (MQ=255)
                                                                                                                     ttGTCTCCCATCTTTTTTCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  2:388495/150‑1 (MQ=255)
                                                                                                                     ctttctcccattttttttcCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  2:62492/137‑1 (MQ=255)
                                                                                                                     cTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTCCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  2:363870/151‑1 (MQ=255)
                                                                                                                                ttgtttccccgcccCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:71594/143‑1 (MQ=255)
                                                                                                                                cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTTTTTTTTGGGGCTCGCTGCTGGGGACAAAAAGGGGAGGGAAGGGGGGGGGATCACATCCCGAACCccgccacttt       >  2:116652/1‑145 (MQ=255)
                                                                                                                                     ttcccgccgAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  2:651123/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357749‑358028

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