Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,761,889 C→A 100% D388Y (GAT→TAT)  fpgS ← Folylpolyglutamate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,761,8890CA100.0% 73.4 / NA 25D388Y (GAT→TAT) fpgSFolylpolyglutamate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/13);  total (12/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGT  >  USA300TCH1516_ALE/1761753‑1762022
                                                                                                                                        |                                                                                                                                     
tcttcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATatac                                                                                                                                    <  1:260459/140‑1 (MQ=255)
   tcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAAAGTAGTCATATACTaa                                                                                                                                 >  1:558211/1‑140 (MQ=255)
     aTGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCAACTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAAtt                                                                                                                               <  1:25681/140‑1 (MQ=255)
                         tttaTAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCTTCATTTTTAAta                                                                                                           <  1:540333/140‑1 (MQ=255)
                           tataATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAAtact                                                                                                         >  1:480013/1‑140 (MQ=255)
                                            tAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTcttctt                                                                                        <  1:471742/140‑1 (MQ=255)
                                                       tAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTAATAGCAAGTTCTTCTTTGCGTAACGCt                                                                             <  1:182702/140‑1 (MQ=255)
                                                        aaCTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCtt                                                                            <  1:591948/140‑1 (MQ=255)
                                                            tCACTTATGAAATACAGACTATCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCgg                                                                        <  1:446128/140‑1 (MQ=255)
                                                                ttATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTCTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGaaaa                                                                    >  1:718260/1‑140 (MQ=255)
                                                                 tATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAAt                                                                   >  1:473561/1‑140 (MQ=255)
                                                                             aCTACCGGTAATTAGTAATCTATCAGCTTGATAATTTTTTATAAGTTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTg                                                       >  1:418751/1‑140 (MQ=255)
                                                                              cTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATAAAATTCAGTGGc                                                      >  1:668310/1‑140 (MQ=255)
                                                                              cTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATAATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGc                                                      >  1:61846/1‑140 (MQ=255)
                                                                                aCCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATATTATCTATAATTTCTTCTTTGAGTAACGCTTTCGGAAAATCAAATTCAGTTGTAt                                                    >  1:257891/1‑140 (MQ=255)
                                                                                       cttACTAATGTAGCACCTTGATAATTTTTTATAAATTCAACTTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCGTTGCGTAACGCGTTCGGAAAATCAAATTCAGTTGCATAAAACGt                                             <  1:72166/139‑1 (MQ=255)
                                                                                         tACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTAt                                           <  1:37264/140‑1 (MQ=255)
                                                                                             aaggTATCATCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATgcgc                                       <  1:357395/137‑1 (MQ=255)
                                                                                                  aTCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTa                                  >  1:372637/1‑140 (MQ=255)
                                                                                                        ttGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTa                            <  1:118590/140‑1 (MQ=255)
                                                                                                         tGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCGACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaa                           >  1:531317/1‑140 (MQ=255)
                                                                                                          gATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaaa                          >  1:631544/1‑140 (MQ=255)
                                                                                                             aaTTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGt                       <  1:350688/140‑1 (MQ=255)
                                                                                                                        aaTTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAtctc            >  1:679856/1‑140 (MQ=255)
                                                                                                                                  aGTCATATACTATTTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGt  <  1:68069/140‑1 (MQ=255)
                                                                                                                                  aGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGt  <  1:330751/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
TCTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGT  >  USA300TCH1516_ALE/1761753‑1762022

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