Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,513,930	A→T	100%	L413F (TTA→TTT)	USA300HOU_RS07375 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,513,930	0	A	T	100.0%	57.7 / NA	19	L413F (TTA→TTT)	USA300HOU_RS07375	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (6/13);  total (6/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAA  >  USA300TCH1516_ALE/1513803‑1514059                                                                                                                                |                                                                                                                                  tatgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                                                                      >  1:190512/1‑141 (MQ=255) tatgtatAGATATATTAGTAGAAAAAGGAATGAGGTTTGTATCTTAAATAAAGCGGACAAAAAATACTGCAATAAAAGTAACAAATTGAGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                                                                      <  2:366693/141‑1 (MQ=255)  atgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatt                                                                                                                     <  2:300365/141‑1 (MQ=255)       tCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTaa                                                                                                                >  1:72431/1‑141 (MQ=255)          cTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGc                                                                                                             >  1:105318/1‑141 (MQ=255)             tattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATAACTCATAATATGTTTTATTGTATTAGGATTATTAATGTCat                                                                                                          >  2:328590/1‑141 (MQ=255)              attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatt                                                                                                         >  2:310480/1‑141 (MQ=255)                             cATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                                                          <  1:133329/141‑1 (MQ=255)                             cATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                                                          <  1:170211/141‑1 (MQ=255)                                   gTTTGGCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCCCTTACAGGACACAGAATCTGCTGTCATGAGGTGGGATTATTAATGCCATTATTTGGATTATCACTTATCCt                                                                                    <  1:117846/141‑1 (MQ=255)                                                 tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCTCTCATTATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                                                                      <  1:377533/141‑1 (MQ=255)                                                                       aaaaaaaATACCAAATCGGGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGa                                                <  2:220344/138‑1 (MQ=255)                                                                             aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATAAAAAAATAACAAGGTATTTATAATTGAATTAATATTATATATTAAAGAtcgtcg                                          >  1:469937/1‑141 (MQ=255)                                                                                aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGc                                       <  1:89737/141‑1 (MQ=255)                                                                                      aagtgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACa                                 <  2:295219/137‑1 (MQ=255)                                                                                                   ccAGCATCTATCTCACTCATAATATGTTTTATTGTAGAAGGAGTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                    <  1:469086/141‑1 (MQ=255)                                                                                                   ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                    <  1:396832/141‑1 (MQ=255)                                                                                                                 aCTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGAc      <  1:215794/141‑1 (MQ=255)                                                                                                                     ataataTGTTTTATTGTATTAGGATTATTAATGCTATTAATTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGaaa  <  2:389912/141‑1 (MQ=255)                                                                                                                                |                                                                                                                                  TATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAA  >  USA300TCH1516_ALE/1513803‑1514059

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