Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,339,049	A→G	Y50Y (TAT→TAC)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,339,049	0	A	G	98.0%	141.9 / ‑3.6	51	Y50Y (TAT→TAC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/21);  minor base C (1/0);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATGATACACACTTTCTCCTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCTTGCGCAAATCTTTCATTACATCTTCATGTGAGTAAGTCTTAATTTTTCCATTTAAGACACCTTGTAAAAAG  >  1/1338901‑1339282                                                                                                                                                     |                                                                                                                                                                                                                                          aaataacatacaccccttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTc                                                                                                                                                     >  1:818686/21‑235 (MQ=255)      acacacaccccttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTAt                                                                                                                                                                                                              <  1:448188/158‑1 (MQ=255)                 ttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAata                                                                                                                                                          >  1:225443/5‑214 (MQ=255)                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                                                                                                     >  1:604529/4‑250 (MQ=255)                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                                                                                                     >  1:545586/4‑250 (MQ=255)                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGATTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTaaa                                                                                                                                                            >  1:79663/4‑211 (MQ=255)                       aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                                                                                               <  1:87745/251‑1 (MQ=255)                       aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                                                                                               <  1:425820/251‑1 (MQ=255)                         gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCtt                                                                                                                                                                      >  1:491282/1‑194 (MQ=255)                         gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                                                                                             >  1:93247/1‑251 (MQ=255)                         gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                                                                                             >  1:629036/1‑251 (MQ=255)                                 ccTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTa                                                                                                                                                                                                               >  1:43191/1‑145 (MQ=255)                                      tGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGC‑GCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatc                                                                                               <  1:193324/251‑4 (MQ=255)                                                 cGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCAt                                                                                      <  1:657619/250‑1 (MQ=255)                                                       tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa                                                                               >  1:400564/1‑251 (MQ=255)                                                       tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa                                                                               >  1:744739/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)                                                          taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTccgc                                                                                                                                                  <  1:269537/181‑4 (MQ=255)                                                          taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACaacgt                                                                                                   >  1:592328/1‑225 (MQ=255)                                                            atatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGc                                                                           >  1:623153/8‑250 (MQ=255)                                                             tatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCt                                                                         <  1:665540/245‑1 (MQ=255)                                                                     cctacaaccaGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTtgtg                                                                                                                                                                  >  1:767337/4‑154 (MQ=255)                                                                     cctacaaccaGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAata                                                                                                                                                          <  1:488859/159‑1 (MQ=255)                                                                                  aCTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAc                                                                                                                >  1:518623/1‑191 (MQ=255)                                                                                   cttcaCTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTAcct                                                    >  1:857298/5‑248 (MQ=255)                                                                                         tGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg                                             <  1:421895/251‑1 (MQ=255)                                                                                         tGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg                                             <  1:410129/251‑1 (MQ=255)                                                                                           aaTGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg                                           <  1:228148/251‑1 (MQ=255)                                                                                           aaTGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg                                           <  1:277073/251‑1 (MQ=255)                                                                                                 cATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCAt                                                         >  1:552000/1‑231 (MQ=255)                                                                                                  aTTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATCTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTaggt                                    >  1:450901/1‑248 (MQ=255)                                                                                                   tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTac                                                                                                                                                                              >  1:144322/1‑111 (MQ=255)                                                                                                   tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTc                                   <  1:188738/251‑1 (MQ=255)                                                                                                   tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTc                                   <  1:669865/251‑1 (MQ=255)                                                                                                        aagtggaaaaaTGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAAt                              <  1:840595/244‑1 (MQ=255)                                                                                                        aagtggaaaaaTGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTCTTAAt                              <  1:10290/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)                                                                                                                                                     |                                                                                                                                                                                                                                          ATGATACACACTTTCTCCTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCTTGCGCAAATCTTTCATTACATCTTCATGTGAGTAAGTCTTAATTTTTCCATTTAAGACACCTTGTAAAAAG  >  1/1338901‑1339282

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