Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,609,750 A→T 28.5% V202V (GTT→GTA eutC ← ethanolamine ammonia‑lyase, small subunit (light chain)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,609,7500AT28.5% 79.5 / 36.1 63V202V (GTT→GTAeutCethanolamine ammonia‑lyase, small subunit (light chain)
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base T (3/15);  total (21/42)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGT  >  minE/1609688‑1609807
                                                              |                                                         
aTCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:884112/69‑1 (MQ=255)
aTCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:1039312/69‑1 (MQ=255)
aTCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:1061427/69‑1 (MQ=255)
aTCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:2936908/69‑1 (MQ=255)
aTCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:1360173/69‑1 (MQ=255)
      actTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:2216927/63‑1 (MQ=255)
      actTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:572778/63‑1 (MQ=255)
      actTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     <  1:2136849/63‑1 (MQ=255)
      actTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAg                                                     >  1:3294398/1‑63 (MQ=255)
                       ctcACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGccc                              >  1:2634039/1‑69 (MQ=255)
                       ctcACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGccc                              >  1:2700588/1‑69 (MQ=255)
                       ctcACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGccc                              >  1:2027748/1‑69 (MQ=255)
                       ctcACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGccc                              >  1:1809663/1‑69 (MQ=255)
                          aCCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGc                           >  1:2479963/1‑69 (MQ=255)
                          aCCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGc                           >  1:2281939/1‑69 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1544919/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1718466/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1492850/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1463393/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:239113/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:882118/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1362597/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:2964260/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1273264/69‑1 (MQ=255)
                            cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:820452/69‑1 (MQ=255)
                             gatctgatctTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:396515/68‑1 (MQ=255)
                                   atcttcaatcttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAg                  >  1:40782/1‑69 (MQ=255)
                                   atcttcaatcttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAg                  >  1:2408053/1‑69 (MQ=255)
                                   atcttcaatcttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAg                  >  1:510371/1‑69 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                       <  1:3065587/62‑1 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                       <  1:374880/62‑1 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                       <  1:1705113/62‑1 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                       <  1:2532493/62‑1 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                       <  1:2610546/62‑1 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAgg                 >  1:919867/1‑68 (MQ=255)
                                     cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGc                >  1:3042205/1‑69 (MQ=255)
                                      ttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                       <  1:2937952/61‑1 (MQ=255)
                                        caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTc                    >  1:2121026/1‑62 (MQ=255)
                                        caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGccc              >  1:628993/1‑68 (MQ=255)
                                        caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGccc              >  1:2182979/1‑68 (MQ=255)
                                        caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGccc              >  1:1259770/1‑68 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:3244234/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:681650/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:657807/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:1297474/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:1559236/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:1175814/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:366460/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:2094780/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:3122675/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:3029553/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:2274686/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:2785520/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:1001378/69‑1 (MQ=255)
                                            cttcaCGCGACCATAACATACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:2685841/69‑1 (MQ=255)
                                             ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc        >  1:1926970/1‑69 (MQ=255)
                                             ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc        >  1:1030650/1‑69 (MQ=255)
                                             ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc        >  1:680004/1‑69 (MQ=255)
                                             ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc        >  1:2290453/1‑69 (MQ=255)
                                              tcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt         <  1:3042524/67‑1 (MQ=255)
                                                       cATAACGAACAAAGAATGGTGTACCACTTTCAGCCCGGCCTGt                       <  1:1458480/43‑1 (MQ=38)
                                                         tAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                         <  1:1638221/40‑1 (MQ=255)
                                                           aCGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGt  >  1:834728/1‑61 (MQ=255)
                                                              |                                                         
ATCACTACTTTCGCGCCGAGGATCTCACCGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGT  >  minE/1609688‑1609807

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: