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%	87.0 / NA	28	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 (11/17);  total (11/17)

GTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATG  >  USA300TCH1516_ALE/2010232‑2010493                                                                                                                                    |                                                                                                                                   gtgtCAATTTTCTTAATCTTCAAAGGAAAACGTACACGGTTATACTGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTTAATTTCACATCTATTAAGAACATCATGTAAAGTGCCATAATAATTCTTATATAATCTAAcc                                                                                                                           >  2:3357/1‑141 (MQ=255)          ttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCa                                                                                                                  >  1:252816/1‑141 (MQ=255)                           aaaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg                                                                                                 >  1:198417/1‑141 (MQ=255)                              aCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaaa                                                                                              <  2:213299/141‑1 (MQ=255)                              aCGTAAACCGTTAAAATGTAAAATAGTAAGATAAACAGTGCTGCTAAAGTTAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaaa                                                                                              <  2:100261/141‑1 (MQ=255)                                     ccGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttg                                                                                       <  2:285590/141‑1 (MQ=255)                                               gtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAATATTAAATCTTTTGTTTGTTCAAc                                                                             <  1:12299/141‑1 (MQ=255)                                                  aaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatc                                                                          <  1:5948/141‑1 (MQ=255)                                                  aaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGCTTAAAATGAAATCTTTTGTTTGTTCAAcatc                                                                          >  1:219576/1‑141 (MQ=255)                                                               aaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTa                                                             >  2:272401/1‑141 (MQ=255)                                                                       ttctAAAGTTGATTTCACATATAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCTGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTAc                                                     <  1:8195/139‑1 (MQ=255)                                                                            aaGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAc                                                <  1:239218/141‑1 (MQ=255)                                                                              ttCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCt                                              <  2:241759/140‑1 (MQ=255)                                                                                 aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                           <  1:44557/141‑1 (MQ=255)                                                                                 aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                           <  1:88024/141‑1 (MQ=255)                                                                                            tAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTCACCTACTTGAttcttacc                                >  2:236807/1‑139 (MQ=255)                                                                                                  gAACATTATGGATAGTACAATAATAATTGTTATATAATCTAACCAAATATTCATTAACACGATTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaaaa                          <  2:97017/141‑1 (MQ=255)                                                                                                   aaCATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTACATCTTCCTTGACCTCCTTGATTATTAACATAAAAt                         >  1:268966/1‑141 (MQ=255)                                                                                                          tGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCAc                  >  1:71533/1‑141 (MQ=255)                                                                                                            gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACtt                <  2:223399/141‑1 (MQ=255)                                                                                                                  gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCt          >  2:26395/1‑141 (MQ=255)                                                                                                                  gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCt          >  2:267723/1‑141 (MQ=255)                                                                                                                   cAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTTCTTGATTATTAGCATAAAATCTCGCACTTAAAGCtt         <  1:272401/141‑1 (MQ=255)                                                                                                                     ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTCCTTGACCTACTTGATTATTAGCATAAAATCTCTCACTTAAAGCTTGt       <  2:279054/141‑1 (MQ=255)                                                                                                                     ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGt       >  1:69592/1‑141 (MQ=255)                                                                                                                      taataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCCCTTAAAGCTTGTa      <  2:251340/141‑1 (MQ=255)                                                                                                                          aatTTTTATATAATCTATTCTAATATTCATTAACCTGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCCCTCTAGATATTACTTTACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATg  <  2:11008/141‑1 (MQ=255)                                                                                                                          aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATg  <  2:77900/141‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   GTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATG  >  USA300TCH1516_ALE/2010232‑2010493

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