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	96.3%	97.7 / NA	27	I216N (ATT→AAT)	oppD_4	Oligopeptide transport ATP‑binding protein OppD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (12/14);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGAT  >  USA300TCH1516_ALE/1416326‑1416583                                                                                                                                     |                                                                                                                              aCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGa                                                                                                                       <  1:410999/141‑1 (MQ=255)   tCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACa                                                                                                                     <  1:277271/141‑1 (MQ=255)      aTCATAGTAATCCCTATTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAata                                                                                                                  <  1:361342/141‑1 (MQ=255)            gTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCaa                                                                                                            >  2:115047/1‑141 (MQ=255)                 cccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACAt                                                                                                       <  2:367116/141‑1 (MQ=255)                 cccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACAt                                                                                                       <  2:222148/141‑1 (MQ=255)                 cccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACAt                                                                                                       <  2:296712/141‑1 (MQ=255)                  ccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTGCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATg                                                                                                      <  2:106446/141‑1 (MQ=255)                    tctTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGAc                                                                                                    >  2:1958/1‑141 (MQ=255)                     ctTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACt                                                                                                   <  1:430920/141‑1 (MQ=255)                       ttAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTa                                                                                                 <  1:157/141‑1 (MQ=255)                                         tAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGAtaataa                                                                               >  1:228940/1‑141 (MQ=255)                                                aCGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGt                                                                        >  1:97233/1‑141 (MQ=255)                                                   ttCAGTATACGGATGCAAATGCTTAAACTTGATATTGTTAATTTTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTg                                                                     <  1:126318/141‑1 (MQ=255)                                                     cAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGAc                                                                   <  1:162788/141‑1 (MQ=255)                                                     cAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGAc                                                                   <  2:181796/141‑1 (MQ=255)                                                                                     tGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATcc                                   <  2:301959/141‑1 (MQ=255)                                                                                      gATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCa                                  >  1:347598/1‑141 (MQ=255)                                                                                                  cTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTAc                      >  1:214158/1‑141 (MQ=255)                                                                                                          aTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCaaa              >  1:44728/1‑141 (MQ=255)                                                                                                               ttGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATc         <  1:1958/141‑1 (MQ=255)                                                                                                                  aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCaa      >  1:145293/1‑141 (MQ=255)                                                                                                                  aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCaa      >  1:54181/1‑141 (MQ=255)                                                                                                                  aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCaa      >  1:365828/1‑141 (MQ=255)                                                                                                                  aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCaa      >  1:184522/1‑141 (MQ=255)                                                                                                                  aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTATAATACAATATGCTGTAATAAAACCATCACTTGATTAACGTTCAAAACATCCaa      >  1:39380/1‑141 (MQ=255)                                                                                                                      tcaaaaaaaaCATAAATAAACTAAAAATATTTCAATAAATAACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAAAAAAAAAGATGTACTAAACCCATCCCTTGCTTACCGTTCCACGCATCCAatgat  <  2:434253/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                              ACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGAT  >  USA300TCH1516_ALE/1416326‑1416583

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: ‑