Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	629,681	C→A	L84I (CTA→ATA)	USA300HOU_RS02990 →	3‑hexulose‑6‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	629,681	0	C	A	100.0%	65.5 / NA	18	L84I (CTA→ATA)	USA300HOU_RS02990	3‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/5);  total (13/5)

GTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATT  >  USA300TCH1516_ALE/629550‑629824                                                                                                                                    |                                                                                                                                                gTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGAt                                                                                                                                >  1:184103/1‑149 (MQ=255)          cagTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAtt                                                                                                                       >  1:31511/1‑149 (MQ=255)          cagTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAtt                                                                                                                       >  1:204814/1‑149 (MQ=255)                         cAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaa                                                                                                        >  1:84484/1‑149 (MQ=255)                          aGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaag                                                                                                       >  1:36145/1‑149 (MQ=255)                            acaacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaagaa                                                                                                     >  1:275840/1‑149 (MQ=255)                               acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCt                                                                                                  <  1:251027/149‑1 (MQ=255)                                caTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTc                                                                                                 <  1:316714/149‑1 (MQ=255)                                          gTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAaataa                                                                                       >  1:39189/1‑149 (MQ=255)                                                   ttAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAAtt                                                                              >  1:75117/1‑149 (MQ=255)                                                               aaaaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGatat                                                                  >  1:206192/1‑149 (MQ=255)                                                                                      tGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGCAGCTAATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCACGATTTAGaaa                                           >  1:287607/1‑149 (MQ=255)                                                                                      tGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGaaa                                           >  1:339201/1‑149 (MQ=255)                                                                                                   gCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGa                              <  1:118204/149‑1 (MQ=255)                                                                                                      aTTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACt                           >  1:27193/1‑149 (MQ=255)                                                                                                                      aTGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCt           <  1:56923/149‑1 (MQ=255)                                                                                                                        gTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGa         <  1:121446/149‑1 (MQ=255)                                                                                                                               cAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTAAGGCATTCGTAGCATGTTATGACGAAAGAGATGAAATTGGTTCTGATTATAtt  >  1:65846/1‑149 (MQ=255)                                                                                                                                    |                                                                                                                                                GTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATT  >  USA300TCH1516_ALE/629550‑629824

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