Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,178,979 G→T 32.8% P143P (CCG→CCT ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,178,9790GT32.8% 84.9 / 46.6 73P143P (CCG→CCTansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base T (21/3);  total (47/26)
Fisher's exact test for biased strand distribution p-value = 4.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATG  >  minE/1178916‑1179046
                                                               |                                                                   
tgAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                                  <  1:1121446/67‑1 (MQ=255)
tgAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                                  <  1:311049/67‑1 (MQ=255)
tgAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                                  <  1:125922/67‑1 (MQ=255)
tgAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                                  <  1:2186482/67‑1 (MQ=255)
tgAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                                  <  1:1611877/67‑1 (MQ=255)
tgAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                                  <  1:145727/67‑1 (MQ=255)
     tACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTATGTTGCGGCGAATTATCCGATCAACGAAg                                                           <  1:1252151/69‑1 (MQ=255)
     tACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                           <  1:1399084/69‑1 (MQ=255)
     tACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                           <  1:846233/69‑1 (MQ=255)
      aCGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGa                                                             <  1:87479/66‑1 (MQ=255)
      aCGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGa                                                             <  1:1373964/66‑1 (MQ=255)
                ggACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGTTCAACGa                                                             <  1:1254928/56‑1 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:439521/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:3386954/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:3222790/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:427031/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:2031449/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:3404197/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:3389227/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:220592/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:2253287/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:3012755/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:505302/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:2709760/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:613071/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:156604/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:1436440/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:1337122/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:1296587/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:93790/1‑68 (MQ=255)
                 gACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCTATCAACGAAGTAACTCTCttt                                                >  1:1033717/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:2304872/1‑69 (MQ=255)
                   cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                             >  1:362395/1‑69 (MQ=255)
                   cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                             >  1:499694/1‑69 (MQ=255)
                   cAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                             >  1:571245/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAAc                                                       >  1:1883083/1‑42 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:3423746/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:3259703/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:988584/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:3020147/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:2921950/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:2894040/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:2624629/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:1037031/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:1519654/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)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:2313797/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:2609692/1‑69 (MQ=255)
                                    tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcg                            >  1:2526320/1‑69 (MQ=255)
                                               ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                  <  1:2710023/68‑1 (MQ=255)
                                               ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                  <  1:2076719/68‑1 (MQ=255)
                                               ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                  <  1:3410863/68‑1 (MQ=255)
                                               ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                  <  1:3431809/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:2005853/69‑2 (MQ=255)
                                                          tATCCTATCAACGAAGTAACTCTCTTTTTTAATAATCGACTGTATCGTGGCAACCGAACCACCAAAGca      <  1:418407/69‑2 (MQ=255)
                                                          tATCCTATCAACGAAGTAACTCTCTTTTTTAATAATCGACTGTATCGTGGCAACCGAACCACCAAAGca      <  1:2435683/69‑2 (MQ=255)
                                                          tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc      >  1:1976982/1‑69 (MQ=255)
                                                          tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc      >  1:1260539/1‑69 (MQ=255)
                                                          tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc      <  1:1731508/69‑1 (MQ=255)
                                                          tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc      <  1:3404969/69‑1 (MQ=255)
                                                          tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc      <  1:2068130/69‑1 (MQ=255)
                                                          tATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc      <  1:2383754/69‑1 (MQ=255)
                                                              cGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATg  >  1:1802168/1‑69 (MQ=255)
                                                               |                                                                   
TGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATG  >  minE/1178916‑1179046

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: