Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,010,363	C→T	C146Y (TGT→TAT)	USA300HOU_RS10125 ←	Putative multidrug export ATP‑binding/permease protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,010,363	0	C	T	100.0%	38.0 / NA	15	C146Y (TGT→TAT)	USA300HOU_RS10125	Putative multidrug export ATP‑binding/permease protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/7);  total (8/7)

AAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGT  >  USA300TCH1516_ALE/2010259‑2010488                                                                                                         |                                                                                                                              aaaCGTACACCGTTTAAATGTAAAATGGGAAAATAAAAAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTa                                                                                                 >  1:41879/1‑135 (MQ=255)  aaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaa                                                                                                <  1:119874/135‑1 (MQ=255)  aaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaa                                                                                                <  1:127786/135‑1 (MQ=255)    cGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGCTTaaaa                                                                                              >  2:139009/1‑135 (MQ=255)                      aaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTACACCGGTTAAAATGAAATCTTTtgtttttt                                                                            >  2:99040/1‑135 (MQ=255)                                gATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAATTTTGTTAAAATGAAATCTTTTGTTTGTTCAACATCAtt                                                                  <  1:56410/135‑1 (MQ=255)                                      cAGTGTTGCTATAGTGTATTTTACTTGTAAGAAGATGATTATGGATATTGTAATTATAATTGTTATATAATCTAACCAAATATTCATTAACCAGGTTAAAATTAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                            <  2:150481/135‑1 (MQ=255)                                        gTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATAAAATCTTTTGTTTGTTCAACATCATTAAGTActc                                                          >  2:154020/1‑135 (MQ=255)                                               tAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat                                                   <  1:82132/135‑1 (MQ=255)                                               tAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat                                                   <  1:96467/135‑1 (MQ=255)                                                               tCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGa                                   >  1:144341/1‑135 (MQ=255)                                                                cTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATATAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGat                                  >  1:34332/1‑135 (MQ=255)                                                                                              taatTGTTATATAATCTAACCAAATATTCATTAAACCGGTTAAAATGAAATCTTTTGTTTGTTCAACATTATTCATAACTATACATATTACTTGCCTTACTTTATTATTAGCATAAAATATAGCACTTATAACtt    >  2:144284/1‑135 (MQ=255)                                                                                                atTGTTATATAATCTAACCAAATATTCATTAACCCGATTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTTTAGATATTACTTTTCTTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGt  <  1:121724/135‑1 (MQ=255)                                                                                                atTGTTATATAATCTAACCAAATATTCATTAACCCGAATAAAAAAAAATCATTAGTTTGTTCAACATCATTAATAACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGATTGt  >  1:106843/1‑135 (MQ=255)                                                                                                         |                                                                                                                              AAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGT  >  USA300TCH1516_ALE/2010259‑2010488

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