Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	275,065	G→A	T18I (ACC→ATC)	insH1 ←	IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	275,065	0	G	A	100.0%	40.4 / NA	16	T18I (ACC→ATC)	insH1	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/7);  total (9/7)

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAC  >  NC_000913_3_vch_tpiA/275005‑275200                                                             |                                                                                                                                        gataaGAAAATCTCTTTTCTGGTCTGACGGCGCCTACTGCTGAATTCACTGTCCGCGAAGATAAGTTGATGACTCATGATCAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGTACCTTCCGGAg                                                 <  2:537490/146‑1 (MQ=11)    aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCGGAGGCg                                              >  2:752054/1‑149 (MQ=25)    aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCg                                              >  1:151889/1‑149 (MQ=25)    aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCg                                              >  1:803827/1‑149 (MQ=25)     aGAAAATCTCTTTTCTGGTCTGTCTGCGCTTACTGCTGAATTCACTATCGGCGAAGATAAATTGATGACTCATGATGAACCCTGTTCTATAGCACCAGATAACAAACATTATCTCATATCCGGGCCTTGTTCGCACCTTCCCGAGGCGt                                             >  2:234489/1‑149 (MQ=11)       aagaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTa                                           <  1:236144/146‑1 (MQ=255)          aTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGa                                        <  1:62400/149‑1 (MQ=255)                       tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATCAGCTGGCGGCCCtt                           >  2:139970/1‑149 (MQ=255)                         tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCtttt                         <  1:616207/149‑1 (MQ=255)                          aacgacgCTTACTTCTGAAGTCACTGTCCGCGAAGATAAGTTGATGACTCAGGATGAACCCTGTTCTATGTCTCCATATCACAAACATCATCTCATATCAGCAACTTTTTCCCACCTTCCGTATGCGTTATGAGCTGGCGGCCCttttt                        <  2:304079/148‑1 (MQ=255)                             ggCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTa                     <  1:649512/149‑1 (MQ=255)                                     ccgctgAATTCACTGGGGGGGAAGATAAGTTGATGACTCATGATGAACTCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatta             <  2:169453/147‑1 (MQ=18)                                        ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGAttatta          >  1:369329/1‑149 (MQ=255)                                           aaTTCACTGTCGGCAAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGCCGTTATGAGCTGGCGGCCCTTTTTGTATCGGATTATTAATc       >  2:515739/1‑149 (MQ=255)                                             ttCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATccc     >  2:117920/1‑149 (MQ=255)                                                aCTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTTTGAGCTGGCGGCCCTTTTTGTTCTTGATTATTATTccccct  >  1:10360/1‑147 (MQ=17)                                                             |                                                                                                                                        GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAC  >  NC_000913_3_vch_tpiA/275005‑275200

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