Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,609,130	A→C	K285Q (AAG→CAG)	insF‑3 →	putative transposase for IS3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,609,130	0	A	C	100.0%	194.2 / NA	47	K285Q (AAG→CAG)	insF‑3	putative transposase for IS3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/21);  total (26/21)

GCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTA  >  REL606/1609104‑1609202                           |                                                                         gcCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCa                         <  2:4791757/76‑1 (MQ=17) gcCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCa                         <  1:800425/76‑1 (MQ=17)   cTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCAcc                       >  2:3420690/1‑76 (MQ=25)   cTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCAcc                       >  2:2902477/1‑76 (MQ=25)   cTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCAcc                       >  2:1392828/1‑76 (MQ=25)    tCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCt                      >  1:1311420/1‑76 (MQ=255)    tCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCt                      <  1:4375739/76‑1 (MQ=255)     cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     <  1:183965/76‑1 (MQ=255)     cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     <  2:3523161/76‑1 (MQ=255)     cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     <  1:1432658/76‑1 (MQ=255)     cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     >  2:2093543/1‑76 (MQ=255)      aGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaa                    >  2:1390524/1‑76 (MQ=255)      aGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaa                    >  1:4877134/1‑76 (MQ=255)      aGTCCGGAACAATATGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaa                    >  2:896584/1‑76 (MQ=255)       gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   >  1:2660266/1‑76 (MQ=255)       gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   <  1:2138526/76‑1 (MQ=255)       gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   >  2:3719069/1‑76 (MQ=255)       gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   >  1:4928458/1‑76 (MQ=255)         ccGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATa                 <  1:806855/76‑1 (MQ=255)         ccGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATa                 <  1:1074996/76‑1 (MQ=255)          cGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATAc                >  1:2866188/1‑76 (MQ=255)            gAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGa              >  2:3228154/1‑76 (MQ=255)            gAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGa              <  2:4961039/76‑1 (MQ=255)             aaCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAg             <  1:389909/76‑1 (MQ=255)             aaCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAg             <  2:2428078/76‑1 (MQ=255)             aaCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAg             >  1:4303021/1‑76 (MQ=255)              aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            <  1:4881451/76‑1 (MQ=255)              aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            <  1:3492342/76‑1 (MQ=255)              aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            >  2:4316468/1‑76 (MQ=255)              aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            <  1:1658793/76‑1 (MQ=255)               cAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGc           <  1:2782381/76‑1 (MQ=255)               cAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGc           >  2:2635261/1‑76 (MQ=255)                 aTTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCAt         <  1:1483732/76‑1 (MQ=255)                  tttGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCAtt        <  2:3044173/76‑1 (MQ=255)                  tttGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCAtt        >  1:4628460/1‑76 (MQ=255)                   ttGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTg       >  2:770591/1‑76 (MQ=255)                    tGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGc      >  2:2394065/1‑76 (MQ=255)                    tGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGc      >  1:143872/1‑76 (MQ=255)                     gAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCt     >  1:1797303/1‑76 (MQ=255)                      aaaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTg    >  2:1293314/1‑76 (MQ=255)                      aaaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTg    >  1:3300429/1‑76 (MQ=255)                       aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   <  1:3061479/76‑1 (MQ=255)                       aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   >  2:4083010/1‑76 (MQ=255)                       aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   <  1:1658227/76‑1 (MQ=255)                       aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   >  1:5016428/1‑76 (MQ=255)                        aaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTa  <  2:3275178/76‑1 (MQ=255)                        aaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTa  <  2:3396815/76‑1 (MQ=255)                        aaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTa  >  2:2701067/1‑76 (MQ=255)                           |                                                                         GCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTA  >  REL606/1609104‑1609202

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