Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_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_pae_tpiA	275,065	0	G	A	100.0%	39.5 / NA	17	T18I (ACC→ATC)	insH1	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/4);  total (13/4)

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCC  >  NC_000913_3_pae_tpiA/275005‑275202                                                             |                                                                                                                                          gACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAg                                                   <  1:278531/149‑1 (MQ=17)  aCAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAATTGGTGACTCATGATGACCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACTTTCCGGAgg                                                  >  2:376891/1‑149 (MQ=11)    aaGAAAATCTCTTTTCTGGTCTGCCGGCGCCTACTGCTGAATTCACTCTCGGCGAGGATGAGTTGCTGACTCATGCTGAACCCTGTCCTATGGCTCCAGGGTACAAACCTGATCTCATATCCGGGACTTGTCCGCCCCTTCCGGAGGCg                                                >  2:208861/1‑149 (MQ=11)    aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCg                                                >  1:550355/1‑149 (MQ=25)             tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATACGTTGATGACTCATGATGAACCGTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCGCAGGCGTTATGAGCt                                       >  2:331995/1‑149 (MQ=14)                      gTCTGACGGCGCTTACTGCTGAATGCACTGTCGGCGAAGATAAGTTGATGACTCAAGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGCGCTGGCGGCCCt                              >  2:132221/1‑149 (MQ=255)                         tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTGGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCtttt                           >  1:307159/1‑149 (MQ=255)                         tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAAACTTTATCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCtttt                           >  1:308211/1‑149 (MQ=17)                          gACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                          <  1:481300/149‑1 (MQ=255)                                  ttACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                          <  1:134856/141‑1 (MQ=255)                                  ttACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                          >  2:134856/1‑141 (MQ=255)                                           aaTTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATc         <  1:647469/149‑1 (MQ=255)                                            aTTCACTGTCGGCGAAGATAAGTTGAAGACTCATGATGAACCCTGTTCTATGGCACCAGATGACAAACATGACCTCATCTCAGGGACTTGCTCGCACCATCCGGAGGCGTTATGAGCTGGGGGCCCTTTTTGTTTCTGATTATTAATcc        >  1:743077/1‑149 (MQ=18)                                               cACTGTCGGCGAAGATAAGTTGATGACTCCTGTTGCACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTCAGCCCCa     >  2:365132/1‑149 (MQ=18)                                                aCTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATGAATCCCCAc    >  1:350421/1‑149 (MQ=255)                                                 cTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTTTATCTGATTATTAATCCCCAcc   >  1:72395/1‑149 (MQ=255)                                                  tGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAACTGCCTGCCCTTTTTGTATCTGATTATTAATCACCAccc  >  2:749702/1‑149 (MQ=18)                                                             |                                                                                                                                          GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCC  >  NC_000913_3_pae_tpiA/275005‑275202

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑