Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,647,521	C→T	G86E (GGA→GAA)	unknown ←	Transposase, IS200 family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,647,521	0	C	T	100.0%	65.1 / NA	26	G86E (GGA→GAA)	unknown	Transposase, IS200 family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/12);  total (14/12)

AACAAACTATAATAATCTTTGTTTACTTACTGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTTGAAAGGGGTCCACATATCCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTCCACTGTCTTTGTATTTAGCCCTACCG  >  1/1647295‑1647608                                                                                                                                                                                                                                   |                                                                                        aacaaacTATAATAATCTTTGTTTACTTACTGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAAt                                                                <  1:631322/251‑1 (MQ=255)         ataataATCTTTGTTTACTTACTGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACCAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCtt                                                        >  1:519128/1‑251 (MQ=255)             taatCTTTGTTTACTTACTGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAAt                                                                          >  1:108958/1‑229 (MQ=255)                tCTTTGTTTACTTACTGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCtttt                                                  >  1:5896/1‑250 (MQ=255)                         acttactGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTc                                          >  1:127330/1‑249 (MQ=255)                              ctGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTtgtat                                   <  1:800039/251‑1 (MQ=255)                              ctGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTtatat                                   >  1:299360/1‑251 (MQ=255)                                gCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTtatata                                  >  1:315518/1‑250 (MQ=255)                                 ccGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATAc                                 <  1:444977/250‑1 (MQ=255)                                  cGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATcc                                                                     >  1:142111/1‑213 (MQ=255)                                    aaCTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATActcc                              <  1:308845/250‑1 (MQ=255)                                    aaCTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATActc                               >  1:830410/1‑249 (MQ=255)                                        aaCAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTacac                         <  1:178769/251‑1 (MQ=255)                                        aaCAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTacac                         <  1:74541/251‑1 (MQ=255)                                          cAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAAt                                                                <  1:171817/210‑1 (MQ=255)                                          cAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTacac                         <  1:654396/249‑1 (MQ=255)                                           agttagttCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTc                                          >  1:359192/1‑231 (MQ=255)                                                 ttCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTACACtacctttg                 <  1:555476/250‑8 (MQ=255)                                                     ggTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTACACTACCTTTGTAt              >  1:91189/1‑249 (MQ=255)                                                         aaaTTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTcc                                                      >  1:525547/1‑205 (MQ=255)                                                         aaaTTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATAct                                >  1:224001/1‑227 (MQ=255)                                                         aaaTTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTACACTACCTTTGTATTTAGcc        <  1:88441/251‑1 (MQ=255)                                                           aTTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATAc                                 >  1:527362/1‑224 (MQ=255)                                                           aTTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTACACTACCTTTGTATTTAGcc        >  1:173696/1‑249 (MQ=255)                                                             tACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTACACTACCTTTGTATTTAGCCCTAc    <  1:239399/251‑1 (MQ=255)                                                               cacCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTT‑AAAGGGGTCCACATATTCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTACACTACCTTTGTATTTAGCCCTACCg  <  1:901509/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                        AACAAACTATAATAATCTTTGTTTACTTACTGCCGAACTAACAGTTAGTTCAGGTAAAATTACACCCTGTACAAAAAAATCATAAAATAACAATCCGCACGTAAAACGTGCAGTTTTTAGGCTACTCTATAAAGACCTAGTACTTCACATGCCCTTCAAAGGGCATATGAAAACCAACAACCGTGAAATAACTCTATGACCTACCCCTTGAAAGGGGTCCACATATCCTTTCTTCGATAAATTATCCTGAATCATGTCTTCCTTTTCCTGATTCCTTATATACTCCTCCACTGTCTTTGTATTTAGCCCTACCG  >  1/1647295‑1647608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑