Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,315,525 G→A 72.3% Q41Q (CAG→CAA 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,5250GA72.3% 33.3 / 11.9 36Q41Q (CAG→CAAhisAN‑(5'‑phospho‑L‑ribosyl‑formimino)‑5‑amino‑1‑ (5'‑phosphoribosyl)‑4‑imidazolecarboxamide isomerase
Reads supporting (aligned to +/- strand):  ref base G (8/2);  new base A (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 = 9.18e-01

ACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAACGTCAAATCCCGCT  >  minE/1315464‑1315593
                                                             |                                                                    
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                                 >  1:2108553/1‑67 (MQ=255)
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                                 >  1:1911531/1‑67 (MQ=255)
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                                 >  1:1736672/1‑67 (MQ=255)
aCGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGc                                                                 >  1:927830/1‑67 (MQ=255)
      aacaacGCGATTACGGTAACAACCCGCTGCCGCGCTTACAGGATTACGCTGCGCAAGGTGCCGAAGtgt                                                         >  1:1028532/1‑68 (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)
      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:3102762/69‑1 (MQ=255)
                                                   aCGCTGCGCAAGGTGCCGAAGTGTTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAGCGTc            <  1:522825/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)
                                                             gggTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAACGTCAAATCCCGCt  <  1:651347/69‑1 (MQ=255)
                                                             gggTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAACGTCAAATCCCGCt  <  1:2160797/69‑1 (MQ=255)
                                                             |                                                                    
ACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCTAAACGTCAAATCCCGCT  >  minE/1315464‑1315593

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: