Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,612,073 C→T 44.1% Q360Q (CAG→CAA eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,612,0730CT44.1% 16.6 / 53.0 71Q360Q (CAG→CAAeutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base C (9/30);  new base T (10/21);  total (19/52)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCA  >  minE/1612007‑1612137
                                                                  |                                                                
gTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTg                                                                <  1:666643/69‑1 (MQ=255)
 tAACGCATAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:2404313/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:9830/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:2579801/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:1307234/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:1856232/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:2143613/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:163095/69‑1 (MQ=255)
 tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                               <  1:413125/69‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:3049725/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:257595/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2179688/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2968649/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2809380/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2766479/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2598080/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2425130/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2587846/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:3083994/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:3243509/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:1621681/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:1520207/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:1462152/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:446669/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:472156/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:474538/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:515934/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:1265533/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:840527/67‑1 (MQ=255)
   aCGCACAGGCAGCGAGCCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                               <  1:2279345/67‑1 (MQ=255)
              gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:3387491/69‑1 (MQ=255)
              gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:2975654/69‑1 (MQ=255)
              gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:2817081/69‑1 (MQ=255)
              gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:985722/69‑1 (MQ=255)
              gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:2489156/69‑1 (MQ=255)
              gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:2204798/69‑1 (MQ=255)
               cGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:2284128/68‑1 (MQ=255)
                 aGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCTTGTTGc                                                <  1:4499/68‑1 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGcc                                               >  1:2777597/1‑68 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:1166626/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:610392/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:1370616/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:202669/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:695394/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:219918/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:2949640/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:886534/1‑69 (MQ=255)
                  ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                              >  1:1389468/1‑69 (MQ=255)
                         aGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccc                                       >  1:3151785/1‑69 (MQ=255)
                         aGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGcccc                                        >  1:1386367/1‑68 (MQ=255)
                                gCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                                  <  1:2155431/51‑1 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCCGCGCCTGTTGCCAGGCCCCCACCAGAt                                >  1:2575353/1‑56 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt                   >  1:491111/1‑69 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt                   >  1:40910/1‑69 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt                   >  1:3398330/1‑69 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt                   >  1:2669320/1‑69 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt                   >  1:1098557/1‑69 (MQ=255)
                                             aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt                   >  1:2494811/1‑69 (MQ=255)
                                                 ttGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCggg               <  1:2420082/69‑1 (MQ=255)
                                                 ttGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCggg               <  1:1389875/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:607864/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:1043488/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:1602166/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:1821658/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:1909804/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:2207739/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:2635437/69‑1 (MQ=255)
                                                              ccAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:244984/69‑1 (MQ=255)
                                                              ccAGATGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:2714844/69‑1 (MQ=255)
                                                               cAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:1449531/68‑1 (MQ=255)
                                                               cAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa  <  1:2128642/68‑1 (MQ=255)
                                                                  |                                                                
GTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCA  >  minE/1612007‑1612137

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: