Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,607,132 C→T 100% E275K (GAA→AAA)  frc ← Formyl‑CoA:oxalate CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,607,1320CT100.0% 40.0 / NA 12E275K (GAA→AAA) frcFormyl‑CoA:oxalate CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (2/10);  total (2/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

AAATTAAGTTTAATTAAAGTTCCTCTTCCAGAACGTGCAATCGAATCTCCTGTCCAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTCTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTTTTATAACCTACATGGCTGTCTATTAAAGTATCTAACCTGGGCACACCATCATTAATGT  >  USA300TCH1516_ALE/1607000‑1607263
                                                                                                                                    |                                                                                                                                   
aaaTTAAGTTTAATTAAAGTTCCTCTTCCAGAACGTGCAATCGAATCTCCTGTCCAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCtt                                                                                                                             >  2:135013/1‑141 (MQ=255)
      aGTTTAATTAAAGTTCCTCTTCCAGAACGTGCAATCGAATCTCCTGTCCAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCg                                                                                                                       <  2:79509/141‑1 (MQ=255)
         ttaattaaAGTTCCTCTTCCAGAACGTGCAATCGAATCTCCTGTCCAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCGACt                                                                                                                    <  2:40051/141‑1 (MQ=255)
                          tCCAGAACGTGCAATCGAATCTCCTGTCCAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGGt                                                                                                   <  2:18174/141‑1 (MQ=255)
                                                      cAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTGTTCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTTATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTAt                                                                       <  1:100416/141‑1 (MQ=255)
                                                             ttATTGGAATTCGCCCATCTATATCTAATATTATGTTAATGGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGCTACTGTCTCGTATGTGTTAGCTTTTCTATTTAACtt                                                                <  2:41423/141‑1 (MQ=255)
                                                                ttGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGCTTAGTGGTACAGTATCGTATGTGTTTCCTTTTCTATTTTACTCCCt                                                             >  1:51868/1‑141 (MQ=255)
                                                                           ccatccatATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTCTTATAACCTa                                                  <  1:3407/141‑1 (MQ=255)
                                                                            catccatATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTTTTATAACCTAc                                                 <  1:46010/141‑1 (MQ=255)
                                                                                                                      tCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTTTTATAACCTACATGGCTGTCTATTAAAGTATCTAACCTGGGCACACCATCAtt       <  1:63539/141‑1 (MQ=255)
                                                                                                                      tCCATTTTTAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTTTTATAACCTACATGGCTGTCTATTAAAGTATCTAACCTGGGCACACCATCAtt       <  1:74413/141‑1 (MQ=255)
                                                                                                                           ttttAATTTTTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTTTTATAACCTACATGGCTGTCTATTAAAGTATCTAACCTGGGCACACCATCATTAATGt  <  2:31010/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
AAATTAAGTTTAATTAAAGTTCCTCTTCCAGAACGTGCAATCGAATCTCCTGTCCAAACATTTATTGGAATTCGCCCATCCATATCTAATGTTATGTTAATAGTCTTTTTTACTTTTTTCCATTTTTAATTTCTGTATCTTTTACCGACTTTATTTTGATTAGTGGTACAGTATCGTATGTGTTATCTTTTCTATTTAACTTCCTTTTATAACCTACATGGCTGTCTATTAAAGTATCTAACCTGGGCACACCATCATTAATGT  >  USA300TCH1516_ALE/1607000‑1607263

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