Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA41.7% 48.7 / 15.5 36L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/4);  new base A (0/15);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 7.21e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357743‑358033
                                                                                                                                                     |                                                                                                                                                  
aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                  <  2:494478/150‑1 (MQ=255)
     cGTAGCGCATCAGACAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                             <  2:258556/150‑1 (MQ=255)
      gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                            <  1:283352/150‑1 (MQ=255)
       tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                          >  1:633252/1‑151 (MQ=255)
                         cgcccGCAGCCCTGTGTGGGCGCTTGCGGGATTTCATTCCCCCAACCACCCGCCCTTTTTCCCACCGCCTTGTAAAAAAATGCTTTCGGGGCGGTTGTCGGCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                       <  1:713337/148‑1 (MQ=255)
                           gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                     >  1:946137/1‑151 (MQ=255)
                             cGCAGACTTGTTCCAGCGCCTACCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCCACCGCCGCGCGAAAAATTGCCCCCGCGGGGGCTGTCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                       <  2:80725/147‑1 (MQ=255)
                               cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                    >  2:344979/1‑148 (MQ=255)
                                      gtgtCAGCGCCAACGGGGTTTAATTCCCCCCAACACCCGCGATTTTACCCACCCCCATGTCAAAAAATCCTTCCGCGGCGGCTTTCGGCCAGCTTTTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                           <  1:398476/149‑1 (MQ=255)
                                       ttccggcgcttcccGTGTTTCATCCCCGCAACCACCCGCGCTTTCCCCCACCGCCATGTGAAAAAAGCCCTGCGGGGCGGTTGTCTGCCAGTTATTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                            <  2:107420/135‑1 (MQ=255)
                                         ccAGCGTCTCCCGTGATTCTTTCCCGCAAACACCCCCGCTTTTACCCACCCCCAGGGAAAAAAATGCCTGCGGGGGGTTTGTCTGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                                      <  2:677918/150‑1 (MQ=255)
                                          cggCGGCTACCGTTTTTCATTCCCCCCAACACCCGCCTTTTTACCCACCCCCATGTCAAAAAACGCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                       <  2:313211/147‑1 (MQ=255)
                                          cAGCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGCATTTACCCAACGCTATGTCAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                          <  1:639819/146‑1 (MQ=255)
                                           ccgccgcctccctttttCCTTTCCCCCCACCACCCCCCCTTTTCCCAACCCCCATGTGAAAAAATCCCTCCGGGGCGGTTGTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                                      <  2:668496/136‑1 (MQ=255)
                                               gcctcccgtatttCTTTCCCCCCACCAACCGCCCATTACCCCAACCCCATGTAAAAAAACCCTTCCGCGGGGGCTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                                   <  2:365774/138‑1 (MQ=255)
                                               gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattt                                                                                                   >  1:14579/1‑141 (MQ=255)
                                                cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGgaaagaaaaaaaatt                                                                                                   >  1:873355/1‑140 (MQ=25)
                                                 tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaatttttga                                                                                               >  1:723870/1‑139 (MQ=25)
                                                 tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattt                                                                                                   >  1:161260/1‑139 (MQ=255)
                                                  cccGTTATTCTTTCCCCCCAACACCCGCGTTTTTTCCCACCGCCATGTACAAAAATGCCTCCGCGGCGGTTTTCGGCCACTTTTTTTCCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                             <  2:797092/150‑1 (MQ=255)
                                                     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGAGAAAAACTTTTAAAcc                                                                                              >  1:893819/1‑148 (MQ=25)
                                                            cATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaacttt                                                                                                  >  2:379016/1‑129 (MQ=25)
                                                              ttccccccAACACCCCCGCTTTTACCCACCCCCATGTAAAAAATTCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                           <  2:895394/136‑1 (MQ=255)
                                                                 ccGCCAACACCCGCGTATTTACCCACCGCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                               <  1:163057/151‑1 (MQ=255)
                                                                         aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAGCCCCTCCTGTTGTGTggggatagc                                                                        >  1:235069/1‑143 (MQ=25)
                                                                           ccGCGCTTTTACCCACCGCCTTGTAAAAAAAGGCTTCCGCGGCGCCTGTCTGCAAGTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                     <  2:763021/151‑1 (MQ=255)
                                                                           ccGCGCTTTTACCCACCCCCTTGTCAAAAAATGCCTCCGCGGCGGTTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                        <  1:941688/148‑1 (MQ=255)
                                                                                    tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCCCCCTGTTGTGTGGGGATCGCTGCTggggaca                                                             >  1:12474/1‑145 (MQ=25)
                                                                                                aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTTCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCCCCTTTTGTGTGGGAAACGCCGCGGGTGAAAAAATGGGAAggg                                                >  1:422698/1‑151 (MQ=255)
                                                                                                  gtgAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTATTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                <  1:26382/149‑1 (MQ=255)
                                                                                                   tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAAATTAAATACGTCCTGTTGTGTggggctcgccgccg                                                                   >  2:135974/1‑118 (MQ=255)
                                                                                                            gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTTAAACCCCCATGTTGTGTGGGGAACGCTGCTGGGGACAAAATGGGAAGGGGAgggggggg                                      >  1:143206/1‑142 (MQ=255)
                                                                                                             ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAATTTTAACCCCCCCTTTTTTTGGGGAATCGCTGGTGGGGACAAAATGGGAAGGGAAGGggggcg                                      >  2:546068/1‑148 (MQ=255)
                                                                                                                 cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAACTTTAAAACCGTCCTGTTTTGTGGGGATTGCTGCGGGGGACAAAATGGGAAGGGaggggggggg                                      >  2:519447/1‑135 (MQ=255)
                                                                                                                                   acttttttcCCGCCGCAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                <  2:150726/144‑1 (MQ=255)
                                                                                                                                          tcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt        <  1:165736/147‑1 (MQ=255)
                                                                                                                                              gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGAATCCCCGCTGGGGACAAATTTGGAAGGGGAAGGGGGGGGATCAAATCCAGAACCACGCCAAATTAGTTTGcccc  >  1:269725/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357743‑358033

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: