Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 642,198 T→G 84.4% I71L (ATC→CTC)  aspC ← aspartate aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE642,1980TG84.4% 180.2 / 18.0 83I71L (ATC→CTC) aspCaspartate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (5/8);  new base G (56/14);  total (61/22)
Fisher's exact test for biased strand distribution p-value = 4.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCTTT  >  minE/642136‑642259
                                                              |                                                             
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:1131192/69‑1 (MQ=255)
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:2079596/69‑1 (MQ=255)
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:687594/69‑1 (MQ=255)
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:1648340/69‑1 (MQ=255)
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:1538608/69‑1 (MQ=255)
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:212379/69‑1 (MQ=255)
aTTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTc                                                         <  1:2152400/69‑1 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGcc                                                            >  1:596964/1‑65 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:1667668/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:2065710/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:1989020/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:3043022/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:241033/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:2432897/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:1849933/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:3353708/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:3420288/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:534814/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:1615415/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:2467988/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:1565861/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:2270396/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:614015/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:1452937/1‑66 (MQ=255)
 ttGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCg                                                           >  1:2535981/1‑66 (MQ=255)
      cAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATg                                                              >  1:489499/1‑58 (MQ=255)
      cAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTCAAt                                                      >  1:2597295/1‑66 (MQ=255)
      cAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTCAATg                                                     >  1:3130669/1‑67 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGCCCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:2822549/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:2457184/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:3299535/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:3271732/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:3225696/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:402424/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:3043826/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:540782/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1205680/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1273538/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1313204/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1318626/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1327459/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1348801/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:760987/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:69605/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1795399/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1898798/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:1928500/1‑69 (MQ=255)
                   ttACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATGCCGAGGTAATTTTtg                                      >  1:643867/1‑69 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:694718/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:691399/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:48484/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:87604/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:2934682/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:1157993/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:1536974/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:1751041/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:2104154/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:2364927/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:2369000/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:237635/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:2486533/48‑1 (MQ=255)
                         aaCAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGAGGCCGTCAATg                                                     <  1:2795700/48‑1 (MQ=255)
                                      gAGTGCAGCGACCAAATTCAGGGATGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCa                   <  1:866182/69‑1 (MQ=255)
                                           cAGCGACCAAATTCAGGGATGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATAc              >  1:83229/1‑69 (MQ=255)
                                           cAGCGAACAAATTCAGGGATGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATAc              >  1:6764/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATa               >  1:526816/1‑56 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:3391714/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:829414/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:1448689/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:318195/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:1602905/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:3220534/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:1973732/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:2084519/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:603053/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:2235705/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:2271693/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:2305381/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:500316/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:2372077/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCttt  >  1:2451312/1‑69 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCtt   >  1:3044027/1‑68 (MQ=255)
                                                       tCAGGGAGGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCtt   >  1:2796290/1‑68 (MQ=255)
                                                              |                                                             
ATTGATCAGGGCGCTACCTTTACCAAACAGCAGTTCCTGAGTGCAGCGACCAAATTCAGGGATGCCGTCAATGCCGAGGTAATTTTTGGTGGTTTCATTTTCGAGCAGATACTGTTCAGCCTTT  >  minE/642136‑642259

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