Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 262,906 G→A 50.9% S328N (AGC→AAC)  ribD → fused diaminohydroxyphosphoribosylaminopyrimidin e deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE262,9060GA50.9% ‑5.6 / 53.1 60S328N (AGC→AAC) ribDfused diaminohydroxyphosphoribosylaminopyrimidin e deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base G (8/21);  new base A (13/17);  total (21/39)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

GGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAAT  >  minE/262853‑262971
                                                     |                                                                 
ggCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     <  1:805570/68‑1 (MQ=255)
ggCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     <  1:681237/68‑1 (MQ=255)
ggCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     <  1:606549/68‑1 (MQ=255)
  cGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTa                                                   >  1:909796/1‑68 (MQ=255)
  cGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTa                                                   >  1:79081/1‑68 (MQ=255)
  cGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTa                                                   >  1:2526757/1‑68 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:105933/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:882929/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:815978/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:699120/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:681949/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:459184/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:2682532/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:2248675/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:1140111/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:1237110/69‑1 (MQ=255)
    ggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:1549726/68‑1 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     >  1:1178411/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     >  1:2053816/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     >  1:651628/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     >  1:3405524/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                     >  1:3344504/1‑60 (MQ=255)
                                  gCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                 <  1:1789541/38‑1 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGTATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:2246704/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCa                                               >  1:2352902/1‑37 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:1555087/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:2968688/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:2953138/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:805775/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:947957/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:98519/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:180823/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:2406096/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct               >  1:167040/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGc                >  1:1275396/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGACAGGGCTTGAGAAATTAGCCGACGct               >  1:3176691/1‑68 (MQ=255)
                                       tAAACTATTAGGCACCGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:394024/52‑1 (MQ=39)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:639252/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:1023297/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:1961004/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:214493/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:130319/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:1349645/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:1419664/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:661170/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:2159928/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:1466425/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:59789/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:2600578/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:2210993/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:3120997/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:2838953/52‑1 (MQ=255)
                                       tAAACTATTAGGAAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:1924370/52‑1 (MQ=255)
                                           cTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCa          <  1:1590534/68‑1 (MQ=255)
                                           cTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCa          <  1:2469397/68‑1 (MQ=255)
                                           cTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCa          <  1:657702/68‑1 (MQ=255)
                                              ttAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                              <  1:2267179/45‑1 (MQ=255)
                                                  gCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa   <  1:2933968/68‑1 (MQ=255)
                                                  gCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa   <  1:249865/68‑1 (MQ=255)
                                                  gCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa   <  1:854961/68‑1 (MQ=255)
                                                   ccgCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa   <  1:1890739/65‑1 (MQ=255)
                                                     gCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAgccgccgcc               >  1:263577/1‑53 (MQ=255)
                                                     gCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAAt  >  1:1995595/1‑66 (MQ=255)
                                                     |                                                                 
GGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAAT  >  minE/262853‑262971

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: