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%	43.4 / NA	15	T18I (ACC→ATC)	insH1	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/6);  total (9/6)

AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAT  >  NC_000913_3_vch_tpiA/275008‑275207                                                          |                                                                                                                                                aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCg                                                     >  1:455339/1‑149 (MQ=25)   gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt                                                   <  2:334879/149‑1 (MQ=37)       aTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTtgtgc                                               >  2:294843/1‑148 (MQ=255)            ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGCGGCGTTATGAGCTgg                                          >  2:24973/1‑149 (MQ=255)               tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCATGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATAAGGGACTTGTTAGCCCCTTCCGGAGGCGTTATGAGCTggcgg                                       >  1:643437/1‑149 (MQ=14)                     cTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttt                                 <  1:346628/149‑1 (MQ=255)                               ttACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGa                       <  1:669705/149‑1 (MQ=255)                               ttACTGCAGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACGTGATCTCATATCAGGGATTTGTTCGCACCTTCAGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATTTGa                       >  1:622952/1‑149 (MQ=17)                                  ctgctgAATTCATTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTATATGGC‑CCACCTAGATCTTCACGATCTCTTATCATGGACTTGTTCGCTCCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatta                    <  2:412624/149‑1 (MQ=18)                                     ctgAATTCACTGTCGGCGAAGATAAGTTGATAACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTGGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGAttatta                 >  1:233064/1‑149 (MQ=255)                                        aaTTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATc              >  2:33212/1‑149 (MQ=255)                                          ttCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATccc            >  1:364762/1‑149 (MQ=255)                                            cACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCa          <  1:454502/149‑1 (MQ=255)                                              cTGTCGGGTAATATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAcc        <  2:492987/149‑1 (MQ=255)                                                    gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAt  >  2:566970/1‑149 (MQ=255)                                                          |                                                                                                                                                AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAT‑GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAT  >  NC_000913_3_vch_tpiA/275008‑275207

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