Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	514,955	T→G	intergenic (+616/‑122)	unknown → / → unknown	Integrase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	514,955	0	T	G	100.0%	128.6 / NA	47	intergenic (+616/‑122)	unknown/unknown	Integrase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/22);  total (25/22)

TTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAG  >  1/514719‑515089                                                                                                                                                                                                                                             |                                                                                                                                       ttACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGa                                                                                                                           <  1:649568/250‑1 (MQ=255)            aTTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTtct                                                                                                                >  1:79872/1‑250 (MQ=255)                      aTTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTa                                                                                                     <  1:253078/251‑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)                                                       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)                                                                             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)                                                                                                             cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGACAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                                           >  1:453298/1‑222 (MQ=255)                                                                                                             cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                                                 >  1:243785/1‑216 (MQ=255)                                                                                                             cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                                                 >  1:169288/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:768899/196‑1 (MQ=255)                                                                                                               ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                   >  1:329445/1‑196 (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:645188/1‑251 (MQ=16)                                                                                                               ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:551074/1‑251 (MQ=16)                                                                                                               ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:292542/1‑251 (MQ=16)                                                                                                               ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg            >  1:249194/1‑251 (MQ=16)                                                                                                                         gCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaaag  >  1:883966/1‑251 (MQ=9)                                                                                                                                                                                                                                             |                                                                                                                                       TTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAG  >  1/514719‑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: ‑