Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,774,221	G→T	A27D (GCC→GAC)	hemA ←	Glutamyl‑tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,774,221	0	G	T	100.0%	50.4 / NA	18	A27D (GCC→GAC)	hemA	Glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (5/13);  total (5/13)

TTGAATATAGTAACGACCTGTGTGAATTTGATCAACAACAGCATATACTTCAGTTCGATTACATGTTGATAATATGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGGCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAACATAAT  >  USA300TCH1516_ALE/1774091‑1774347                                                                                                                                   |                                                                                                                               acgAATATAGTACCGACCTGTTTGAATTTGATCAACAACAGCATCTACTTCAGTTCGATTACATGTTGATAATATGACATTTTCTAATATAGATTTAGTTTCAAATAAATCTTCCTAGGTAAGTCGTAAGTCATCctctct                                                                                                                      <  2:191369/139‑1 (MQ=255)   gAATATAGTAACGACCTGTGTGAATTTGATCAACAACAGCATATACTTCAGTTCGATTACATGTTGATAATATGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTaa                                                                                                                    >  1:111998/1‑141 (MQ=255)                                         gCATATACTTCAGTTCGATTACATGTTGATAATATGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCga                                                                              >  2:363341/1‑141 (MQ=255)                                                                  ttGATAATATGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAAt                                                     <  2:372433/141‑1 (MQ=255)                                                                          atGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGa                                             <  1:349803/141‑1 (MQ=255)                                                                          atGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGa                                             <  2:111998/141‑1 (MQ=255)                                                                             aCATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAAc                                          <  1:137061/141‑1 (MQ=255)                                                                               aTTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGc                                        <  1:225434/141‑1 (MQ=255)                                                                               aTTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGc                                        <  1:302038/141‑1 (MQ=255)                                                                                                 aGTTTCATATAAATCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACaa                      <  1:272920/141‑1 (MQ=255)                                                                                                            aaTCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAAtatta           >  1:9105/1‑141 (MQ=255)                                                                                                              tCTTCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATtata         <  2:47698/141‑1 (MQ=255)                                                                                                                ttCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAAc       <  2:87862/141‑1 (MQ=255)                                                                                                                ttCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAAc       <  2:4359/141‑1 (MQ=255)                                                                                                                ttCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAAc       <  2:334578/141‑1 (MQ=255)                                                                                                                 tCATGGGCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAACa      <  1:11800/141‑1 (MQ=255)                                                                                                                     gggCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAACATAAt  >  2:202806/1‑141 (MQ=255)                                                                                                                     gggCAATTCGTAAGTCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAAAGCCCCCATATATACTAACAATACAAATATTATAACATAAt  >  1:327405/1‑141 (MQ=255)                                                                                                                                   |                                                                                                                               TTGAATATAGTAACGACCTGTGTGAATTTGATCAACAACAGCATATACTTCAGTTCGATTACATGTTGATAATATGACATTTTCTAAAATAGATTTAGTTTCATATAAATCTTCATGGGCAATTCGTAAGGCATCATCTCTAAAAGCAACTTGCTCTCTTAGTGCAACATCAGCTGTGCGATGATTTATACTAATTGCAATAAAATGCATTTGAAACGCCCCCATATATACTAACAATACAAATATTATAACATAAT  >  USA300TCH1516_ALE/1774091‑1774347

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