Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,515,086 G→A 74.5% I135I (ATC→ATT purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,515,0860GA74.5% 97.8 / 30.3 82I135I (ATC→ATTpurFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (16/5);  new base A (6/55);  total (22/60)
Fisher's exact test for biased strand distribution p-value = 2.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-03

ATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTT  >  minE/1515021‑1515135
                                                                 |                                                 
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGGTTGCGAAAata                                                <  1:2650322/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1027692/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1831299/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:2918623/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1893354/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1894572/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:3344128/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:194461/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:3386327/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:341039/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:360076/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1406928/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1375410/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1373223/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:532387/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:818529/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1260100/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:2537015/69‑1 (MQ=255)
aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1234445/69‑1 (MQ=255)
  ggcggcGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:3091806/67‑1 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:2091478/1‑53 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:2771308/1‑53 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:2332282/1‑53 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:3145195/1‑53 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:1853129/1‑53 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:918423/1‑53 (MQ=255)
              atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa                                                  >  1:1656683/1‑53 (MQ=255)
                      cTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:2092724/47‑1 (MQ=25)
                               ggTACTGGCGGAAGTTGTCCAGTTCTCTTGCGAAAata                                                <  1:84760/38‑1 (MQ=16)
                                  aGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              <  1:176064/69‑1 (MQ=255)
                                  aGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              <  1:365576/69‑1 (MQ=255)
                                  aGTGGCGGAAGGTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              <  1:495147/69‑1 (MQ=255)
                                   gTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              <  1:2142539/68‑1 (MQ=255)
                                   gTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              <  1:2167487/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2344193/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3066869/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:24598/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2282202/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2561580/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:126527/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1626131/68‑1 (MQ=255)
                                     ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1885209/68‑1 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAgt                      >  1:10641/1‑57 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              >  1:1239033/1‑65 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              >  1:1271811/1‑65 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              >  1:1566859/1‑65 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa              >  1:2983261/1‑65 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg          >  1:2544573/1‑69 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg          >  1:1428367/1‑69 (MQ=255)
                                      gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg          >  1:236374/1‑69 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1804484/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:964651/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1171863/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1264711/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1504802/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3299601/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3229136/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3173778/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1846310/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2567625/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2448392/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2203646/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1882832/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2191977/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2172000/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2888177/65‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc        >  1:164612/1‑69 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc        >  1:25598/1‑69 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc        >  1:597161/1‑69 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:398712/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:3040862/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:2319720/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:2220680/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:783440/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:1371423/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         >  1:2565512/1‑68 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         >  1:2813559/1‑68 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:2567770/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:1803357/68‑1 (MQ=255)
                                        ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         >  1:2717476/1‑68 (MQ=255)
                                              tGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGtt  >  1:3154759/1‑69 (MQ=255)
                                                      tCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:1339525/54‑1 (MQ=255)
                                                                 |                                                 
ATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTT  >  minE/1515021‑1515135

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: