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

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357741‑358038
                                                                                                                                                       |                                                                                                                                                      
tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                        >  2:219232/1‑150 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                      <  2:579069/151‑1 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                        <  2:341052/149‑1 (MQ=255)
   tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                     >  2:29825/1‑150 (MQ=255)
      gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                 <  1:370891/151‑1 (MQ=255)
            cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                           >  2:570210/1‑151 (MQ=255)
             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                          >  2:199486/1‑151 (MQ=255)
                             gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                           >  1:169651/1‑150 (MQ=255)
                              ccGCAGCCTTGTGTCAGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                            <  2:254116/148‑1 (MQ=255)
                                   gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                        >  2:325109/1‑146 (MQ=255)
                                     ccTGTGTCGGCGCCTCCCGTGTTTCATTCCCCCCAACACCCGCCTTTTTACCCACCCCCATGTGAAAAAATGCCTGCGGGGCGGCTGTCTGCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                  <  2:278963/150‑1 (MQ=255)
                                     ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAag                                                                                                                     >  2:215879/1‑148 (MQ=255)
                                         gtggcgccgcctccGGGATTCCTTCCCCCCAACAACCGCCCAATTTACCCAACCCCATGCAAAAAAACGCCTGCGGGCCGGCTGTCCCCAACCTATTTGCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                              <  1:29821/139‑1 (MQ=255)
                                            cggCGGCTACCGTTATTCACTCCCCCCAACACCCGCTTATTTACCCACCCCCTTGAGAAAAAACCCCTCCGCGGCGGTTTTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                       <  1:466669/137‑1 (MQ=255)
                                                   tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGCGaaaaaaacttaaaa                                                                                                    >  1:6814/1‑139 (MQ=25)
                                                                  ccGCCAACACCCGCGTTTTTACCCACCCCTTTGTAAAAAATTCCCTCGGCGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                     <  2:454333/150‑1 (MQ=255)
                                                                             cgcgCTTTTACCCACCCCCGTGTCAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGTTATTTCCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                          <  1:53622/150‑1 (MQ=255)
                                                                               cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCCCCATGTTGTGTGGGGATCgc                                                                             >  2:355375/1‑145 (MQ=25)
                                                                                             ccccATGTAAAAAAAGGCTTTCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                         <  2:214872/149‑1 (MQ=255)
                                                                                                   gaacAAAATTGCCTCCGCGGGGTCTGCCCCCAGTTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                           <  1:509728/140‑1 (MQ=255)
                                                                                                        aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATAGCCGCTGGGGAAAAAGTGGGAAGGGGAgaggg                                              >  2:37579/1‑147 (MQ=255)
                                                                                                              ctttcgggcgggcTGTCTCCCTCTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                            <  2:496042/137‑1 (MQ=255)
                                                                                                                   gcggcgGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                                    <  2:251938/150‑1 (MQ=255)
                                                                                                                              cggCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                         <  2:2382/148‑1 (MQ=255)
                                                                                                                                 ccattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                         <  2:92946/140‑1 (MQ=255)
                                                                                                                                      ttttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                      <  1:570091/143‑1 (MQ=255)
                                                                                                                                      tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGCGAGAAATATCTTTAAAACCCCCCTGTTGTGTGGGTATCTCTGCTTGGTACCAAGTGGGAAGGGGAGGGAGGGGAGTCAGCGCCCGaaacccccgcatt                  >  1:286973/1‑138 (MQ=255)
                                                                                                                                            tCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGGAAGCGGAAAAAAATTAAACCCCCCGATCTTTGTGGGAAAAGGCGTTGGTGAGAAAGCGGGAAGGGGAGAGCGGGGAATTAACTCCAGAACcgcgggcgct                  >  1:70982/1‑137 (MQ=255)
                                                                                                                                                     aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:180415/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357741‑358038

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