Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,462,758	G→A	A70T (GCA→ACA)	USA300HOU_RS12680 →	Putative 3‑methyladenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,462,758	0	G	A	100.0%	46.6 / NA	17	A70T (GCA→ACA)	USA300HOU_RS12680	Putative 3‑methyladenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/6);  total (11/6)

CAAACGTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCT  >  USA300TCH1516_ALE/2462629‑2462893                                                                                                                                  |                                                                                                                                        cAAACGTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGATGATACCATTCAAAAACATGTCAt                                                                                                                               >  2:140731/1‑140 (MQ=255)                                                   tGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTg                                                                             >  1:143779/1‑140 (MQ=255)                                                        gTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGc                                                                        >  2:4361/1‑140 (MQ=255)                                                            ggCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTTAAGGTATAACTAAAGGAGTAc                                                                    >  2:265144/1‑140 (MQ=255)                                                             gCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATATGAAGGTATACCTGAAGGCGTaat                                                                   >  1:248188/1‑138 (MQ=255)                                                                               aaaaTAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAAcc                                                 <  1:229141/140‑1 (MQ=255)                                                                               aaaaTAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAAcc                                                 <  1:247906/140‑1 (MQ=255)                                                                                 aaTAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCATACTTATCCGAGCAATTGAACCag                                               >  2:53406/1‑140 (MQ=255)                                                                                  aTAACACCTAAAAACACGTCATTATATAAACGTGGTAGTACAAATTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTAAAGGTATACCTAAAGGCATACTTATCCGCGCAATTGAACCaga                                              >  2:25632/1‑140 (MQ=255)                                                                                        ccTAAAGTCACGTCATTATATAAATGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAgg                                        <  1:140731/140‑1 (MQ=255)                                                                                                     cATTATATAAACGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATg                           <  1:186127/140‑1 (MQ=255)                                                                                                      aTTATATAAACATAGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTGTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTTAACCAGAAGAAGGTTTATCCGCTATGt                          <  1:35631/140‑1 (MQ=255)                                                                                                             aaaCGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCGATGTTCCGTaa                   >  1:101582/1‑140 (MQ=255)                                                                                                              aaCGTGGTGGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAAc                  >  1:247797/1‑140 (MQ=255)                                                                                                                 gtggtgGTACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACaga               >  1:223975/1‑140 (MQ=255)                                                                                                                        tACAATTTATACACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTCTCCGCTATTCTCCGTAACagatctaata        >  1:80495/1‑133 (MQ=255)                                                                                                                              ttATACACATGTCATGCATACGCATTTACTCATTAATTGTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCt  <  1:25859/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        CAAACGTATACAGGCTACATCGTGGAAACGGAAGCTTACTTAGGTTTGAATGATCGTGCGGCTCATGGCTATGGCGGTAAAATAACACCTAAAGTCACGTCATTATATAAACGTGGTGGTACAATTTATGCACATGTCATGCATACGCATTTACTCATTAATTTTGTAACAAAATCTGAAGGTATACCTGAAGGCGTACTTATCCGCGCAATTGAACCAGAAGAAGGTTTATCCGCTATGTTCCGTAACAGAGGTAAGAAAGGCT  >  USA300TCH1516_ALE/2462629‑2462893

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