Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,081,329	A→G	F401F (TTT→TTC)	unknown ←	Membrane alanine aminopeptidase N (EC 3.4.11.2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,081,329	0	A	G	100.0%	180.6 / NA	65	F401F (TTT→TTC)	unknown	Membrane alanine aminopeptidase N (EC 3.4.11.2)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/40);  total (25/40)

CATTTGGTCGGCGGTCGGTGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGAAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGAAAAACCGTCAAGCC  >  1/1081101‑1081568                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                cATTTGGTCGGCGGTCGGTGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAgg                                                                                                                                                                                                                           <  1:417884/251‑1 (MQ=255) cATTTGGTCGGCGGTCGGTGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAgg                                                                                                                                                                                                                           <  1:552629/251‑1 (MQ=255)                cggtGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGaa                                                                                                                                                                                                            <  1:96373/251‑1 (MQ=255)                cggtGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGaa                                                                                                                                                                                                            <  1:247846/251‑1 (MQ=255)                cggtGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGaa                                                                                                                                                                                                            <  1:603180/251‑1 (MQ=255)                        ttAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTc                                                                                                                                                                                                     <  1:336144/250‑1 (MQ=255)                           aTTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGACCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGAACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTgcg                                                                                                                                                                                                 >  1:881400/1‑251 (MQ=255)                              gagaAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTgcgc                                                                                                                                                                                                >  1:617561/1‑249 (MQ=255)                                          aTTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGaa                                                                                                                                                                                                            >  1:712120/1‑225 (MQ=255)                                           ttGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACgg                                                                                                                                                                                 <  1:263484/251‑1 (MQ=255)                                              aTACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTc                                                                                                                                                                               <  1:420573/250‑1 (MQ=255)                                               tACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCAc                                                                                                                                                                             <  1:119959/251‑1 (MQ=255)                                                 cacaTGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCg                                                                                                                                                                           <  1:561987/251‑1 (MQ=255)                                                       aTTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGa                                                                                                                                                                     >  1:834183/1‑251 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:80591/251‑1 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:791359/251‑1 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:83561/251‑1 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:57860/251‑1 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:501980/251‑1 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:459885/251‑1 (MQ=255)                                                         ttttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:164865/251‑1 (MQ=255)                                                          tttCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTGCCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCAATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAgg                                                                                                                                                                                                                           <  1:274127/194‑1 (MQ=255)                                                               cGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTc                                                                                                                                                             <  1:891941/251‑1 (MQ=255)                                                                  cATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCTAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGa                                                                                                                                                                     <  1:428083/240‑1 (MQ=255)                                                                  cATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGgtg                                                                                                                                                                                                                         <  1:173767/188‑1 (MQ=255)                                                                          tCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACCCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   >  1:734648/1‑234 (MQ=255)                                                                                       aTTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGc                                                                                                                                                                                                                      >  1:258461/1‑170 (MQ=255)                                                                                       aTTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGaa                                                                                                                                                                                                            >  1:692143/1‑180 (MQ=255)                                                                                       aTTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGa                                                                                                                                     >  1:87067/1‑251 (MQ=255)                                                                                       aTTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGa                                                                                                                                     >  1:177398/1‑251 (MQ=255)                                                                                                tACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGAtg                                                                                                                             <  1:649237/250‑1 (MQ=255)                                                                                                  ccGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAAc                                                                                                                                                                                                                   <  1:130104/162‑1 (MQ=255)                                                                                                  ccGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   <  1:451545/210‑1 (MQ=255)                                                                                                  ccGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGAtgtg                                                                                                                           >  1:375265/1‑250 (MQ=255)                                                                                                  ccGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGAtg                                                                                                                             <  1:333069/248‑1 (MQ=255)                                                                                                         ggCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTgc                                                                                                                                                                                                  <  1:230727/172‑1 (MQ=255)                                                                                                         ggCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCC‑TTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   >  1:212983/1‑202 (MQ=255)                                                                                                          gCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCggg                                                                                                                   <  1:283628/250‑1 (MQ=255)                                                                                                          gCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCggg                                                                                                                   <  1:781984/250‑1 (MQ=255)                                                                                                           cTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCCTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGc                                                                                                                  <  1:447935/250‑1 (MQ=255)                                                                                                                    cGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTgcgc                                                                                                                                                                                                <  1:422965/163‑1 (MQ=255)                                                                                                                              gggTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGaaa                                                                                                                                                                   >  1:679094/1‑182 (MQ=255)                                                                                                                              gggTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTg                                                                                               >  1:487586/1‑250 (MQ=255)                                                                                                                                     aTCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGaa                                                                                                                                                                                                            >  1:5504/1‑134 (MQ=255)                                                                                                                                          aCTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  <  1:822203/251‑1 (MQ=255)                                                                                                                                          aCTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  <  1:230825/251‑1 (MQ=255)                                                                                                                                          aCTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  <  1:250824/251‑1 (MQ=255)                                                                                                                                          aCTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  <  1:552479/251‑1 (MQ=255)                                                                                                                                             cATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  <  1:358371/248‑1 (MQ=255)                                                                                                                                             cATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGatt                                                                               >  1:374024/1‑249 (MQ=255)                                                                                                                                             cATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAAt                                                                               >  1:717199/1‑251 (MQ=255)                                                                                                                                                    aCGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCggg                                                                                                                   >  1:863260/1‑208 (MQ=255)                                                                                                                                                       ttCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTc                                                                                                                                                             <  1:722076/163‑1 (MQ=255)                                                                                                                                                       ttCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  <  1:192712/238‑1 (MQ=255)                                                                                                                                                       ttCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGttat                                                                     >  1:350694/1‑251 (MQ=255)                                                                                                                                                                                       ccgcTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGa                                                                                                                               >  1:28916/1‑161 (MQ=255)                                                                                                                                                                                                       ttCCGCGATGTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGAt                     <  1:553169/251‑1 (MQ=255)                                                                                                                                                                                                                gTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCGTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGaaaaa            >  1:892787/1‑251 (MQ=255)                                                                                                                                                                                                                gTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAg                                                                                  >  1:330642/1‑181 (MQ=255)                                                                                                                                                                                                                gTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTg                       >  1:175730/1‑240 (MQ=255)                                                                                                                                                                                                                gTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGaaaaa            >  1:59989/1‑251 (MQ=255)                                                                                                                                                                                                                gTATAACTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGaaaaa            >  1:288666/1‑251 (MQ=255)                                                                                                                                                                                                                      cTTCATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGAAAAACCGTCa      >  1:820314/1‑251 (MQ=255)                                                                                                                                                                                                                         cATGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGAAAAACCGTCAAGc   <  1:34869/251‑1 (MQ=255)                                                                                                                                                                                                                          aTGCCTTTTTGGAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGAAAAACCGTCAAGcc  <  1:841475/251‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                CATTTGGTCGGCGGTCGGTGGAGTTAATTGAGAAACCTGCAATTGATACACATGATTTTTCTCGTCATAACGATCGCTGATGGTCAATTCAGGCGTACCGGATTGGCTGTACCAACGGCGAAATTGGGTTAAATCCAACTCATTGGCACGTTCCATCGCAGAAACAAAATCTTCACAGGTCGCCGCTTTGCCGTCATTTTCCGCGATGTATAACTTCATGCCTTTTTGAAAGCCTTTTTCGCCCAATAAGGTGTGCAACATACGAATGACTTCTGCGCCTTTTTCATACACGGTCACCGTATAGAAATTATTCATTTCAATCACTTTTTCCGGACGAATCGGATGTGCCATCGGGCCGGCATCCTCGGCGAATTGCACGGTGCGTAAGAATTTCACGTTATTAATGCGATTTACCGGGCGTGAGCCGGTGTCGGAGGAAAATTCCTGATCGCGAAAAACCGTCAAGCC  >  1/1081101‑1081568

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