Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	638,614	G→T	100%	T14K (ACA→AAA)	pdxK ←	Pyridoxine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	638,614	0	G	T	100.0%	59.7 / NA	21	T14K (ACA→AAA)	pdxK	Pyridoxine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/8);  total (13/8)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  USA300TCH1516_ALE/638493‑638749                                                                                                                          |                                                                                                                                        tgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATAACTGCGCCAACACTTTTGTCAGAACCGGCAATTg                                                                                                                       >  1:180215/1‑140 (MQ=255)     tcttCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaa                                                                                                                   <  1:28362/137‑1 (MQ=255)                    gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGc                                                                                                    >  1:183395/1‑140 (MQ=255)                         tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACTTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTg                                                                                               <  1:16396/140‑1 (MQ=255)                                         aTGGCAGTTAAAGCGACCATGGCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTATTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACat                                                                               >  1:150245/1‑140 (MQ=255)                                             cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataata                                                                           >  2:183530/1‑140 (MQ=255)                                             cAGTTAAAGCGAACATGACATACGTATCTAATTCATCGAACGTATAGAAATCAGCGAATATACCGGCTACAGCACTTTTTTCTGATCCGGCTATTTTTAAAACTTTCGTTAAAGCCATTGAGCTTCACTCCTACataata                                                                           <  2:287370/140‑1 (MQ=255)                                              aGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                          >  1:259210/1‑140 (MQ=255)                                                    aGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGt                                                                    >  1:221416/1‑140 (MQ=255)                                                     gCGACAATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTCAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                                                                   <  1:185966/140‑1 (MQ=255)                                                              ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCCGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                          >  1:143608/1‑140 (MQ=255)                                                                        cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                >  1:146539/1‑140 (MQ=255)                                                                                 ggAACGATTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCGTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                       >  2:160413/1‑140 (MQ=255)                                                                                    aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaa                                    >  2:275161/1‑140 (MQ=255)                                                                                       tctgcAAATCTACTTGCATACCGGCTCCTGCACTTTTGTCATAACCTTCTATTGTTAAAAGCTTCTTTAAATCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTTATTGTAAtata                                 <  2:254397/136‑1 (MQ=255)                                                                                                cTGCTTGCATACCTGCGCAAGCACGTTTGTCACAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTAACAAATATTAAGAAtt                        >  2:33193/1‑140 (MQ=255)                                                                                                 tGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTc                       <  2:36953/140‑1 (MQ=255)                                                                                                       cATACCTGCGCCAGCACTTTTGTCAGAGCCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATAtt                 <  1:164214/140‑1 (MQ=255)                                                                                                                gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg        >  2:99226/1‑140 (MQ=255)                                                                                                                 ccAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCACTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt       <  1:177824/140‑1 (MQ=255)                                                                                                                      aCTTTTGTCAGAACCAGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa  >  2:200328/1‑140 (MQ=255)                                                                                                                          |                                                                                                                                        TGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  USA300TCH1516_ALE/638493‑638749

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