Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	763,202	C→T	I9I (ATC→ATT)	ECD_RS03830 →	metallophosphoesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	763,202	0	C	T	100.0%	115.2 / NA	48	I9V (ATC→GTT)	ECD_RS03830	metallophosphoesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/26);  total (22/26)

GAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762921‑763311                                                                                                                                                                                                                                                                                          |                                                                                                              gcgaATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCg                                                                                              <  2:544757/297‑1 (MQ=37) gcgaATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCg                                                                                              <  1:482091/297‑1 (MQ=37)     aTTAACCGTGTTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTTAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGTGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACgg                                                                                         <  2:85766/300‑1 (MQ=25)     aTTAACCGTGGTGAGGCGATACCTGAACCAATAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAACTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGCGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACgg                                                                                         <  2:85781/300‑1 (MQ=25)            gtgttgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTAt                                                                                    <  1:22035/294‑1 (MQ=25)            gtggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTAt                                                                                    <  1:22004/298‑1 (MQ=25)            gtggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTa                                                                                  >  1:100087/1‑300 (MQ=25)             tggttAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:498462/300‑1 (MQ=25)             tggtgAGGCGATACCTGNACCAGTAAAACAACTTCCNGTCATGTGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCCAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTNGCGAACGGTGTATTAc                                                                                 <  1:688920/300‑1 (MQ=25)             tggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:543158/300‑1 (MQ=25)             tggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGTGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:457599/300‑1 (MQ=25)             tggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:435108/300‑1 (MQ=25)             tggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCCCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAGCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGGGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:114713/300‑1 (MQ=25)               gtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCg                                                                               <  1:472840/300‑1 (MQ=25)                   ggCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGTACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGTCTAAAATGTCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCg                                                                               <  2:239119/296‑1 (MQ=25)                   ggCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGttc                                                                           >  1:458541/1‑299 (MQ=25)                         accTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTGTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAACCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTaac                                                                     >  2:419262/2‑298 (MQ=21)                                ccAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAgaa                                                               >  2:525190/1‑299 (MQ=21)                                   gTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCCGGCGAACGGTGTATTATCGGTTCGCTACCAGAGAAGAAc                                                           <  1:147227/300‑1 (MQ=21)                                       aacaacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAATGGTGTATTACCGGTTCGCTACCAGagaagagtggga                                                       >  1:284179/1‑294 (MQ=21)                                       aacaacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAATGGTGTATTACCGGTTCGCTACCAGAGAAGAAGGgga                                                       >  1:284686/1‑300 (MQ=21)                                        acaacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGggaa                                                      <  1:619734/300‑1 (MQ=21)                                                      tGGGCGGTAGACCTCTGAATCGTGCACAGGCTCTGGCGAAGATCGTAGAAATCAAGGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCAc                                        <  2:804241/300‑1 (MQ=18)                                                      tGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCAc                                        <  1:757118/300‑1 (MQ=18)                                                        ggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAGATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCa                                         >  2:730892/1‑297 (MQ=18)                                                         gCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGGAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGCGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGaa                                     >  2:550400/1‑300 (MQ=18)                                                                 ccTCTACATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATCAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAAACGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGATTCGCTAACAGAGAAGAACGGGAAGGAAAGGTGAGCAAGAATCtgactt                             >  2:96840/1‑297 (MQ=17)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:803012/248‑1 (MQ=21)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 >  2:803012/1‑248 (MQ=21)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:452157/248‑1 (MQ=21)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  2:173917/1‑298 (MQ=17)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  1:247917/1‑298 (MQ=17)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCATGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  2:99894/1‑298 (MQ=17)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTaa                                                                                 >  2:452157/1‑247 (MQ=21)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTAACAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  2:353505/1‑298 (MQ=17)                                                                 ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGATCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAACCtgattt                             >  2:97716/1‑298 (MQ=17)                                                                       aaTCGTGCACGGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCATAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaagg                       <  2:100087/300‑9 (MQ=17)                                                                        aTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaagga                      >  1:70750/1‑291 (MQ=17)                                                                            tGCACAGGCTCTGTCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggag                     <  1:451194/297‑11 (MQ=14)                                                                                  tgttctGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGTCAAAGAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCATAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccag            <  2:182556/297‑20 (MQ=11)                                                                                  ggCTCTGGCGAAGATCGCAGAAATCAGAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccag            <  2:182546/300‑20 (MQ=14)                                                                                   gCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccaga           >  2:731314/1‑280 (MQ=14)                                                                                    ctctGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgacaga           >  2:385370/1‑279 (MQ=14)                                                                                      ctGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagt         >  2:677579/1‑277 (MQ=14)                                                                                       ttgCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagtgc       <  2:479506/298‑25 (MQ=11)                                                                                        ggCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagtgcc      >  1:493382/1‑275 (MQ=11)                                                                                        ggCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTANTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagtgcc      >  1:594673/1‑275 (MQ=11)                                                                                            aaGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagtgccgcga  <  1:110562/300‑30 (MQ=11)                                                                                                                                                                                                                                                                                          |                                                                                                              GAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762921‑763311

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