Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,339,136	T→C	A21A (GCA→GCG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,339,136	0	T	C	100.0%	204.3 / NA	74	A21A (GCA→GCG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/36);  total (38/36)

CTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCTTGCGCAAATCTTTCATTACATCTTCATGTGAGTAAGTCTTAATTTTTCCATTTAAGACACCTTGTAAAAAGCGCTGCTCTTCTGTTAAGGATTGCTCAAACATTGATAAAACTCTTTCTTGAAGATGAGTTTGCCCATTTTTTGATTTAGCTTGATTTAACATATT  >  1/1338918‑1339377                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                                                                                                                                                                                                    >  1:604529/4‑250 (MQ=255) tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                                                                                                                                                                                                    >  1:545586/4‑250 (MQ=255)      aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                                                                                                                                                                                              <  1:87745/251‑1 (MQ=255)      aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                                                                                                                                                                                              <  1:425820/251‑1 (MQ=255)        gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                                                                                                                                                                                            >  1:93247/1‑251 (MQ=255)        gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                                                                                                                                                                                            >  1:629036/1‑251 (MQ=255)                     tGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGC‑GCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatc                                                                                                                                                                                              <  1:193324/251‑4 (MQ=255)                                cGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCAt                                                                                                                                                                                     <  1:657619/250‑1 (MQ=255)                                      tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa                                                                                                                                                                              >  1:744739/1‑251 (MQ=255)                                      tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa                                                                                                                                                                              >  1:400564/1‑251 (MQ=255)                                       cATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTc                                                                                                                                                                                                                   >  1:792671/1‑213 (MQ=255)                                        ataataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatcc                                                                                                                                                                                             >  1:39444/1‑230 (MQ=255)                                        ataataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAAt                                                                                                                                                                            >  1:491142/1‑251 (MQ=255)                                         taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACaacgt                                                                                                                                                                                                  >  1:592328/1‑225 (MQ=255)                                           atatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGc                                                                                                                                                                          >  1:623153/8‑250 (MQ=255)                                            tatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCt                                                                                                                                                                        <  1:665540/245‑1 (MQ=255)                                                                 aCTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAc                                                                                                                                                                                                               >  1:518623/1‑191 (MQ=255)                                                                  cttcaCTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTAcct                                                                                                                                                   >  1:857298/5‑248 (MQ=255)                                                                        tGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg                                                                                                                                            <  1:410129/251‑1 (MQ=255)                                                                        tGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg                                                                                                                                            <  1:421895/251‑1 (MQ=255)                                                                          aaTGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg                                                                                                                                          <  1:277073/251‑1 (MQ=255)                                                                          aaTGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg                                                                                                                                          <  1:228148/251‑1 (MQ=255)                                                                                cATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCAt                                                                                                                                                        >  1:552000/1‑231 (MQ=255)                                                                                 aTTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATCTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTaggt                                                                                                                                   >  1:450901/1‑248 (MQ=255)                                                                                  tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTc                                                                                                                                  <  1:669865/251‑1 (MQ=255)                                                                                  tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTc                                                                                                                                  <  1:188738/251‑1 (MQ=255)                                                                                       aagtggaaaaaTGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAAt                                                                                                                             <  1:10290/244‑1 (MQ=255)                                                                                       aagtggaaaaaTGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAAt                                                                                                                             <  1:840595/244‑1 (MQ=255)                                                                                                   gtgaaatctCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCAttcaa                                                                                                                 >  1:50896/9‑248 (MQ=255)                                                                                                         tctCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatcc                                                                                                                                                                                             <  1:176733/167‑5 (MQ=255)                                                                                                         tctCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg                                                                                                                                          <  1:467730/218‑1 (MQ=255)                                                                                                            cccGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTg                                                                                                        >  1:826244/1‑251 (MQ=255)                                                                                                              cGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTa                                                                                                      >  1:508055/1‑251 (MQ=255)                                                                                                                 cagTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGc                                                                                                                                                                          >  1:144242/3‑180 (MQ=255)                                                                                                                 cagTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGac                                                                                                              <  1:109825/238‑1 (MQ=255)                                                                                                                  agTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTaaaaa                                                                                                  >  1:2064/2‑251 (MQ=255)                                                                                                                   gTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGTGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAg                                                                                                 >  1:504438/1‑251 (MQ=255)                                                                                                                   gTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAg                                                                                                 >  1:790983/1‑251 (MQ=255)                                                                                                                   gTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAg                                                                                                 >  1:883757/1‑251 (MQ=255)                                                                                                                    tataGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTaaaa                                                                                                   <  1:265917/248‑1 (MQ=255)                                                                                                                        ggcgtCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg                                                                                                                                            <  1:821727/199‑1 (MQ=255)                                                                                                                                  atgtaACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCaa                                                                                                                                                                             <  1:859494/157‑1 (MQ=255)                                                                                                                                      aaCTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAgg                                                                              <  1:839028/251‑1 (MQ=255)                                                                                                                                                   tGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACac                                                                 <  1:267620/251‑2 (MQ=255)                                                                                                                                                   tGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACac                                                                 <  1:282192/251‑2 (MQ=255)                                                                                                                                                        gacgctctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACac                                                                 <  1:24566/242‑2 (MQ=255)                                                                                                                                                            ctctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg                                                                                                                                          <  1:430895/169‑1 (MQ=255)                                                                                                                                                            ctctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACAcc                                                                                                           >  1:814858/1‑200 (MQ=255)                                                                                                                                                            ctctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATaatac                                                        <  1:686141/251‑4 (MQ=255)                                                                                                                                                            ctctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATaatac                                                        <  1:380722/251‑4 (MQ=255)                                                                                                                                                            ctctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATaata                                                         >  1:87181/1‑248 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAgg                                                                              >  1:560529/1‑227 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGac                                                                            <  1:142115/229‑2 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGac                                                                            <  1:166077/229‑2 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATActc                                                      <  1:690848/251‑1 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATActc                                                      <  1:848303/251‑1 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATActc                                                      <  1:776313/251‑1 (MQ=255)                                                                                                                                                              ctATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATAct                                                       >  1:345053/1‑250 (MQ=255)                                                                                                                                                                      tCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATaatac                                                        >  1:574803/1‑238 (MQ=255)                                                                                                                                                                         aTGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAg                                                                                                                                        >  1:135188/1‑158 (MQ=255)                                                                                                                                                                           gAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCtcttcag                                                                                    >  1:417789/1‑206 (MQ=255)                                                                                                                                                                           gAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGac                                                                            <  1:615472/216‑2 (MQ=255)                                                                                                                                                                           gAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTg                                               <  1:547591/245‑1 (MQ=255)                                                                                                                                                                               ttGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTGAAAGTGAGtt                                     >  1:847551/1‑251 (MQ=255)                                                                                                                                                                                 gTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatcc                                                                                                                                                                                             <  1:412830/97‑5 (MQ=255)                                                                                                                                                                                         cTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAAtt                                                                                                                            >  1:546373/1‑154 (MQ=255)                                                                                                                                                                                           tCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTGAAAGTGAGTTTGCTCAtttttt                         >  1:664728/1‑251 (MQ=255)                                                                                                                                                                                               tactttGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATaatac                                                        >  1:565035/4‑213 (MQ=255)                                                                                                                                                                                                     ggCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCAt                                                                                                                                                                                     >  1:439411/1‑85 (MQ=255)                                                                                                                                                                                                           gtgtgtAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTa                                                                                 <  1:7649/179‑1 (MQ=255)                                                                                                                                                                                                                tAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTGAAAGTGAGTTTGCTCATTTTTTAATTTAGCTTGATTTAACat     <  1:482085/250‑1 (MQ=255)                                                                                                                                                                                                                  aaTATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTGAAAGTGAGTTTGCTCATTTTTTAATTTAGCTTGATTTAACat     >  1:667739/1‑248 (MQ=255)                                                                                                                                                                                                                  aaTATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTGAAAGTGAGTTTGCTCATTTTTTAATTTAGCTTGATTTAACATAtt  >  1:312984/1‑251 (MQ=255)                                                                                                                                                                                                                  aaTATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAATTTTGCCATTCAAGACACCTTGTAAAAAGCGTTGCTCTTCAGTTAAGGACTGTTCAAACACTGATAATACTCTTTCTTGAAAGTGAGTTTGCTCATTTTTTAATTTAGCTTGATTTAACATAtt  >  1:451611/1‑251 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  CTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCTTGCGCAAATCTTTCATTACATCTTCATGTGAGTAAGTCTTAATTTTTCCATTTAAGACACCTTGTAAAAAGCGCTGCTCTTCTGTTAAGGATTGCTCAAACATTGATAAAACTCTTTCTTGAAGATGAGTTTGCCCATTTTTTGATTTAGCTTGATTTAACATATT  >  1/1338918‑1339377

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