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

AGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357746‑358033
                                                                                                                                                   |                                                                                                                                                 
aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                             <  1:396226/151‑1 (MQ=255)
 gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                            <  2:792969/151‑1 (MQ=255)
  cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                           >  1:809644/1‑151 (MQ=255)
            cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                 >  1:165323/1‑151 (MQ=255)
                   ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                          >  2:342744/1‑151 (MQ=255)
                        gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                    >  2:695955/1‑151 (MQ=255)
                                     gtCAGGGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                               <  1:311373/143‑1 (MQ=255)
                                        aGCGGCTAACGTGATTCAAATCCGCCAACAATCGCGCATTGATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATTTTTCGGCCGCAACAAAAACGCCGTTCTGTCCAGTTGTGGCGGGGCAggacacaaagtg                                                                                                     >  2:722223/1‑141 (MQ=18)
                                           ggCTCCCGTTATTCATTCCCCCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCCTCCGCGGCGTCTGTCGCCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                                       <  2:775215/144‑1 (MQ=255)
                                                  gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaaatttaacccc                                                                                           >  1:68215/1‑135 (MQ=25)
                                                   tGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGGAAAAAACTTTatctc                                                                                             >  2:383828/1‑143 (MQ=25)
                                                      ttCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAATTCCCTCGGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                       <  1:784696/149‑1 (MQ=255)
                                                            cccccaaacaccCCCGCTTTTACCCACCCCAATGTCAAAAAAGCCCTCCGCGGCGTCTTTCTCCAAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                               <  2:217394/145‑1 (MQ=255)
                                                              ccgcaacagCCGCCTATTTTACCAACGCCATGTTAAAAAATGCCTTCGGGGGGGGTGTCTTGCCGGCTTTTTTTCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                               <  2:876475/147‑1 (MQ=255)
                                                                     aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACCCCCCCTGTTGTGTgggg                                                                            >  1:796708/1‑143 (MQ=25)
                                                                        cacgcATTTACCCACCGCCGTGAACAAAATTCCTCCCGCGGCGCCTGTCTCCGAGCTGTTTCCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                   <  2:178621/149‑1 (MQ=255)
                                                                                   ccAACGCCATTTCAAAAATTGCTTCCGCGGCGGGTGCCTCCCAGTTATTTTCCCGCCCCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                            <  1:275398/147‑1 (MQ=255)
                                                                                       ccccttgtgAAAAAAGCCCTCCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                    <  1:137734/146‑1 (MQ=255)
                                                                                            tgtcaaaaaaTCCCTGCGCGGCGCCTTTCTCCCGGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                 <  2:332029/144‑1 (MQ=255)
                                                                                             gcaaaaaaaGGCCTTCGCGCCGGCTGCTGCCAGGCTTTTTTCCCGCGGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                     <  2:271228/140‑1 (MQ=255)
                                                                                                 aaaaaTGCTTGCGCGGCGGCTGTCGGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                           <  1:661949/150‑1 (MQ=255)
                                                                                                 aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTTAAACCCCCACGTTTTGTGGGGATCGCTGCTGGGGAAAAAGTGGGAAGGGGAgggg                                          >  1:895985/1‑148 (MQ=255)
                                                                                                  aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAAATTTAACCCCTCCTGTTGTGTGGGGAACGCCGCTGGGGAAAATATGGGAAGGGGAggggg                                         >  1:832003/1‑147 (MQ=255)
                                                                                                               ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTAAACCCCCCATGTTGTGTGGGGAGAGCTGGTGGTGAAAATGTGGGAAGGGGAGAgcgggggaacaacc                              >  1:770323/1‑139 (MQ=255)
                                                                                                                 gcgcTGTCCGCCCTCTTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                 <  2:816091/143‑1 (MQ=255)
                                                                                                                            ccAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                 <  1:294247/150‑1 (MQ=255)
                                                                                                                                             ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc  <  1:45382/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357746‑358033

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