Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 262,916 C→T 48.4% R331R (CGC→CGT 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,9160CT48.4% 1.7 / 54.5 64R331R (CGC→CGTribDfused diaminohydroxyphosphoribosylaminopyrimidin e deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base C (10/23);  new base T (13/18);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATAC  >  minE/262853‑262983
                                                               |                                                                   
ggCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 <  1:606549/68‑1 (MQ=255)
ggCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 <  1:681237/68‑1 (MQ=255)
ggCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 <  1:805570/68‑1 (MQ=255)
  cGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTa                                                               >  1:79081/1‑68 (MQ=255)
  cGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTa                                                               >  1:909796/1‑68 (MQ=255)
  cGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTa                                                               >  1:2526757/1‑68 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:2248675/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:681949/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:699120/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:459184/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:815978/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:882929/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:2682532/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:105933/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:1237110/69‑1 (MQ=255)
   gggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:1140111/69‑1 (MQ=255)
    ggTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                             <  1:1549726/68‑1 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 >  1:3405524/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 >  1:651628/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 >  1:2053816/1‑60 (MQ=255)
        tAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                 >  1:1178411/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:2406096/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct                           >  1:2953138/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct                           >  1:2968688/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct                           >  1:180823/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct                           >  1:167040/1‑68 (MQ=255)
                                     ccTAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGct                           >  1:1555087/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)
                                     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:661170/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:59789/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:3120997/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:2600578/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:2838953/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:1466425/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:1961004/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:214493/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:2159928/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:2210993/52‑1 (MQ=255)
                                       tAAACTATTAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:1023297/52‑1 (MQ=255)
                                       tAAACTATTAGGAAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:1924370/52‑1 (MQ=255)
                                           cTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCa                      <  1:2469397/68‑1 (MQ=255)
                                           cTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCa                      <  1:657702/68‑1 (MQ=255)
                                           cTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCa                      <  1:1590534/68‑1 (MQ=255)
                                              ttAGGCAACGACGCCCGTGGATTATGCACGCTGCCAGGGCTTgag                                          <  1:2267179/45‑1 (MQ=255)
                                                  gCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa               <  1:249865/68‑1 (MQ=255)
                                                  gCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa               <  1:854961/68‑1 (MQ=255)
                                                  gCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa               <  1:2933968/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)
                                                           cccGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAgaga     <  1:356625/69‑1 (MQ=255)
                                                              gcgGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTaaa               <  1:940075/56‑1 (MQ=255)
                                                               cgGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATAc  >  1:1028870/1‑68 (MQ=255)
                                                               cgGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATAc  >  1:2639067/1‑68 (MQ=255)
                                                               cgGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATAc  >  1:140197/1‑68 (MQ=255)
                                                               cgGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATAc  >  1:642864/1‑68 (MQ=255)
                                                               |                                                                   
GGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATAC  >  minE/262853‑262983

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: