Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	265,333	G→A	P499P (CCG→CCA)	ugpQ_2 →	Glycerophosphodiester phosphodiesterase, cytoplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	265,333	0	G	A	100.0%	67.8 / NA	23	P499P (CCG→CCA)	ugpQ_2	Glycerophosphodiester phosphodiesterase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/11);  total (12/11)

GATTATACACAACGTGTTATTGATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCGTTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAATAAAGAAGTCTATACTTGGACT  >  USA300TCH1516_ALE/265196‑265453                                                                                                                                          |                                                                                                                         gATCATAAACAACGTGTTATTGATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTg                                                                                                                       <  2:226360/141‑1 (MQ=255)       acacAACGTGTTATTGATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGACAATGGTATCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGttt                                                                                                                 >  2:342915/1‑141 (MQ=255)                    ttGATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAg                                                                                                    <  2:96885/141‑1 (MQ=255)                    ttGATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAg                                                                                                    >  1:223926/1‑141 (MQ=255)                      gATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGaa                                                                                                  >  1:194498/1‑141 (MQ=255)                             tGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTCGTTATATCATTCCATTGCAGTTTGGTGATTTTAAAGAAAcatcat                                                                                           >  1:21567/1‑141 (MQ=255)                                                         tCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTc                                                               <  2:156065/141‑1 (MQ=255)                                                                        ggATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATCTTTCTTATTCGCCAAGACt                                                >  2:219138/1‑141 (MQ=255)                                                                              tgatgtgatgACTAAGTTGAAAAAAGAAGGGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTaa                                          <  2:63821/141‑1 (MQ=255)                                                                                    gatgaCTAAGTATAAAAAAGAAGCGTCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAgc                                    <  1:294135/141‑1 (MQ=255)                                                                                           aaGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTg                             >  1:353423/1‑141 (MQ=255)                                                                                              ttGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGaa                          >  1:313545/1‑141 (MQ=255)                                                                                                gAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGaaaa                        >  1:96271/1‑141 (MQ=255)                                                                                                 aaaaaaGAAGCGCCATATCTCAAGTGTGGTTATATTATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAAt                       <  2:246308/141‑1 (MQ=255)                                                                                                 aaaaaaGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAAt                       <  2:7873/141‑1 (MQ=255)                                                                                                 aaaaaaGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAAt                       <  2:290049/141‑1 (MQ=255)                                                                                                   aaaaGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAATaa                     <  2:215936/141‑1 (MQ=255)                                                                                                      agaagCGCCATATCTCAAGTGTGGTTATCTCATTCCATTGCAGTTTGGTCATTTTAAAGACACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCCTGACTTGTTAATCAAGCGCACTTGCAAAATAcaga                  >  2:153900/1‑141 (MQ=255)                                                                                                           ggcCATAGCGCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGTAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAATAAAGAAGTCt             <  2:57953/140‑1 (MQ=255)                                                                                                               atatCACAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACAGCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAATAAAGAAGTCTATAc         <  1:185658/141‑1 (MQ=255)                                                                                                                tatCTCAAGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAACCGCACTTGGAAGAGAGAGAACACTATACt        >  2:55200/1‑141 (MQ=255)                                                                                                                      aaGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAATAAAGAAGTCTATACTTGGACt  >  2:240571/1‑141 (MQ=255)                                                                                                                      aaGTGTGGTTATATCATTCCATTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCCCACTTGCAACATACAGAAGTCTATACTTGGACt  >  2:240570/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                         GATTATACACAACGTGTTATTGATATTTTGAAAAAGCATGGTGTTGAACATCAATATCGTGTGATGTCTTTGGATTATGATGTGATGACTAAGTTGAAAAAAGAAGCGCCATATCTCAAGTGTGGTTATATCATTCCGTTGCAGTTTGGTCATTTTAAAGAAACATCATTAGATTTCTTTGTCATCGAAGATTTTTCTTATTCGCCAAGACTTGTTAATCAAGCGCACTTGGAAAATAAAGAAGTCTATACTTGGACT  >  USA300TCH1516_ALE/265196‑265453

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