Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 773,262 G→C 66.2% T94T (ACC→ACG ymfB ← bifunctional thiamin pyrimidine pyrophosphate hydrolase and thiamin pyrophosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE773,2620GC66.2% 45.9 / 34.3 49T94T (ACC→ACGymfBbifunctional thiamin pyrimidine pyrophosphate hydrolase and thiamin pyrophosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (10/6);  new base C (18/14);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGG  >  minE/773201‑773324
                                                             |                                                              
ccTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCggg                                                          >  1:865386/1‑68 (MQ=255)
      aaaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAttt                                                   >  1:2369576/1‑69 (MQ=255)
      aaaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAttt                                                   >  1:1177546/1‑69 (MQ=255)
      aaaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAttt                                                   >  1:742889/1‑69 (MQ=255)
      aaaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAttt                                                   >  1:742162/1‑69 (MQ=255)
      aaaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAttt                                                   >  1:1777017/1‑69 (MQ=255)
      aaaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAtt                                                    >  1:2207512/1‑68 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTtctca                                              <  1:3293766/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:1701137/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:2494680/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:2750325/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:2754399/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:335960/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:1584202/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:630850/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:2878828/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:733677/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:786174/65‑1 (MQ=255)
               cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgatca                                              <  1:2873422/65‑1 (MQ=255)
                 gcgcTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:2286789/63‑1 (MQ=255)
                            aCGGCAGCAGTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctca                                              <  1:1941613/52‑1 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTgctcaagctcaa                                       >  1:3139162/1‑52 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAAc                        >  1:910562/1‑67 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:321959/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:2959853/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:3318071/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:3331435/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:685514/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:2618549/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1762690/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1116552/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1383486/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:141898/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:2817581/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1938840/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:209013/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:2758542/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGGCTGAGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1407909/1‑68 (MQ=255)
                                   cagTCGATATCGCTGTCATGAGCCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1097253/1‑68 (MQ=255)
                                   cagTCGACATCGCTGTCATGAGGCTGCGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACg                       >  1:1857664/1‑68 (MQ=255)
                                              gCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGt           >  1:2215806/1‑69 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  >  1:2019321/1‑69 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  <  1:318579/69‑1 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  <  1:419197/69‑1 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  <  1:1548059/69‑1 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  <  1:2896164/69‑1 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  <  1:1334319/69‑1 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTgg  <  1:119405/69‑1 (MQ=255)
                                                       aGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTg   >  1:1272177/1‑68 (MQ=255)
                                                             |                                                              
CCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGG  >  minE/773201‑773324

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