Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA50.0% 15.8 / 24.3 34L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (14/3);  new base A (0/17);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 9.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC  >  NC_000913/357742‑358031
                                                                                                                                                     |                                                                                                                                                
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                 <  2:161219/150‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                              <  2:122449/150‑1 (MQ=255)
            gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                    >  1:617630/1‑151 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                     >  1:788951/1‑151 (MQ=255)
                             ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                   >  2:389249/1‑150 (MQ=255)
                              cgcgacccTGTGTCGCCGGCTACCTTGTTTCTTTCCCGCAACCACCCGCGCATTAACCAACCCCCAGGTCAAAAAATCCCTTCGGGGCGGCTGCCGCCCAGTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                                 <  2:407285/145‑1 (MQ=255)
                                  gACCTGTGTCAGCGGCTCCGGTGTTTATTTCCCGCCAACACCCGCCCTTTTCCCCACCCCAATGTGAAAAAATCCCTCCGCGGCGGCTGTCTCCCGGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                  <  1:196613/146‑1 (MQ=255)
                                    cttttgtCAGGGGCTACCGTGATTCTTTTCCCCCAACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATGCCTTCGCGGGGGCTTTCTGCCAGTTATTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                           <  2:758470/146‑1 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGaaa                                                                                                       >  2:407010/1‑150 (MQ=255)
                                         gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                  >  1:845445/1‑139 (MQ=255)
                                           cAGCGGCCACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAACGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGGAACAAAACCGCCGTTCTCTCCAGTACTGGCGGTCGAAgagacgaaga                                                                                                   >  1:392705/1‑144 (MQ=21)
                                                     cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAagagggaaaagct                                                                                                  >  1:364790/1‑133 (MQ=255)
                                                     cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGCGaaaaaaa                                                                                                   >  2:612424/1‑136 (MQ=255)
                                                               ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGGAAAAAACTTTAAACCCgcc                                                                                       >  2:259955/1‑141 (MQ=25)
                                                                         aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGGAAAAAATTTTATCCAGTCAtgttgttttggggac                                                                       >  1:619834/1‑139 (MQ=25)
                                                                           ccGCGCTTTCCTCCACCGCCATGTAAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                    <  1:99916/149‑1 (MQ=255)
                                                                                  tttccccACCGCCATGTGAAAAATTGCTTTCGGGGCGGTTGTCTGCGAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                             <  2:766787/144‑1 (MQ=255)
                                                                                  tttACCCACCCCCATTTAAAAAAAGCCCTCCGGGGCGGCTTTCGCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                              <  2:145730/148‑1 (MQ=255)
                                                                                       ccAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAGCTTTAAACCCTCCTGTTTTGTGGGGATCggctgtgggtaaa                                                          >  2:420469/1‑135 (MQ=25)
                                                                                                     aaaaaTCCTTCCGCGGCGGTTGTCTCCCACCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                           <  2:439233/148‑1 (MQ=255)
                                                                                                     aaaaaTCCCTCCGCGGGGCCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                          <  2:3894/149‑1 (MQ=255)
                                                                                                       aaaaGCCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                       <  2:243066/151‑1 (MQ=255)
                                                                                                       aaaTGCCTCCGCGGCGCCTTTCGCCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                      <  2:294702/151‑1 (MQ=255)
                                                                                                        aaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGACGAGAAAAAGATTTTAACCCCCCATGTTGTGTGGCCATTGCCGCTGGTGAAAAAGTGggaagggaagggggg                                      >  2:406343/1‑145 (MQ=255)
                                                                                                        aaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                     <  2:618075/151‑1 (MQ=255)
                                                                                                                 cgcgGCGGCTGTCTCCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                             <  1:215249/151‑1 (MQ=255)
                                                                                                                  gcggcgGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                             <  2:306914/150‑1 (MQ=255)
                                                                                                                   cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCTCCTTTTGTGTGGGGATCCCTGCTGGTGAAAAAGTGGGAAGGGGAgggggggggagc                                >  2:215273/1‑135 (MQ=255)
                                                                                                                    gggggcttcctcccAGTTTTTTTCCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                             <  1:673827/136‑1 (MQ=255)
                                                                                                                    ggcggcTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                           <  2:458246/150‑1 (MQ=255)
                                                                                                                       ggcTGTCTCCCAGTTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                        <  1:660753/150‑1 (MQ=255)
                                                                                                                               cccgttaaTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                      <  2:362261/137‑1 (MQ=255)
                                                                                                                                ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                     <  2:756897/144‑1 (MQ=255)
                                                                                                                                     tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAAGCCCTCCTGTTGTTTGGGGAACGCTGCTGGGGAAAAAGGGGGAAGGGGAGGGGGGGGGAGCAGCCCCCCAACCCCGCCAAtt           >  2:348568/1‑149 (MQ=255)
                                                                                                                                        tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt       <  1:844925/150‑1 (MQ=255)
                                                                                                                                            cccccacaaaaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc  <  2:701197/145‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC  >  NC_000913/357742‑358031

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: