Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,018,582	T→A	K113* (AAA→TAA)	mgt ←	Monofunctional glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,018,582	0	T	A	100.0%	59.1 / NA	22	K113* (AAA→TAA)	mgt	Monofunctional glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (10/12);  total (10/12)

TTTACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCC  >  USA300TCH1516_ALE/2018455‑2018709                                                                                                                                |                                                                                                                                tttACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTATTCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCg                                                                                                                          <  1:172693/135‑1 (MQ=255)            aGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCatgatg                                                                                                               <  1:104055/135‑1 (MQ=255)                  tGAACGATCATTATAATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACAGCTCTGTCTCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTa                                                                                                         <  2:149729/135‑1 (MQ=255)                    aaCGATCATTATCATAAAAATAATTTTTGACAATTTGTTTTGTAATGGTACTACCACGTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGa                                                                                                       <  2:197888/135‑1 (MQ=255)                    aaCGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGa                                                                                                       >  1:136299/1‑135 (MQ=255)                                       ataatTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTACACATAACAGACACTAATAGTTGAAAATAAAGATCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTACAATCGTTCATCTTCCATTGa                                                                                    >  1:41839/1‑135 (MQ=255)                                          atTTTTTACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAACAAAGCTCTAGTTGTACCTTACAAATTTATTCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAAt                                                                                 <  1:106064/135‑1 (MQ=255)                                                     cTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATTATTGTAGATTCGTACATCTTTCTATGAATTAAAGGcccgtt                                                                      >  2:178706/1‑130 (MQ=255)                                                          tgtgTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACa                                                                 >  2:177008/1‑135 (MQ=255)                                                            tgtAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATa                                                               <  2:136299/135‑1 (MQ=255)                                                               aaTGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATTTTCCATTGAAATAAAGGCACCTTTAACATActc                                                            >  2:1766/1‑135 (MQ=255)                                                                          caccTTGCACATCTCTGTCGCTAATCGTTTAAAATAAAGCTCTAGTTGTACCTTACAAATCGATTTCATGATGATTGTAGAATCGTTCATTTTCTATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTAt                                                 <  1:50368/135‑1 (MQ=255)                                                                           accTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCTATGATGATTGTATAATCGTTTATTTTTCATTTAAATAAAGGCACCTTTAACATACTCTGGCATGTTATc                                                <  1:36021/135‑1 (MQ=255)                                                                              ttGCACATCTATGTCGCTAATTGTTGAAAATAAAGCTCTAGTTGTACATTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACGTTGAACATACTTTGGCATGTTATCAGc                                             <  2:36277/135‑1 (MQ=255)                                                                                          gTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTTCATTGAAATAAAGGCTCCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAact                                 <  1:164352/135‑1 (MQ=255)                                                                                          gTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAact                                 >  2:114153/1‑135 (MQ=255)                                                                                               tAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACttt                            >  1:194442/1‑135 (MQ=255)                                                                                                               aGCTCTAGTTGTACCTGACAAATCGATTTTCTGATGAGTGAAGAAACGTTCATCTTCCATTGAAATAAAGACACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCt            >  1:202660/1‑135 (MQ=255)                                                                                                                   atAGTTCTACCTTATAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTCATTTTCAATTTTTCTTAGt        >  2:25553/2‑135 (MQ=255)                                                                                                                    tAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAATCTACTTTTATTTTCAATTTTTCTTAGtt       <  1:1766/135‑1 (MQ=255)                                                                                                                        tGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACACAAGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATc   >  2:86054/1‑135 (MQ=255)                                                                                                                         gTATCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATTAAGGCACTTTTAACTTTCTCTAGTATTTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATcc  <  1:46313/135‑1 (MQ=255)                                                                                                                                |                                                                                                                                TTTACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCC  >  USA300TCH1516_ALE/2018455‑2018709

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