Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	885,975	T→G	D125E (GAT→GAG)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	885,975	0	T	G	100.0%	153.1 / NA	53	D125E (GAT→GAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/29);  total (24/29)

ATGCCGTTTTCAGCATCCGCCGCACGCATGAACCGAAATATCAGGATCCAAGATACAAAATTTCCACGGAAAACATGGAAAAATTCATTTTAGCCAAAAATAACCTGGAATATTGTTTGTTACCGGAATTGGGCAAGCAATCGGATAAACGTGTTACCGCCCTTGAAAAACGGCTATTAAAAGATATGATCAATGACGAGTTGTCGAAGATTGTGGGCAAAGCGTCGGCGCGTACGATTTTTGACGATCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCACTTTCGCGATCGTCTCAACGAACGCAAACGCCGGCAAGGCTTGGCGGGCGACAGTCGTTCTGTTGAAAACCAAGTGCAAAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAA  >  1/885728‑886219                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     aTGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAg                                                                                                                                                                                                                                                   <  1:219933/251‑1 (MQ=255)     ctttttCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACAAGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGtctt                                                                                                                                                                                                                                               <  1:186546/249‑1 (MQ=255)       ttttCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGtcttc                                                                                                                                                                                                                                              <  1:22906/250‑1 (MQ=255)        tttCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTcttctt                                                                                                                                                                                                                                            >  1:78169/1‑251 (MQ=255)              cacccgccgCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTa                                                                                                                                                                                                                                           >  1:816875/4‑246 (MQ=255)                         cgcaTGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTa                                                                                                                                                                                                                                           >  1:536522/1‑235 (MQ=255)                         cgcaTGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAAttt                                                                                                                                                                                                                           >  1:195905/1‑251 (MQ=255)                         cgcaTGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAAttt                                                                                                                                                                                                                           >  1:869396/1‑251 (MQ=255)                                      atatCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATcaacaa                                                                                                                                                                                                               >  1:610534/1‑250 (MQ=255)                                              aTCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCAc                                                                                                                                                                                                      <  1:676342/251‑1 (MQ=255)                                                         aaaaTTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATaa                                                                                                                                                                                           <  1:823312/251‑1 (MQ=255)                                                           aaTTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTa                                                                                                                                                                                                                                           <  1:623237/201‑1 (MQ=255)                                                           aaTTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACa                                                                                                                                                                                         <  1:306521/251‑1 (MQ=255)                                                                                          ttAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACg                                                                                                                                                          <  1:49425/251‑1 (MQ=255)                                                                                                  aaaTAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCAtttt                                                                                                                                                   >  1:156926/1‑250 (MQ=255)                                                                                                   aaTAACTTGGAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTcg                                                                                                                                                 <  1:254921/251‑1 (MQ=255)                                                                                                           ggAATATTGTTTGTTGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACg                                                                                                                                                                                                                                         <  1:597034/155‑1 (MQ=255)                                                                                                                   gtttgttGCCGGAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCaacgaa                                                                                                                                 >  1:347133/1‑251 (MQ=255)                                                                                                                              gAATTGGGCAAGCAATCGGACGAACGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCgcc                                                                                                                                         <  1:827179/232‑3 (MQ=255)                                                                                                                                                gacgaaCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGc                                                                                                                                                                        >  1:239162/5‑183 (MQ=255)                                                                                                                                                    aaCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGAt                                                                                                                                             <  1:119193/206‑1 (MQ=255)                                                                                                                                                    aaCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGt                                                                                                >  1:624029/1‑251 (MQ=255)                                                                                                                                                     aCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAAt                                                                                                                                                                          <  1:421098/176‑1 (MQ=255)                                                                                                                                                     aCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGc                                                                                                                                                                        <  1:129612/178‑1 (MQ=255)                                                                                                                                                     aCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGcctc                                                                                                                                       <  1:707499/211‑1 (MQ=255)                                                                                                                                                     aCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGt                                                                                                >  1:226126/1‑250 (MQ=255)                                                                                                                                                     aCGTGTTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGGGACAGt                                                                                                >  1:558556/1‑250 (MQ=255)                                                                                                                                                       gtgtTACTGCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGt                                                                                             >  1:562384/1‑251 (MQ=255)                                                                                                                                                               gCCCTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGa                                                                                     >  1:350079/1‑251 (MQ=255)                                                                                                                                                                cccTTGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAAc                                                                                                                                                                                          <  1:63035/149‑1 (MQ=255)                                                                                                                                                                   ttGAAAAACGCTTATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAAc                                                                                 <  1:826907/251‑1 (MQ=255)                                                                                                                                                                               tATTAAAAGATATGATCAATGACGAGTTTGCGAAGATTGTGGGCAAAACTTCGGCGCGCATGATTTTTGACGTGCAGACTTCTTACGATCATTTTCAATTTAAGAATCAACAACTCAATCACGGCATGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGCAAACGc                                                                                                                        >  1:27893/1‑200 (MQ=255)                                                                                                                                                                               tATTAAAAGATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGc                                                                                                                                                                        <  1:427334/152‑1 (MQ=255)                                                                                                                                                                                       gATATGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACa                                                             >  1:14135/1‑251 (MQ=255)                                                                                                                                                                                        atatGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCa                                                                               <  1:90385/232‑1 (MQ=255)                                                                                                                                                                                        atatGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACaa                                                            <  1:634074/251‑1 (MQ=255)                                                                                                                                                                                          atGATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACAt                                                                                                                                                                                        >  1:465626/1‑125 (MQ=255)                                                                                                                                                                                            gATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGa                                                        >  1:839317/1‑251 (MQ=255)                                                                                                                                                                                            gATCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGa                                                        >  1:797376/1‑251 (MQ=255)                                                                                                                                                                                             aTCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGAt                                                       >  1:796636/1‑251 (MQ=255)                                                                                                                                                                                              tCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCAc                                                                                                                                                                                                      <  1:422230/107‑1 (MQ=255)                                                                                                                                                                                              tCAATGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATa                                                      >  1:652731/1‑251 (MQ=255)                                                                                                                                                                                                 aTGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGAc                                                   <  1:561636/251‑1 (MQ=255)                                                                                                                                                                                                 aTGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGAc                                                   <  1:839783/251‑1 (MQ=255)                                                                                                                                                                                                  tGACGAGTTGGCGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACg                                                  <  1:696428/251‑1 (MQ=255)                                                                                                                                                                                                         ttggcGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGt                                            <  1:379285/246‑1 (MQ=255)                                                                                                                                                                                                           ggcGAAGATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAAc                                                                                                                          <  1:848133/168‑1 (MQ=255)                                                                                                                                                                                                                 gATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGAc                                                   >  1:849452/1‑235 (MQ=255)                                                                                                                                                                                                                 gATTGTGGGCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGtt                                   >  1:691709/1‑251 (MQ=255)                                                                                                                                                                                                                       gggCAAACCATCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCaacga                                                                                                                                  >  1:863529/1‑150 (MQ=255)                                                                                                                                                                                                                                atCGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACgag                    <  1:169170/250‑1 (MQ=255)                                                                                                                                                                                                                                  cGGCGCGCACGATTTTTGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGaa                  <  1:93155/251‑1 (MQ=255)                                                                                                                                                                                                                                                  tGACGAGCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCATTTTCGCGATCGCCTCAACGAACGTAAACGCCGTCAAGGCTTGGCGGGTGACAGTCGTTCCGTTGAAAACCAAATGCAGAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCaaa  <  1:103833/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     ATGCCGTTTTCAGCATCCGCCGCACGCATGAACCGAAATATCAGGATCCAAGATACAAAATTTCCACGGAAAACATGGAAAAATTCATTTTAGCCAAAAATAACCTGGAATATTGTTTGTTACCGGAATTGGGCAAGCAATCGGATAAACGTGTTACCGCCCTTGAAAAACGGCTATTAAAAGATATGATCAATGACGAGTTGTCGAAGATTGTGGGCAAAGCGTCGGCGCGTACGATTTTTGACGATCAGTCTTCTTACGATTATTTTCAATTTAAGTATCAACAACTCAATCACGGTGTGAATAACATTCGTGAATTTGAATGCCGCGATTTGAAACGGCACTTTCGCGATCGTCTCAACGAACGCAAACGCCGGCAAGGCTTGGCGGGCGACAGTCGTTCTGTTGAAAACCAAGTGCAAAAGCAACGACAAACGATAGACGATGCGTTTCGCCGTTTTCACACGCAACGAGAATCTATTGAACGCCAAA  >  1/885728‑886219

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