Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA61.9% ‑6.1 / 19.6 21L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (0/13);  total (4/17)
Fisher's exact test for biased strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-02

GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357749‑358028
                                                                                                                                                     |                                                                                                                                          
gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                   <  2:29102/151‑1 (MQ=255)
            gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                       >  2:547480/1‑151 (MQ=255)
                tcccACCCGCAGCCCTGTGTCGGCGCCTCCCGTGTTTCTTTCCCCCAACAACCCGCGCTTTTACCCAACGCCATGTCAAAAAATCCTTTCGCGGCGGCTTTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  1:486966/146‑1 (MQ=255)
                  ccAGCCGCAGACTTGTGTAGGCGGCTCCCGTGATTCTTTCCCCCAAACACCCGCGATTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   <  1:601154/149‑1 (MQ=255)
                      ccGCAGACTTGTTCCAGCGCCTACCGTGTTTCTTTCCCGCAACCACCCCCGAATTTACCCACCCCCTTGGGAAAAAATGCCTTCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                               <  1:472741/149‑1 (MQ=255)
                       cGCAGCCCTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTGAAAAAATCCCTCCGGGGCGGCTGTCGGCCGCCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                            <  1:453450/151‑1 (MQ=255)
                                gggtCGGCGGCACCCGTATTTATTTCCCCCCAACACCCGCGCTTTAACCCAACGCAATGTCAAAAATGGCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                    <  1:547525/148‑1 (MQ=255)
                                   tCAGCGCCTCCGGTGTTTCTTTCCCCCAACCACCCGCCCATTTCCCCCACCCTTTTTGAAAAACTCCCTGCGGGGGGGTTCTCCGCCAGTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                          <  2:315094/141‑1 (MQ=255)
                                      gcgcCTCCCGTAATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCGCCTTGTAAAAAAATGCCTCCGCGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                     <  2:816482/143‑1 (MQ=255)
                                         gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGgaaaaaatttta                                                                                           >  1:803257/1‑141 (MQ=25)
                                          cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGagggcgagaa                                                                                                 >  2:222198/1‑135 (MQ=255)
                                           caCCGTTATTTGTTTCCCCCCAACACCCGGGCATTTCTCAACGCCAAGGACAAAAATGCCTTGGGGGGGGGCGCTAGCCCGGTTGTTGACCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                         <  2:546411/150‑1 (MQ=255)
                                                            cccccAACACCCGCGCTTTCTCCCACGGCAATGTGAAAAATTCCCTCCGGGGCGTCTTTCCCCTGATTTTTTTCACGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                          <  2:530433/147‑1 (MQ=255)
                                                                                     caaccgctgtgtgAAAAAACGCTTTCGCGGCGTCTGTCTGCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                    <  1:376743/141‑1 (MQ=255)
                                                                                            aaTGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                            <  1:152752/150‑1 (MQ=255)
                                                                                                    aaaaaggcttgggcggcgGTTGTCGCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                      <  1:50201/138‑1 (MQ=255)
                                                                                                    aaaaTTCCTTCCGGGGCGGTTTTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                       <  1:302722/149‑1 (MQ=255)
                                                                                                        tGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAACATTTAACCCCCCATTTTTTTTGGGGAACGCCGCTGGGGAAAAAGTGGGAAGGGGAgggggggg                                 >  2:54284/1‑143 (MQ=255)
                                                                                                                          tCTGCAAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                 <  2:803190/149‑1 (MQ=255)
                                                                                                                                   ttttttcCCCCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:583093/148‑1 (MQ=255)
                                                                                                                                           ccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  1:270224/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 < 17 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: