Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	763,166	C→T	intergenic (+104/‑10)	ECD_RS03820 → / → ECD_RS03830	replication protein P/metallophosphoesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	763,166	0	C	T	100.0%	177.1 / NA	66	intergenic (+104/‑10)	ECD_RS03820/ECD_RS03830	replication protein P/metallophosphoesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/33);  total (33/33)

CCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762885‑763311                                                                                                                                                                                                                                                                                          |                                                                                                                                                  gcgcaaggctgccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGcgc                                                                                                                                   >  1:615550/28‑296 (MQ=255) gcgcaaggctgccgatgaactgtcctgtatGACCGCGCGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGAAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGAAGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCATAGGACGTTGCCGCGGTAACAGGTGCAAACGTAAACAGCCTAAATCAGGCTGCGGCTAAAATGGc                                                                                                                                     >  2:453105/28‑296 (MQ=255)          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                             >  1:247917/1‑298 (MQ=17)                                                                                                     ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGATCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGTGCAACCGTAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAACGGGAAGGAAAGGTGAGCACGAATCtgattt                             >  2:173917/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)                                                                                                                                                                                                                                                                                          |                                                                                                                                                  CCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGATGAACAAACTGGATACGATTGGATTCGACA  >  NC_012971/762885‑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: ‑