Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,942,424	A→G	R175R (CGT→CGC)	unknown ←	Tail‑specific protease precursor (EC 3.4.21.102)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,942,424	0	A	G	100.0%	250.3 / NA	75	R175R (CGT→CGC)	unknown	Tail‑specific protease precursor (EC 3.4.21.102)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/38);  total (37/38)

GTGGCACCGACCCCTTCTAAAGATAAATTCATGCTTTCATTGAAACTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACACGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAATTACCTGCATTCAACTCGTCATCAATTTTATTGCCATATTT  >  1/1942179‑1942648                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 gtgGCACCGACCCCTTCTAAAGATAAATTCATGCTTTCATTGAAACTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAAcgcgcg                                                                                                                                                                                                                              >  1:400072/1‑250 (MQ=255)           ccccTTCTAAAGATAAATTCATGCTTTCATTGAAACTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAAc                                                                                                                                                                                                                   <  1:50874/251‑1 (MQ=255)              cTTCTAAAGATAAATTCATGCTTTCATTGAAACTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAAc                                                                                                                                                                                                                   >  1:895982/1‑248 (MQ=255)                tCTAAAGATAAATTCATGCTTTCATTGAAACTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTAtt                                                                                                                                                                                                              >  1:769873/1‑251 (MQ=255)                                              cTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGAGAATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACAACATCATCCGCTTTAGTTTGGGTTAAGCGAGGTATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGCTGTCGTTTTTAACGCGAGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCTAgg                                                                                                                                                                                    >  1:769890/1‑247 (MQ=255)                                                     ccGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCa                                                                                                                                                                         <  1:809295/251‑1 (MQ=255)                                                          cGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAAc                                                                                                                                                                                                                   >  1:466643/1‑204 (MQ=255)                                                              gggCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAAc                                                                                                                                                                                                                   >  1:274312/1‑200 (MQ=255)                                                                      tAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTg                                                                                                                                                        >  1:684509/1‑251 (MQ=255)                                                                         cTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGat                                                                                                                                                      >  1:537782/1‑250 (MQ=255)                                                                           ggTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGatcat                                                                                                                                                   >  1:563606/1‑251 (MQ=255)                                                                              aTGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGatcatcat                                                                                                                                                <  1:712447/251‑1 (MQ=255)                                                                                gtgGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCg                                                                                                                                              >  1:732680/1‑251 (MQ=255)                                                                                   gATCAATTTCACGGGCAAAAGCGTTTAATTAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATCTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCTCATAACTTCTTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGtt                                                                                                                                            >  1:365600/1‑250 (MQ=255)                                                                                   gATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGtt                                                                                                                                            >  1:676742/1‑250 (MQ=255)                                                                                           tCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCgg                                                                                                                                   <  1:389383/251‑1 (MQ=255)                                                                                                ggCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCtt                                                                                                                              >  1:740420/1‑251 (MQ=255)                                                                                                    aaaGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATc                                                                                                                          <  1:357886/251‑1 (MQ=255)                                                                                                     aaGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATcc                                                                                                                         <  1:758357/251‑1 (MQ=255)                                                                                                                  ggTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACtt                                                                                                                                                                                                                 <  1:639309/150‑1 (MQ=255)                                                                                                                    tttGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGc                                                                                                           <  1:24029/250‑1 (MQ=255)                                                                                                                       gTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGatcatcat                                                                                                                                                >  1:653071/1‑210 (MQ=255)                                                                                                                       gTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTAt                                                                                                                           <  1:325386/231‑1 (MQ=255)                                                                                                                       gTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTAt                                                                                                                           >  1:131329/1‑231 (MQ=255)                                                                                                                       gTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATcc                                                                                                                         >  1:560360/1‑233 (MQ=255)                                                                                                                       gTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataa                                                                                                       >  1:334844/1‑251 (MQ=255)                                                                                                                       gTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataa                                                                                                       >  1:406975/1‑251 (MQ=255)                                                                                                                        tACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCgtaac                                                                                                      <  1:636661/251‑1 (MQ=255)                                                                                                                         ccaccatcatCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTAT‑TTTAAATTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataacga                                                                                                    <  1:216277/247‑1 (MQ=255)                                                                                                                          ccACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataacga                                                                                                    <  1:216306/251‑1 (MQ=255)                                                                                                                               tcatcCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataacgataacg                                                                                               <  1:838765/251‑1 (MQ=255)                                                                                                                                catcCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGt                                                                                              <  1:114251/251‑1 (MQ=255)                                                                                                                                catcCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGt                                                                                              <  1:270434/251‑1 (MQ=255)                                                                                                                                   ccGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGAtcaatttcaatt                                                                                                                                                          <  1:592495/188‑1 (MQ=255)                                                                                                                                     gCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataac                                                                                                      >  1:797676/1‑238 (MQ=255)                                                                                                                                     gCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATa                                                                                         >  1:427407/1‑251 (MQ=255)                                                                                                                                     gCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATa                                                                                         >  1:45072/1‑251 (MQ=255)                                                                                                                                      cTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATaa                                                                                        >  1:46203/1‑251 (MQ=255)                                                                                                                                       tttAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAAc                                                                                       <  1:127645/251‑1 (MQ=255)                                                                                                                                          aGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGat                                                                                                                                                      <  1:368110/185‑1 (MQ=255)                                                                                                                                          aGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACgg                                                                                     >  1:632483/1‑250 (MQ=255)                                                                                                                                               gggTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataacga                                                                                                    >  1:161132/1‑230 (MQ=255)                                                                                                                                                ggTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTtca                                                                              >  1:689157/1‑251 (MQ=255)                                                                                                                                                    aaGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGc                                                                                                           <  1:901338/218‑1 (MQ=255)                                                                                                                                                         acgaATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACt                                                                     <  1:852834/251‑1 (MQ=255)                                                                                                                                                               cGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAAc                                                                                       <  1:86461/227‑1 (MQ=255)                                                                                                                                                                 ccAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTcatacat                                                              <  1:238604/250‑1 (MQ=255)                                                                                                                                                                 ccAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTcatacat                                                              <  1:386063/250‑1 (MQ=255)                                                                                                                                                                     aTTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataacga                                                                                                    <  1:666931/208‑1 (MQ=255)                                                                                                                                                                     aTTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAAc                                                                                       >  1:214373/1‑221 (MQ=255)                                                                                                                                                                      ttATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCAc                                                                                                                                                                        >  1:195147/1‑139 (MQ=255)                                                                                                                                                                      ttATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAAc                                                                                       >  1:880492/1‑220 (MQ=255)                                                                                                                                                                      ttATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTcatac                                                                <  1:875150/243‑1 (MQ=255)                                                                                                                                                                      ttATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTaaa                                                        >  1:392479/1‑251 (MQ=255)                                                                                                                                                                        ataACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTtcatca                                                                           <  1:336339/230‑1 (MQ=255)                                                                                                                                                                         taACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGc                                                     <  1:17392/251‑1 (MQ=255)                                                                                                                                                                         taACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGc                                                     <  1:386392/251‑1 (MQ=255)                                                                                                                                                                            cGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAg                                                      <  1:413578/247‑1 (MQ=255)                                                                                                                                                                            cGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAat                                                  >  1:709642/1‑251 (MQ=255)                                                                                                                                                                            cGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGATTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAat                                                  >  1:446275/1‑251 (MQ=255)                                                                                                                                                                              cTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGatcatcat                                                                                                                                                >  1:164400/1‑155 (MQ=255)                                                                                                                                                                                tAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCt                                              <  1:574314/251‑1 (MQ=255)                                                                                                                                                                                   cTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTaaa                                           <  1:813865/251‑1 (MQ=255)                                                                                                                                                                                    tAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAAt                                          <  1:230430/251‑1 (MQ=255)                                                                                                                                                                                          aTCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTCAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAg                                                      >  1:741593/1‑233 (MQ=255)                                                                                                                                                                                            cTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAATTACCTGCa                                  >  1:249072/1‑251 (MQ=255)                                                                                                                                                                                             tttAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGc                                                     <  1:617554/231‑1 (MQ=255)                                                                                                                                                                                                  tttCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAATTACCTGCATTCAAc                            <  1:537790/251‑1 (MQ=255)                                                                                                                                                                                                            ttttttATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATg                                                                                                                                      <  1:626326/135‑1 (MQ=255)                                                                                                                                                                                                            ttttttATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACgg                                                                                     >  1:259713/1‑184 (MQ=255)                                                                                                                                                                                                            ttttttATCTTTTAATTTCAAGCTGATGAGGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAAc                                                                      <  1:447877/199‑1 (MQ=255)                                                                                                                                                                                                                          aTTTCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTc                                                                    <  1:694970/187‑1 (MQ=255)                                                                                                                                                                                                                            ttCAAGCTGATGATGTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAATTACCTGCATTCAACTCGTCATCAATTTTATTGCCATAttt  <  1:83759/251‑1 (MQ=255)                                                                                                                                                                                                                                  cTGATGATGTCGTTTTTAACGCGCGCTTCCCCTAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGatcatcatc                                                                                                                                               >  1:723581/1‑104 (MQ=255)                                                                                                                                                                                                                                        atgTCGTTTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAATTACCTGc                                   >  1:307392/1‑206 (MQ=255)                                                                                                                                                                                                                                          gtcgtcTTTAACGCGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCataacgataa                                                                                                 >  1:127616/1‑142 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 GTGGCACCGACCCCTTCTAAAGATAAATTCATGCTTTCATTGAAACTTTTTGCCGTACGTGGGGCTAAATAACTGGTATGTGGATCAATTTCACGGGCAAAAGCGTTTAAATAGGTTTGTACCACATCATCCGCTTTAGTTTGGGTTAAGCGACGAATCGCCAAATTATAACGCTTAACTAATTTATCTTTAATTTCCGACCATTTTTTATCTTTTAATTTCAAGCTGATGATGTCGTTTTTAACACGCGCTTCCCATAACTTATTCGCTTCTTTTTCGGTTTTCGGCCAGGCTGCTTTTTCACGATCAATTTCAATTTGATCATCATCGTTCAAATGCGGTTCTTTATCCAGCAAAGACAAGGCATAACGATAACGTTCATAACGGCGTTTCATCATTAACTCATACATATTAAAAGCAATATCTAAATTACCTGCATTCAACTCGTCATCAATTTTATTGCCATATTT  >  1/1942179‑1942648

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