Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,633,288	C→T	G32S (GGC→AGC)	unknown ←	Mobile element protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,633,288	0	C	T	100.0%	177.2 / NA	64	G32S (GGC→AGC)	unknown	Mobile element protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/36);  total (28/36)

CGGTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCCTCGTTTTCCATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAAGGCTTTTAGCCAATAAATGACGAGCGAATGGG  >  1/1633143‑1633531                                                                                                                                                  |                                                                                                                                                                                                                                                    cGGTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTGGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGttt                                                                                                                                            <  1:442241/251‑1 (MQ=9)   gTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTg                                                                                                                                           <  1:234831/250‑1 (MQ=14)            aTCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGAtt                                                                                                                                 >  1:10336/1‑251 (MQ=255)             tCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGAttt                                                                                                                                >  1:692983/1‑251 (MQ=255)                    aCGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGc                                                                                                                         <  1:863768/251‑1 (MQ=255)                    aCGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGa                                                                                                                                   <  1:285616/241‑1 (MQ=25)                     cGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGct                                                                                                                        <  1:553444/251‑1 (MQ=255)                            aaaTAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGtt                                                                                                                                             >  1:30635/1‑223 (MQ=14)                            aaaTAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCtttt                                                                                                                 >  1:899815/1‑251 (MQ=255)                                  agCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATc                                                                                                            <  1:382515/250‑1 (MQ=255)                                      tttCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTgg                                                                                                       >  1:463951/1‑251 (MQ=255)                                       ttcttcCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATc                                                                                                            >  1:464094/1‑245 (MQ=255)                                                 aCTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTAttt                                                                                            >  1:219501/1‑251 (MQ=255)                                                    cGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAg                                                                                         <  1:359153/251‑1 (MQ=255)                                                      aaaCTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAggg                                                                                       <  1:62761/251‑1 (MQ=255)                                                           cAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGa                                                                                                                                   <  1:573677/202‑1 (MQ=255)                                                           cAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTAtttt                                                                                           >  1:607384/1‑242 (MQ=255)                                                             gTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAg                                                                                         >  1:554589/1‑242 (MQ=255)                                                                aCATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAggg                                                                              <  1:481328/250‑1 (MQ=255)                                                                    cGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtat                                                                         >  1:663511/1‑251 (MQ=255)                                                                       ttACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:851427/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:222304/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:901454/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:366529/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:365770/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:771553/250‑1 (MQ=255)                                                                         aCCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAg                                                                    >  1:669858/1‑251 (MQ=255)                                                                          cccATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAgg                                                                                        <  1:37554/230‑1 (MQ=255)                                                                          cccACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCTCATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAgg                                                                   >  1:218177/1‑251 (MQ=255)                                                                          cccACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAgg                                                                   >  1:510499/1‑251 (MQ=255)                                                                             acttttGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTc                                                                                                                                        <  1:738673/177‑1 (MQ=19)                                                                             aTTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATCCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGa                                                                >  1:558898/1‑251 (MQ=255)                                                                              cttttGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAggg                                                                              <  1:578190/235‑1 (MQ=255)                                                                              cttttGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGAt                                                               >  1:465967/2‑251 (MQ=255)                                                                                    attgaGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTAc                                                          <  1:882187/250‑1 (MQ=255)                                                                                                  cTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCct                                             >  1:744166/1‑249 (MQ=255)                                                                                                    gTAAAATCACGATTTAGCTCATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATc                                                                                                            <  1:460359/184‑1 (MQ=255)                                                                                                       aaaTCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCAc                                      <  1:725765/251‑1 (MQ=255)                                                                                                             cGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAGTCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAAgg                                <  1:48009/251‑1 (MQ=255)                                                                                                             cGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAggg                                                                              <  1:262570/205‑1 (MQ=255)                                                                                                               aTTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAggg                                                                  <  1:246880/215‑1 (MQ=255)                                                                                                               aTTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGaaatg                                >  1:638162/1‑247 (MQ=255)                                                                                                               aTTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGaaatg                                >  1:526596/1‑247 (MQ=255)                                                                                                               aTTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGaaatg                                >  1:737050/1‑247 (MQ=255)                                                                                                               aTTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAAgg                                >  1:210197/1‑249 (MQ=255)                                                                                                                  tAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGc                                                                                                      <  1:714138/176‑1 (MQ=255)                                                                                                                   aGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAggg                                                                  >  1:90513/1‑211 (MQ=255)                                                                                                                   aGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTa                          >  1:775324/1‑251 (MQ=255)                                                                                                                     cacaTTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGc                        <  1:492684/251‑1 (MQ=255)                                                                                                                      acaTTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGcc                       <  1:155350/251‑1 (MQ=255)                                                                                                                       caTTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCa                      >  1:732084/1‑251 (MQ=255)                                                                                                                        aTTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       >  1:310059/1‑201 (MQ=255)                                                                                                                        aTTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCaa                     <  1:519184/251‑1 (MQ=255)                                                                                                                             gggCAATATGCGATGTTTTTCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGc                                                <  1:436411/219‑1 (MQ=255)                                                                                                                             gggCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAAt                >  1:668204/1‑251 (MQ=255)                                                                                                                               gCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATg               <  1:493687/250‑1 (MQ=255)                                                                                                                               gCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGa              >  1:225965/1‑251 (MQ=255)                                                                                                                                 aaTATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:771067/192‑1 (MQ=255)                                                                                                                                   tatGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:218485/190‑1 (MQ=255)                                                                                                                                    atgtgATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAggg                                                                                       <  1:443540/169‑1 (MQ=255)                                                                                                                                    atGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATGCCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat                                                                       <  1:868606/189‑1 (MQ=255)                                                                                                                                         aTGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTAtttt                                                                                           >  1:706359/1‑164 (MQ=255)                                                                                                                                           gTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAggg                                                                  <  1:784933/187‑1 (MQ=255)                                                                                                                                           gTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGACGAGCGAATggg  >  1:333872/1‑251 (MQ=255)                                                                                                                                                  |                                                                                                                                                                                                                                                    CGGTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCCTCGTTTTCCATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAAGGCTTTTAGCCAATAAATGACGAGCGAATGGG  >  1/1633143‑1633531

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: ‑