Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,444,608 T→A 40.0% V176V (GTA→GTT ilvD ← dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,444,6080TA40.0% 12.3 / 25.3 30V176V (GTA→GTTilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base T (15/3);  new base A (5/7);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 4.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CCGGAGCAGGAACCGCAGGTCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGAAAGTTT  >  minE/2444546‑2444676
                                                              |                                                                    
ccGGAGCAGGAACCGCAGGTCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTc                                                                >  1:902647/1‑69 (MQ=255)
ccGGAGCAGGAACCGCAGGTCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTc                                                                >  1:3282669/1‑69 (MQ=255)
ccGGAGCAGGAACCGCAGGTCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTc                                                                >  1:1682827/1‑69 (MQ=255)
  ggAGCAGGAACCGCAGGTTGGGCAGGCGGATCGTTCAACCTGAGCGCTCTGGGAGTCAGAAACTTTCgg                                                              >  1:1245394/1‑69 (MQ=255)
  ggAGCAGGAACCGCAGGTTGGGCAGGCGGATCGTTCAACCTGAGCGCTCTGGGAGTCAGAAACTTTCgg                                                              >  1:3176722/1‑69 (MQ=255)
  ggAGCAGGAACCGCAGGTTGGGCAGGCGGATCGTTCAACCTGAGCGCTCTGGGAGTCAGAAACTTTCgg                                                              >  1:271683/1‑69 (MQ=255)
  ggAGCAGGAACCGCAGGTTGGGCAGGCGGATCGTTCAACCTGAGCGCTCTGGGAGTCAGAAACTTTCgg                                                              >  1:2416727/1‑69 (MQ=255)
  ggAGCAGGAACCGCAGGTTGGGCAGGCGGATCGTTCAACCTGAGCGCTCTGGGAGTCAGAAACTTTCgg                                                              >  1:2265696/1‑69 (MQ=255)
                   tCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTc                                                           >  1:1452729/1‑55 (MQ=255)
                   tCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGatca                                             >  1:1731877/1‑69 (MQ=255)
                   tCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGatca                                             >  1:43439/1‑69 (MQ=255)
                   tCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGatca                                             >  1:3399657/1‑69 (MQ=255)
                   tCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGatca                                             >  1:2602619/1‑69 (MQ=255)
                   tCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGatca                                             >  1:2553748/1‑69 (MQ=255)
                                              cGCTCTGGGAGTCAGATACTTTCTGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtgatg                   <  1:2282724/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:122517/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:879034/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:836442/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:1445280/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:241124/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:2063240/68‑1 (MQ=255)
                                              cGCTCTGGGAGTCAGAAACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTtaatg                   <  1:1863988/68‑1 (MQ=255)
                                               gCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTTGATCATCGCATCAACCAGATCGAGCtt                       >  1:2234391/1‑63 (MQ=255)
                                               gCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATc                             >  1:1516307/1‑57 (MQ=255)
                                               gCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCtt                       >  1:1435002/1‑63 (MQ=255)
                                               gCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCtt                       >  1:1275628/1‑63 (MQ=255)
                                               gCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCtt                       >  1:1765266/1‑63 (MQ=255)
                                               gCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCtt                       >  1:884290/1‑63 (MQ=255)
                                                        gTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGa        <  1:2475833/69‑1 (MQ=255)
                                                        gTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGa        <  1:2115481/69‑1 (MQ=255)
                                                              tACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGAAAGttt  >  1:2339127/1‑69 (MQ=255)
                                                              tACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGAAAGttt  >  1:1463482/1‑69 (MQ=255)
                                                              tACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGAAAGttt  >  1:2258561/1‑69 (MQ=255)
                                                              tACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGAAAGttt  >  1:457697/1‑69 (MQ=255)
                                                              |                                                                    
CCGGAGCAGGAACCGCAGGTCGGACACGCGGAACGTTCAACCTGATCGCTCTGGGAGTCAGATACTTTCGGGTCTGCGCCCTGGATCATCGCATCAACCAGATCGAGCTTGATGATCTGATCGGAAAGTTT  >  minE/2444546‑2444676

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: