Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	514,909	C→T	intergenic (+570/‑168)	unknown → / → unknown	Integrase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	514,909	0	C	T	100.0%	192.3 / NA	69	intergenic (+570/‑168)	unknown/unknown	Integrase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/32);  total (37/32)

GCCGACTATCAGCTTCGGTTGCAGCGTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAG  >  1/514665‑515089                                                                                                                                                                                                                                                     |                                                                                                                                                                                     gCCGACTATCAGCTTCGGTTGCAGCGTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCaa                                                                                                                                                                                >  1:698424/1‑251 (MQ=255)                  gTTGCAGCGTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAAc                                                                                                                                                                              <  1:72046/236‑1 (MQ=255)                          gTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGt                                                                                                                                                       >  1:442439/1‑251 (MQ=255)                                aCCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGt                                                                                                                                                 >  1:490786/1‑251 (MQ=255)                                 cccGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                >  1:656487/1‑251 (MQ=255)                                 cccGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                >  1:689712/1‑251 (MQ=255)                                 cccGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                >  1:845976/1‑251 (MQ=255)                                   cGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGa                                                                                                                                                                    <  1:639863/229‑1 (MQ=255)                                   cGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAg                                                                                                                                              >  1:199425/1‑251 (MQ=255)                                    gATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                             >  1:783682/1‑251 (MQ=255)                                     aTAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGa                                                                                                                                                                    <  1:304500/227‑1 (MQ=255)                                     aTAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATc                                                                                                                                                    <  1:711753/243‑1 (MQ=255)                                                 gCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAg                                                                                                                                                                            <  1:333025/207‑1 (MQ=255)                                                    tGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAgg                                                                                                                                                                           >  1:41959/1‑205 (MQ=255)                                                       ttACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGa                                                                                                                           <  1:649568/250‑1 (MQ=255)                                                          cATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTa                                                                                                                                                      >  1:682684/1‑220 (MQ=255)                                                                  aTTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTtct                                                                                                                >  1:79872/1‑250 (MQ=255)                                                                            aTTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTa                                                                                                     <  1:253078/251‑1 (MQ=255)                                                                                    aTTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTgg                                                                                                                                                        >  1:49601/1‑192 (MQ=255)                                                                                    aTTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                             <  1:723290/203‑1 (MQ=255)                                                                                         gttgATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAAc                                                                                                                                                                              <  1:457569/165‑1 (MQ=255)                                                                                                 cAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGaa                                                                                 <  1:619853/250‑1 (MQ=255)                                                                                                  aaCGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                                               <  1:54247/251‑1 (MQ=255)                                                                                                    cGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAagag                                                                             <  1:342170/251‑1 (MQ=255)                                                                                                      ggCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAggg                                                                           >  1:244119/1‑251 (MQ=255)                                                                                                           cAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGa                                                                                                                                                                    <  1:499887/157‑1 (MQ=255)                                                                                                             ggTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCg                                                                                                                    <  1:598584/203‑1 (MQ=255)                                                                                                               tAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                   <  1:359912/250‑1 (MQ=255)                                                                                                                 ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAg                                                                                                            >  1:177415/1‑207 (MQ=255)                                                                                                                 ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGt                                                                >  1:40265/1‑251 (MQ=255)                                                                                                                 ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGt                                                                >  1:349809/1‑251 (MQ=255)                                                                                                                 ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGt                                                                >  1:34491/1‑251 (MQ=255)                                                                                                                     aTTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTaaaa                                                            <  1:381233/251‑1 (MQ=255)                                                                                                                             aCTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGcgaa                                                                                                                                      <  1:675338/169‑4 (MQ=255)                                                                                                                                   aaTTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATaa                                               <  1:369192/250‑1 (MQ=255)                                                                                                                                     ttAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc                                            <  1:529715/251‑1 (MQ=255)                                                                                                                                      tAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                                           <  1:8185/251‑1 (MQ=255)                                                                                                                                                     aTTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt                            <  1:330172/251‑1 (MQ=255)                                                                                                                                                     aTTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt                            <  1:764886/251‑1 (MQ=255)                                                                                                                                                       tttCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                          <  1:367703/251‑1 (MQ=255)                                                                                                                                                         tCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg                        <  1:126220/251‑1 (MQ=255)                                                                                                                                                         tCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg                        <  1:173318/251‑1 (MQ=255)                                                                                                                                                          ccTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAAt                                                                     <  1:809433/205‑1 (MQ=255)                                                                                                                                                          ccTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc                                            >  1:65219/1‑230 (MQ=255)                                                                                                                                                          ccTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAc                                             >  1:700199/1‑229 (MQ=255)                                                                                                                                                            ttATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAg                     <  1:735380/251‑1 (MQ=255)                                                                                                                                                             tataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                                           <  1:148805/228‑1 (MQ=255)                                                                                                                                                             tataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc                                            >  1:723118/1‑227 (MQ=255)                                                                                                                                                             tataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAg                     >  1:563030/1‑250 (MQ=255)                                                                                                                                                              ataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTgg                                                                                                                                                        <  1:824448/118‑1 (MQ=255)                                                                                                                                                                   cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                             >  1:755469/1‑124 (MQ=255)                                                                                                                                                                   cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGACAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                                           >  1:453298/1‑222 (MQ=255)                                                                                                                                                                   cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                                                 >  1:169288/1‑216 (MQ=255)                                                                                                                                                                   cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                                                 >  1:243785/1‑216 (MQ=255)                                                                                                                                                                   cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTACGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGa              >  1:89271/1‑251 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCAtt                                                                                                       >  1:210003/1‑160 (MQ=11)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGcc                                                                                          <  1:506908/173‑1 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCt                                                                                         >  1:528298/1‑174 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                                               >  1:90337/1‑184 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                   >  1:329445/1‑196 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                   <  1:768899/196‑1 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAAt                                                           >  1:635466/1‑204 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                                                 >  1:413223/1‑214 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAc                                             <  1:586638/218‑1 (MQ=255)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:292542/1‑251 (MQ=16)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:249194/1‑251 (MQ=16)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:645188/1‑251 (MQ=16)                                                                                                                                                                     ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:551074/1‑251 (MQ=16)                                                                                                                                                                               gCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaaag  >  1:883966/1‑251 (MQ=9)                                                                                                                                                                                                                                                     |                                                                                                                                                                                     GCCGACTATCAGCTTCGGTTGCAGCGTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAG  >  1/514665‑515089

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