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%	126.8 / NA	46	A2A (GCA→GCG)	ECD_RS03830	metallophosphoesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/18);  total (28/18)

CGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762894‑763311                                                                                                                                                                                                                                                                                                |                                                                                                                                   tgccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaat                                                                                                                         <  2:169414/281‑4 (MQ=255)   ccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAAACAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcct                                                                                                                      >  2:54710/17‑294 (MQ=255)        gaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtc                                                                                                                  >  2:535533/12‑289 (MQ=37)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGCCCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  2:7666/7‑284 (MQ=37)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  1:13624/7‑284 (MQ=37)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAAACAGCATAAATCAGGCTGCGGCTAAAATGACGCGGGCAggaatcctggtcgttgat                                                                                                            >  1:302732/7‑284 (MQ=255)             gtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGAAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAAAAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  2:7518/7‑284 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCCGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggt                                                                                                         >  2:342684/4‑281 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGTGCGGGCAggaatcctggtcgttgatggt                                                                                                         >  2:402585/4‑281 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggt                                                                                                         >  1:19504/4‑281 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGAGCGGGCAggaatcctggtcattgatggt                                                                                                         >  1:230673/4‑281 (MQ=37)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGAAGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAAACGTAAACAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtagttgatggt                                                                                                         >  2:342700/4‑281 (MQ=255)                ctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTAGTGGGCGGTAGAATTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATAGCGCGGGCAggaatcctggtcgttgat                                                                                                            >  1:231451/4‑281 (MQ=37)                 tgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGTCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggt                                                                                                         <  2:687667/297‑20 (MQ=37)                     tGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAggaatcctggtcgttgatggtaaggt                                                                                                    <  1:504824/300‑25 (MQ=37)                               aaTTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGCTCGAAGAAATCAAAGCTAAGTTCGGACTGAAAGGATTAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTCACAGGTGCAACCGTAACCCGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGAAGGAATCCTGGTCATTTATGGTAAGGTCTGGCGAAc                                                                                           >  2:204451/1‑299 (MQ=25)                                       gtggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTa                                                                                  >  1:230382/1‑300 (MQ=25)                                        ttgtgatgtgataccTGAACCAGTAAAACAGCTTCCTGTCATGGGCGGGAGACCTCTAAATCATGCACAGGCTCTGTCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTTGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGTTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:46906/287‑1 (MQ=21)                                        tggtgATGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:375347/300‑1 (MQ=25)                                        tggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAATAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGAGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:89392/300‑1 (MQ=25)                                        tggtgAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  1:3798/300‑1 (MQ=25)                                        tggggAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTTTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:9164/296‑1 (MQ=25)                                         ggtgAGGCGATACCTGAACCAGTAAAACAACGTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGTCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAc                                                                                 <  2:754731/299‑1 (MQ=25)                                                    accTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAAACGTAACCAGCATAAATCAAGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGttcgttac                                                                      >  2:283937/2‑293 (MQ=25)                                                    accTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTACAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATCACCGGTTCGCtccc                                                                     >  2:283172/2‑297 (MQ=21)                                                              gTAAAACAACTTCCTGTCATGTGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAATAGGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAAc                                                           <  2:416888/300‑1 (MQ=21)                                                                             gTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGa                                            >  1:178461/1‑300 (MQ=18)                                                                              tCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAg                                           <  1:199502/300‑1 (MQ=18)                                                                                 tGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCAc                                        <  2:57940/300‑1 (MQ=18)                                                                                    gCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGaa                                     >  1:201004/1‑300 (MQ=18)                                                                                            ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTNCCTGGGGACGCATAAGAGCTTCTGTGCACAGGNCGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  2:601846/1‑298 (MQ=17)                                                                                            ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGCACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGaatctaattt                             >  2:299776/1‑292 (MQ=17)                                                                                            ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGaagcacaaccttgctt                             >  1:265942/1‑285 (MQ=18)                                                                                            ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGTACGaatttgattt                             >  2:309829/1‑292 (MQ=17)                                                                                            ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  2:699355/1‑298 (MQ=17)                                                                                              tctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattttt                           >  1:751664/1‑296 (MQ=17)                                                                                               ctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattttta                          >  1:458442/1‑295 (MQ=17)                                                                                                tAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAGCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaa                         <  2:476892/300‑7 (MQ=17)                                                                                                   aTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaagga                      >  1:136731/1‑291 (MQ=17)                                                                                                    tCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggag                     <  2:668458/300‑11 (MQ=17)                                                                                                            aGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGTACTGAAAGGATCAAGTGTATGACGGTCAAATATGCAATTATTCATTACCTGTGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgcc              <  2:426295/299‑18 (MQ=11)                                                                                                            aGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgcc              <  1:472531/299‑18 (MQ=14)                                                                                                             ggCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgcca             >  1:324644/1‑281 (MQ=14)                                                                                                                tctGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagt         <  2:555097/300‑23 (MQ=14)                                                                                                                 ctGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTATGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatctttaaggagtgtcgccagagt         >  2:500263/1‑277 (MQ=14)                                                                                                                       aaGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgatttttaaggagtgtcgccagagtgccgcga  <  1:473659/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 < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑