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%	80.0 / NA	24	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 (9/15);  total (9/15)

TTTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTT  >  USA300TCH1516_ALE/2010226‑2010496                                                                                                                                          |                                                                                                                                      tttCACGTGTCAATTTTCTTAAGCTTCCAGAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAAt                                                                                                                                    <  2:165791/141‑1 (MQ=255)    cACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTa                                                                                                                                 >  2:218206/1‑141 (MQ=255)        cgtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGCTAGTGCAATAATAATTGTTATATAATCTAACCa                                                                                                                             <  1:168416/140‑1 (MQ=255)                   ttAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAACTATTCATTa                                                                                                                  >  2:217135/1‑141 (MQ=255)                       tCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAccc                                                                                                              <  1:224977/141‑1 (MQ=255)                              aaGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaaa                                                                                                       >  2:128588/1‑141 (MQ=255)                              aaGAAAACGTACACCGTTAAAATGTAAAAGGGTAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaaa                                                                                                       >  2:181378/1‑141 (MQ=255)                                 aacaagTACTCCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg                                                                                                    <  2:176160/136‑1 (MQ=255)                                 aaaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg                                                                                                    <  1:263431/141‑1 (MQ=255)                                  aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGa                                                                                                   >  1:292989/1‑141 (MQ=255)                                  aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGa                                                                                                   >  1:39200/1‑141 (MQ=255)                                                  aatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTc                                                                                   <  1:76161/141‑1 (MQ=255)                                                                                 aaaGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                    <  2:39200/141‑1 (MQ=255)                                                                                 aaaGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                    <  2:292989/141‑1 (MQ=255)                                                                                       aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                              <  1:50328/141‑1 (MQ=255)                                                                                        aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACtt                                             <  2:36581/141‑1 (MQ=255)                                                                                                      aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                               <  1:218206/141‑1 (MQ=255)                                                                                                         aaCATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAt                            >  1:146001/1‑141 (MQ=255)                                                                                                            aTTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtctc                         <  1:195182/141‑1 (MQ=255)                                                                                                                  gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGTTCTTTTTGATTATTAGCATAAAATCTCGCACtt                   <  1:128588/141‑1 (MQ=255)                                                                                                                  gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACtt                   <  1:181378/141‑1 (MQ=255)                                                                                                                   aTAGTGCAATAATAATTTTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTa                  <  2:125823/141‑1 (MQ=255)                                                                                                                          aataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCttt           >  2:14899/1‑140 (MQ=255)                                                                                                                                   tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGtt  >  1:2923/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                      TTTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTT  >  USA300TCH1516_ALE/2010226‑2010496

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