Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,178,994 G→T 31.9% T148T (ACG→ACT ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,178,9940GT31.9% 87.9 / 45.0 72T148T (ACG→ACTansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base G (35/14);  new base T (20/3);  total (55/17)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAT  >  minE/1178933‑1179060
                                                             |                                                                  
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:3386954/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:1033717/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:2709760/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:2253287/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:220592/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:3012755/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:2031449/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:3222790/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:3389227/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:3404197/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:427031/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:156604/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:93790/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:613071/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:1296587/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:1337122/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:1436440/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:505302/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:439521/1‑68 (MQ=255)
gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGAAActct                                                                  >  1:2081516/1‑64 (MQ=255)
gACAAATTAATCTGATGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                              >  1:2386042/1‑68 (MQ=255)
  cAAATTAATCTGCTGAATGCGTTGTTCGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                                           >  1:2380877/1‑69 (MQ=255)
  cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                                           >  1:1959763/1‑69 (MQ=255)
  cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                                           >  1:571245/1‑69 (MQ=255)
  cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                                           >  1:499694/1‑69 (MQ=255)
  cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                                           >  1:362395/1‑69 (MQ=255)
  cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                                           >  1:2304872/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:988584/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:2609692/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:2624629/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:2894040/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:2921950/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:3423746/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:3020147/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:3259703/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:2526320/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:1037031/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:1519654/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:2313797/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:1578024/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:1664667/1‑69 (MQ=255)
                   tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                                          >  1:1794697/1‑69 (MQ=255)
                              ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:2710023/68‑1 (MQ=255)
                              ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:2076719/68‑1 (MQ=255)
                              ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:3431809/68‑1 (MQ=255)
                              ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:3410863/68‑1 (MQ=255)
                              ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:1991410/68‑1 (MQ=255)
                              ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:3383709/68‑1 (MQ=255)
                              gtGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCa                               <  1:1175690/68‑1 (MQ=255)
                                     gAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaa                        >  1:2488199/1‑69 (MQ=255)
                                         tATCCTATCAACGAAGTAACTCTCTTTTTTAATAATCGACTGTATCGTGGCAACCGAACCACCAAAGca                    <  1:2435683/69‑2 (MQ=255)
                                         tATCCTATCAACGAAGTAACTCTCTTTTTTAATAATCGACTGTATCGTGGCAACCGAACCACCAAAGca                    <  1:418407/69‑2 (MQ=255)
                                         tATCCTATCAACGAAGTAACTCTCTTTTTTAATAATCGACTGTATCGTGGCAACCGAACCACCAAAGca                    <  1:2005853/69‑2 (MQ=255)
                                         tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                    >  1:1976982/1‑69 (MQ=255)
                                         tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                    <  1:2383754/69‑1 (MQ=255)
                                         tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                    >  1:1260539/1‑69 (MQ=255)
                                         tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                    <  1:2068130/69‑1 (MQ=255)
                                         tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                    <  1:1731508/69‑1 (MQ=255)
                                         tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                    <  1:3404969/69‑1 (MQ=255)
                                             cGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATg                >  1:1802168/1‑69 (MQ=255)
                                                tCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCg             >  1:178651/1‑69 (MQ=255)
                                                tCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCg             >  1:3051717/1‑69 (MQ=255)
                                                tCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCg             >  1:2309243/1‑69 (MQ=255)
                                                tCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCg             >  1:565764/1‑69 (MQ=255)
                                                tCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCg             >  1:3298445/1‑69 (MQ=255)
                                                  aaCGAATTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGAt           <  1:3060223/69‑1 (MQ=255)
                                                  aaCGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGAt           <  1:2385513/69‑1 (MQ=255)
                                                    cGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                <  1:1402008/46‑1 (MQ=255)
                                                           aCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt  >  1:3234812/1‑69 (MQ=255)
                                                           aCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt  >  1:1540304/1‑69 (MQ=255)
                                                           aCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt  >  1:2773442/1‑69 (MQ=255)
                                                           aCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt  >  1:693960/1‑69 (MQ=255)
                                                           aCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt  >  1:88482/1‑69 (MQ=255)
                                                             |                                                                  
GACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAT  >  minE/1178933‑1179060

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: