Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	763,181	A→G	A2A (GCA→GCG)	ECD_RS03830 →	metallophosphoesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	763,181	0	A	G	100.0%	174.1 / NA	64	A2A (GCA→GCG)	ECD_RS03830	metallophosphoesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/33);  total (31/33)

CGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762894‑763311                                                                                                                                                                                                                                                                                                |                                                                                                                                   tgccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAATAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaat                                                                                                                         <  2:726736/281‑4 (MQ=255)  gccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcc                                                                                                                       <  1:106839/283‑6 (MQ=255)   ccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcct                                                                                                                      >  2:757118/17‑294 (MQ=255)      gtgaactgtcctgtatGACCGCGCGAATTAACCTTGGTTAGTCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGTTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACGGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggt                                                                                                                   <  2:321134/287‑10 (MQ=25)       tgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtc                                                                                                                  <  1:556104/288‑11 (MQ=37)        gaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGTCTAAAATGTCGCGGGCAggaatcctggtc                                                                                                                  >  1:239119/12‑289 (MQ=255)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCCTAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAAAAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  1:335595/7‑284 (MQ=37)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  1:115390/7‑284 (MQ=37)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  2:201671/7‑284 (MQ=37)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCAGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtagttgat                                                                                                            >  2:33979/7‑284 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggt                                                                                                         >  1:6758/4‑281 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGATAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCNGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggt                                                                                                         >  1:151375/4‑281 (MQ=37)                 tgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggta                                                                                                        <  2:400451/298‑21 (MQ=37)                 tgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGAAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggta                                                                                                        <  1:502828/298‑21 (MQ=37)                     tGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggtaaggt                                                                                                    <  1:790236/300‑25 (MQ=37)                            gcgaATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCg                                                                                              <  2:544757/297‑1 (MQ=37)                            gcgaATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCg                                                                                              <  1:482091/297‑1 (MQ=37)                             cgaATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCTTGTCATGGGCGGTAGAACTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGGCGTTGCCGCGGTAACAGGTGCAACCGTAACCAACATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcatggtcgttgatggtaaggtct                                                                                                  >  2:106577/2‑268 (MQ=25)                                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                                                                                 >  2:803012/1‑248 (MQ=21)                                                                                            ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:803012/248‑1 (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)                                                                                                                                                                                                                                                                                                |                                                                                                                                   CGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762894‑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: ‑