Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,247,208 C→T 65.3% A307A (GCC→GCT thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,247,2080CT65.3% 45.2 / 38.6 69A307A (GCC→GCTthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (9/15);  new base T (23/22);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACTT  >  minE/2247149‑2247276
                                                           |                                                                    
ccGGTGCCGATCGGTACAGTGCCGATCTACAAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCtt                                                             <  1:919551/69‑1 (MQ=255)
     gCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGaa                                                                   >  1:3325659/1‑58 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTTGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:429437/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAATGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3394317/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAATGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3398121/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:636367/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:1054879/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:224455/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:2293447/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:462792/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:344877/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:2634740/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3392697/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:2850322/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3029456/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3059565/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3134693/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3149009/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:1558868/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:1808382/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:854296/68‑1 (MQ=255)
        gATCGGTACAGTGCCGATCTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:1808094/68‑1 (MQ=255)
           cGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGc                                                  >  1:2986796/1‑69 (MQ=255)
           cGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGc                                                  >  1:80103/1‑69 (MQ=255)
                  gTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCcgc                                           >  1:1277055/1‑69 (MQ=255)
                  gTGCCGATCTACAAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAg                                                                  >  1:609717/1‑46 (MQ=255)
                  gTGCCGATCTAACAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCcgc                                           >  1:3306537/1‑69 (MQ=255)
                      cGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCtt                                                             <  1:1271413/47‑1 (MQ=255)
                          cTACCAGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:1457261/50‑1 (MQ=255)
                               aGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTgg                              <  1:1854687/69‑1 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTcc                                             >  1:3276893/1‑54 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTTCGCGACACCCTGCTgg                              >  1:1388572/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACcctgctcg                              >  1:1174307/1‑67 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:3196795/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:3045386/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:3245168/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:3173541/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:3012144/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:1899298/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:2897907/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:137709/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:2872015/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:3385221/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:1465871/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:1781720/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:437203/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:2414550/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:568201/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:1827147/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:2221891/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:1990573/1‑69 (MQ=255)
                               aGGCACTGGAGAAGGTGAACGGGATCGCTGAAGATATTACCTGGGAAGTCTTCCGCGACACCctgctg                               >  1:970993/1‑68 (MQ=255)
                               aGGCACAGGAGAAGGTGAACGGGATCGCTGAAGATCTTACCTGGGAAGTCTTCCGCGACACCCTGCTgg                              >  1:435419/1‑69 (MQ=255)
                                      ggAGAAGGTGAACGGGATCGCTGAAGATCTTACCTggg                                                      <  1:3170813/38‑1 (MQ=38)
                                               tAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAgg              <  1:1745843/69‑1 (MQ=255)
                                               tAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAgg              <  1:2567208/69‑1 (MQ=255)
                                                aaCGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGgt             <  1:1099921/69‑1 (MQ=255)
                                                aaCGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGgt             <  1:1328570/69‑1 (MQ=255)
                                                    ggATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCGGAGCAAGGTGTgg         >  1:1708128/1‑69 (MQ=255)
                                                    ggATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTgg         >  1:2740351/1‑69 (MQ=255)
                                                    ggATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTgg         >  1:2618528/1‑69 (MQ=255)
                                                      aTCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGgt             <  1:321846/63‑1 (MQ=255)
                                                      aTCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGAt       <  1:729372/69‑1 (MQ=255)
                                                        cGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTa     <  1:586217/69‑1 (MQ=255)
                                                        cGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTa     <  1:3335003/69‑1 (MQ=255)
                                                        cGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTa     <  1:3071361/69‑1 (MQ=255)
                                                        cGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTa     <  1:2031727/69‑1 (MQ=255)
                                                        cGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTa     <  1:3330244/69‑1 (MQ=255)
                                                        cGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTa     <  1:3416920/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAAGAGGCCGAGCAAGGTGTGGATTACtt  <  1:264543/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:2420610/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:2353472/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:2307281/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:2893574/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:196868/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:2567085/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:760407/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:29074/69‑1 (MQ=255)
                                                           cGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACtt  <  1:2706350/69‑1 (MQ=255)
                                                           |                                                                    
CCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAGATCTTACCTGGGAAGCGTTCCGCGACACGCTGCTGGAACAGGCCGAGCAAGGTGTGGATTACTT  >  minE/2247149‑2247276

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: