Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,515,074 C→T 66.7% E139E (GAG→GAA purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,515,0740CT66.7% 63.8 / 48.7 90E139E (GAG→GAApurFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/5);  new base T (6/54);  total (31/59)
Fisher's exact test for biased strand distribution p-value = 5.59e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

TGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTT  >  minE/1515011‑1515135
                                                               |                                                             
tgtgGCAGCAATGGCAGCGAAACTATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt                                                          >  1:506017/1‑69 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGc                                                           >  1:1692178/1‑68 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGc                                                           >  1:400334/1‑68 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGc                                                           >  1:2568088/1‑68 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt                                                          >  1:2002063/1‑69 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt                                                          >  1:2717117/1‑69 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt                                                          >  1:3323577/1‑69 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt                                                          >  1:2119379/1‑69 (MQ=255)
tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt                                                          >  1:1070412/1‑69 (MQ=255)
          aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGGTTGCGAAAata                                                <  1:2650322/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:1831299/69‑1 (MQ=255)
          aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:3386327/69‑1 (MQ=255)
          aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:194461/69‑1 (MQ=255)
          aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:341039/69‑1 (MQ=255)
          aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:1027692/69‑1 (MQ=255)
          aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata                                                <  1:2918623/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: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:2091478/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:2561580/68‑1 (MQ=255)
                                               ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:24598/68‑1 (MQ=255)
                                               ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:126527/68‑1 (MQ=255)
                                               ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2344193/68‑1 (MQ=255)
                                               ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1626131/68‑1 (MQ=255)
                                               ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2282202/68‑1 (MQ=255)
                                               ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3066869/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:236374/1‑69 (MQ=255)
                                                gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg          >  1:1428367/1‑69 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3173778/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3229136/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:3299601/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1882832/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1846310/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1804484/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1504802/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1264711/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:1171863/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:964651/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2172000/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2191977/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2203646/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2567625/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2448392/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg            <  1:2888177/65‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc        >  1:25598/1‑69 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc        >  1:164612/1‑69 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc        >  1:597161/1‑69 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         >  1:2565512/1‑68 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:2567770/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:783440/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         >  1:2717476/1‑68 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:1371423/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:2319720/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         >  1:2813559/1‑68 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:398712/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:1803357/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:2220680/68‑1 (MQ=255)
                                                  ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg         <  1:3040862/68‑1 (MQ=255)
                                                        tGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGtt  >  1:3154759/1‑69 (MQ=255)
                                                               |                                                             
TGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTT  >  minE/1515011‑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: