Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,668,958 T→G 51.5% G112G (GGT→GGG malP → maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,668,9580TG51.5% 10.0 / 71.3 66G112G (GGT→GGGmalPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (21/11);  new base G (16/18);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTATGACATCAATCTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTC  >  minE/2668903‑2669012
                                                       |                                                      
cTTATGACATCAATCTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTgg                                           <  1:2984545/69‑1 (MQ=255)
cTTATGACATCAATCTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTgg                                           <  1:2440514/69‑1 (MQ=255)
cTTATGACATCAATCTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTgg                                           <  1:1696829/69‑1 (MQ=255)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACgg                                                   >  1:3399072/1‑52 (MQ=255)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTgg                                           >  1:394976/1‑60 (MQ=255)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:2316240/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:2788214/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:228868/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:1727355/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:3302335/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:2305251/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:2774253/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:2396242/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:146719/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:415492/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:775628/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:1327124/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTgg                                  >  1:908233/1‑69 (MQ=38)
         tCAACCTGACTGACCTGTTGGAAGAAGAGATCGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTg                                   >  1:2032967/1‑68 (MQ=255)
              cTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCt                                             >  1:2093861/1‑53 (MQ=255)
              cTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTggcggc                              >  1:2080996/1‑68 (MQ=255)
              cTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGACTGGGACGTCTgg                                  >  1:268265/1‑64 (MQ=255)
                      ccTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCt                         >  1:3329031/1‑65 (MQ=255)
                      ccTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCt                         >  1:2490389/1‑65 (MQ=255)
                      ccTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCt                         >  1:1864458/1‑65 (MQ=255)
                            ggAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTc               >  1:3272042/1‑69 (MQ=255)
                            ggAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTc               >  1:2157809/1‑69 (MQ=255)
                               agaagaGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAAt            <  1:1578326/69‑1 (MQ=255)
                               agaagaGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAAt            <  1:802237/69‑1 (MQ=255)
                               agaagaGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAAt            <  1:1268625/69‑1 (MQ=255)
                                gaagaGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTc               >  1:2080344/1‑65 (MQ=255)
                                gaagaGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTc               <  1:3354730/65‑1 (MQ=255)
                                gaagaGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTc               >  1:3400941/1‑65 (MQ=255)
                                 aagaGATCGACCCGGCGCTGGGTAACGGTTGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTc               <  1:1607216/64‑1 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:2943862/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:1125423/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:2717881/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:925839/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:2358094/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:1305457/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:1379072/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:1508591/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:1689459/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:1603952/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGt                            >  1:2564093/1‑49 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCa        <  1:1010379/69‑1 (MQ=255)
                                   gagaTCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCa        <  1:2113785/69‑1 (MQ=255)
                                    agaTCGACCCGGCGCTGGGTAACGGTGTTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCa        <  1:348721/68‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:1544793/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:990753/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:696264/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:1379961/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:374988/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:1672163/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:185356/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:2856922/69‑1 (MQ=255)
                                        cGACCCGGCGCTGGGGAACGGTGGGCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGt   <  1:3149847/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCACTGGCAACTGTc  <  1:1164305/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:1724003/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:3159480/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:644968/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:700714/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:1971166/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:215778/69‑1 (MQ=255)
                                         gACCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:994335/69‑1 (MQ=255)
                                          aCCCGGCGCTGGGGAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTc  <  1:3087095/68‑1 (MQ=255)
                                                       |                                                      
CTTATGACATCAATCTGACGGACCTGCTGGAAGAAGAGATCGACCCGGCGCTGGGTAACGGTGGTCTGGGACGTCTGGCGGCGTGCTTCCTCGACTCAATGGCAACTGTC  >  minE/2668903‑2669012

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