Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,744,362	G→A	100%	H104Y (CAC→TAC)	apt ←	Adenine phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,744,362	0	G	A	100.0%	94.9 / NA	32	H104Y (CAC→TAC)	apt	Adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/19);  total (13/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAA  >  USA300TCH1516_ALE/1744225‑1744494                                                                                                                                          |                                                                                                                                     ttcaattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACAt                                                                                                                                   <  1:812971/141‑1 (MQ=255)                 aTACCTACTACGATACCGCCTAATTTTTCACCTTAGTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTa                                                                                                                   >  1:194182/1‑141 (MQ=255)                   aCCTACTACGATACCGCCTAATTTTTCAACTTACCCAGCAGCTACTTAAATTCTCCCACCATTAGCTAATAAATCATCTGTAATTAACACACCTTGACCTGGTTTAATTGCATCTTGGTACATTCTTAAAACAGTTGTGcc                                                                                                                 >  1:802187/1‑141 (MQ=255)                   aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAcc                                                                                                                 >  1:800341/1‑141 (MQ=255)                           cGATACCTCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTa                                                                                                         <  1:684102/141‑1 (MQ=255)                              tACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGt                                                                                                      >  1:481417/1‑141 (MQ=255)                                 cGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTCGGTCaa                                                                                                   >  1:302783/1‑140 (MQ=255)                                 cGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCAt                                                                                                   >  1:428672/1‑141 (MQ=255)                                              aaCTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAAt                                                                                      <  1:18682/141‑1 (MQ=255)                                                          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                                                          <  1:442471/141‑1 (MQ=255)                                                          tgcggcCTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                                                          <  1:770777/141‑1 (MQ=255)                                                              gctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAt                                                                      >  1:614983/1‑141 (MQ=255)                                                                ttCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAttt                                                                    <  1:440967/141‑1 (MQ=255)                                                                 tCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTc                                                                   <  1:712379/141‑1 (MQ=255)                                                                 tCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTc                                                                   <  1:262967/141‑1 (MQ=255)                                                                 tCAATCGTACCACCAGGAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTc                                                                   <  1:100569/141‑1 (MQ=255)                                                                  cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                                                                  >  1:50336/1‑141 (MQ=255)                                                                   aaTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTccc                                                                 <  1:48249/141‑1 (MQ=255)                                                                        tACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAATATTTGTACCATATTCTAGGTCATACTCATTACGAATGACTTCACGAGGTAATTTCCcttctt                                                            <  1:584473/141‑1 (MQ=255)                                                                        tACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttctt                                                            <  1:520912/141‑1 (MQ=255)                                                                         accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                           >  1:32650/1‑141 (MQ=255)                                                                           caccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTc                                                         <  1:643762/141‑1 (MQ=255)                                                                               aGTAGCTAATAAATCATCGTTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc                                                     <  1:183396/141‑1 (MQ=255)                                                                                         aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaaaa                                           >  1:610520/1‑140 (MQ=255)                                                                                            tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCa                                        <  1:185988/141‑1 (MQ=255)                                                                                                        aaCACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGaa                            <  1:386651/141‑1 (MQ=255)                                                                                                         acacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAACCCAATCCCCATTGAAt                           >  1:760428/1‑141 (MQ=255)                                                                                                            cacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTTTAACAAGTGCAAAGCCAATCCTCATTGAATAAg                        >  1:156617/1‑141 (MQ=255)                                                                                                            cacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAg                        >  1:329267/1‑141 (MQ=255)                                                                                                             gcGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTCCTTTTCTAACAGGTGCAAAGCCAATCCCCAGTGAATAAGc                       <  1:545743/140‑1 (MQ=255)                                                                                                                      tGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGCACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTACCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCGTTGAATAAGCTACAGGACa              <  1:529029/141‑1 (MQ=255)                                                                                                                                  aTCTTTGTACATTGTTACAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa  <  1:220688/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     TTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAA  >  USA300TCH1516_ALE/1744225‑1744494

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