Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,315,519 C→T 72.3% A39A (GCC→GCT hisA → N‑(5'‑phospho‑L‑ribosyl‑formimino)‑5‑amino‑1‑ (5'‑phosphoribosyl)‑4‑imidazolecarboxamide isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,315,5190CT72.3% 34.8 / 11.6 36A39A (GCC→GCThisAN‑(5'‑phospho‑L‑ribosyl‑formimino)‑5‑amino‑1‑ (5'‑phosphoribosyl)‑4‑imidazolecarboxamide isomerase
Reads supporting (aligned to +/- strand):  ref base C (8/2);  new base T (18/8);  total (26/10)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAACGTC  >  minE/1315464‑1315583
                                                       |                                                                
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                       >  1:2108553/1‑67 (MQ=255)
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                       >  1:1911531/1‑67 (MQ=255)
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                       >  1:927830/1‑67 (MQ=255)
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                       >  1:1736672/1‑67 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:3404808/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:649802/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:882186/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:2377863/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:2272251/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:2051924/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1886435/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1854385/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:946040/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1674289/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1542165/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:144146/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1327858/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1269837/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1028532/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAAtgt                                               >  1:2525995/1‑68 (MQ=255)
      aacaacGCGATTACGGTAACAAACCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:1683965/1‑68 (MQ=255)
      aacaacCCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                               >  1:356262/1‑68 (MQ=255)
        cAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGtgctg                                             >  1:2060813/1‑69 (MQ=255)
        cAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGtgctg                                             >  1:1274825/1‑69 (MQ=255)
        cAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGtgctg                                             >  1:3318572/1‑69 (MQ=255)
                       aaCGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTg                              <  1:924048/69‑1 (MQ=255)
                       aaCGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTg                              <  1:12347/69‑1 (MQ=255)
                                          tGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGAt                                 >  1:1747396/1‑47 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:3126171/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:522825/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:3102762/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:3046267/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:2606742/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:1892801/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:972617/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTAACCGGGGCAAAAGATCCGGCTAAGCGTc  <  1:444218/69‑1 (MQ=255)
                                                       |                                                                
ACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAACGTC  >  minE/1315464‑1315583

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: