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%	86.2 / NA	25	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/17);  total (8/17)

GTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGTATA  >  USA300TCH1516_ALE/2010232‑2010501                                                                                                                                    |                                                                                                                                           gtttctaTTTTCTTTATCTTTCATAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATAtt                                                                                                                            <  1:45560/142‑1 (MQ=255)  tgtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTc                                                                                                                           <  1:50427/148‑1 (MQ=255)   gttaaTTGTCTTAATCTGCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCa                                                                                                                          <  1:615306/145‑1 (MQ=255)   gtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCa                                                                                                                          >  1:369400/1‑148 (MQ=255)            cTTAATCTTCCAAAGAAAACGTACAACGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGCTAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCgg                                                                                                                 <  1:624355/148‑1 (MQ=255)                                           aaatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCACTTTCACAATTACATATGTGATTATGGATATTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGCTTGTTCAAcatc                                                                                  <  1:431866/148‑1 (MQ=255)                                                     tgGGAAGATAACCCGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAct                                                                        <  1:212055/148‑1 (MQ=255)                                                         aaGATAAACAGTGCTGCAGAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTATGTTCAACATCATTAATCACTCTat                                                                    >  1:464503/1‑147 (MQ=255)                                                             tAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATAACTCTAGatat                                                                >  1:511934/1‑148 (MQ=255)                                                               aaCAGTGCTGCTAAAGTCAAGTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTa                                                              >  1:403898/1‑148 (MQ=255)                                                                    tgctgcTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                         <  1:614682/148‑1 (MQ=255)                                                                    tgctgcTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                         <  1:579373/148‑1 (MQ=255)                                                                      ctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGGcc                                                       >  1:188728/1‑148 (MQ=255)                                                                      ctgctAAAGTCAATGTCACATCTAAAAAGAACATTATGGATAGGGCAATAATAATTGTTATATAATCTGACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACGTCATTAATCACTCTAGATATTACTTGAcc                                                       <  1:442232/148‑1 (MQ=255)                                                                        gctAAAGTAAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTa                                                     <  1:44316/148‑1 (MQ=255)                                                                            aaGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTg                                                 <  1:543611/148‑1 (MQ=255)                                                                             aGTCAATTCCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGa                                                <  1:333402/148‑1 (MQ=255)                                                                                 aaTTTCACATCTCAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGattat                                            <  1:554966/148‑1 (MQ=255)                                                                                 aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGattat                                            <  1:569015/148‑1 (MQ=255)                                                                                    ttCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAg                                         <  1:612451/148‑1 (MQ=255)                                                                                                   aaCATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCAc                          <  1:340141/148‑1 (MQ=255)                                                                                                            gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCtt                 >  1:496014/1‑148 (MQ=255)                                                                                                                  gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGGTAAAAGGAAATCTTTTGTTTGTTCAACATCATTAATCACACTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAATATTGTAAAt           >  1:477994/1‑148 (MQ=255)                                                                                                                         tagttGTTATATAATCTAATCAAATATTCATTAACCCGGTTAATATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGta    <  1:183809/145‑1 (MQ=255)                                                                                                                           atTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGtata  >  1:29593/1‑148 (MQ=255)                                                                                                                                    |                                                                                                                                           GTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGTATA  >  USA300TCH1516_ALE/2010232‑2010501

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