Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,416,458	A→T	I216N (ATT→AAT)	oppD_4 ←	Oligopeptide transport ATP‑binding protein OppD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,416,458	0	A	T	86.7%	34.3 / ‑5.0	15	I216N (ATT→AAT)	oppD_4	Oligopeptide transport ATP‑binding protein OppD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (8/5);  minor base C (2/0);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGG  >  USA300TCH1516_ALE/1416330‑1416591                                                                                                                                 |                                                                                                                                      aaTCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTAACaa                                                                                                                            >  1:261806/1‑140 (MQ=255)        gTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCa                                                                                                                     >  1:128723/1‑140 (MQ=255)          aaTCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCGTCCTTGAAATGATTAATATTACGTAGTTCACTGATTTGACCTTCTTTTAGAAAGTAACTGTACTGACAATATTTCAAt                                                                                                                   >  2:274918/1‑140 (MQ=255)              ccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTTACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATAAAt                                                                                                               >  1:62840/1‑140 (MQ=255)                                        ttAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                                                                     <  1:130565/140‑1 (MQ=255)                                          aGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATTATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAACTGAACTGACTATATTTCTATACATGACTTAAGTTATGTGTGATAataaataa                                                                                   >  2:132502/1‑140 (MQ=255)                                                                     tCATACTTGAAATGATTAATATTACCTCGTTCAATGTTTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCCTGTTCTAATACATTATg                                                        >  2:73878/1‑140 (MQ=255)                                                                          cTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATTTGTGATAATAAATAATATTTGACCATTTTCTAATAAAATATGCTTta                                                   >  1:145137/1‑140 (MQ=255)                                                                           ttGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGtaa                                                  <  2:46964/140‑1 (MQ=255)                                                                                attaTTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTGCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAAt                                             <  1:132502/137‑1 (MQ=255)                                                                                              ttcttttaaTGATTTGACCTTCTGTTAAAACATAATTGTACTGACAAGATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATATATCCATCACTTGATTa                               <  2:268786/133‑1 (MQ=255)                                                                                                        gATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAg                     >  1:81228/1‑140 (MQ=255)                                                                                                                  tcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAatga           <  2:238236/140‑1 (MQ=255)                                                                                                                         aaaCATAATTGTACTGACACTATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACt    >  1:286609/1‑140 (MQ=255)                                                                                                                           aCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTgg  >  1:100428/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                      AATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGG  >  USA300TCH1516_ALE/1416330‑1416591

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