Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	910,485	G→A	K98K (AAG→AAA)	unknown →	phage‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	910,485	0	G	A	100.0%	181.5 / NA	64	K98K (AAG→AAA)	unknown	phage‑related protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/27);  total (37/27)

CAAGGTTTTTCGAACTACGCCACCGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATCTGGAACATGGCTTGCCGTGGCTTGAACAAATGGGGCGCGGTGTGAATCTGGCGGATTGGGAAATTAAAATAAAGCGGCATGTACTGGGAACGGAAGGAGTGTCAAAAATTACCGATTATCAGGCGGATTTTGAACCCGATTCGCGCAAGCTGACAATCTCAATTGACTATCAAGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGGCGCGCCATTTCAGCTTTTCGCACAC  >  1/910243‑910725                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 cAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACaa                                                                                                                                                                                                                                           <  1:653799/250‑1 (MQ=255)             aaCTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTCTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGatt                                                                                                                                                                                                                               <  1:292221/250‑3 (MQ=255)               cTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTAt                                                                                                                                                                                                                             >  1:828515/1‑250 (MQ=255)                 aCGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCag                                                                                                                                                                                                                          >  1:288930/1‑250 (MQ=255)                          aaTCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTa                                                                                                                                                                                                                  <  1:463689/250‑1 (MQ=255)                           aTCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACg                                                                                                                                                                                                                <  1:242264/251‑1 (MQ=255)                                   cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGatt                                                                                                                                                                                                                               <  1:173288/228‑3 (MQ=255)                                   cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGcaacaa                                                                                                                                                                                                        >  1:466254/1‑251 (MQ=255)                                   cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGcaacaa                                                                                                                                                                                                        >  1:342967/1‑251 (MQ=255)                                   cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGcaacaa                                                                                                                                                                                                        >  1:30810/1‑251 (MQ=255)                                   cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTC‑AAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGcaacaa                                                                                                                                                                                                        >  1:187175/1‑250 (MQ=255)                                    tATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCaacaac                                                                                                                                                                                                       <  1:296157/251‑1 (MQ=255)                                     aTCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAacaaca                                                                                                                                                                                                      >  1:184280/1‑251 (MQ=255)                                     aTCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAacaaca                                                                                                                                                                                                      >  1:399625/1‑251 (MQ=255)                                         cacaGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGAc                                                                                                                                                                                                   <  1:654092/250‑1 (MQ=255)                                             gAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGAc                                                                                                                                                                                                                       >  1:30634/1‑226 (MQ=255)                                             gAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAg                                                                                                                                                                                                     >  1:95316/1‑244 (MQ=255)                                             gAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCa                                                                                                                                                                                               >  1:273577/1‑250 (MQ=255)                                               aTGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAAGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACg                                                                                                                                                                                             >  1:423446/1‑250 (MQ=255)                                                 gTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGAc                                                                                                                                                                                                                       >  1:601833/1‑222 (MQ=255)                                                 gTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGAc                                                                                                                                                                                                   <  1:825599/242‑1 (MQ=255)                                                 gTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCta                                                                                                                                                                                          >  1:838375/1‑251 (MQ=255)                                                    cAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATAc                                                                                                                                                                                       <  1:422224/251‑1 (MQ=255)                                                       aCCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGc                                                                                                                                                                                    <  1:872055/251‑1 (MQ=255)                                                       aCCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGc                                                                                                                                                                                    <  1:537712/251‑1 (MQ=255)                                                           gCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATAc                                                                                                                                                                                       <  1:397764/244‑1 (MQ=255)                                                             cTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAg                                                                                                                                                                              <  1:746853/251‑1 (MQ=255)                                                             cTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAg                                                                                                                                                                              <  1:677169/251‑1 (MQ=255)                                                                         gCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAAt                                                                                                                                                                   >  1:191585/1‑250 (MQ=255)                                                                          ccAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATc                                                                                                                                                                  <  1:568297/250‑1 (MQ=255)                                                                                     tttCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCAttttt                                                                                                                                                      >  1:693062/1‑251 (MQ=255)                                                                                           gatttagAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGa                                                                                                                                                >  1:106774/7‑251 (MQ=255)                                                                                                 gAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGAt                                                                                                                                          >  1:500570/1‑251 (MQ=255)                                                                                                           tGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAAttt                                                                                                                                <  1:711722/251‑1 (MQ=255)                                                                                                           tGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAAttt                                                                                                                                <  1:453999/251‑1 (MQ=255)                                                                                                            gccatGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAAttt                                                                                                                                >  1:391621/5‑250 (MQ=255)                                                                                                                 ggCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATAc                                                                                                                                                                                       <  1:129661/190‑1 (MQ=255)                                                                                                                 ggCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGc                                                                                                                                                                                    >  1:692418/1‑193 (MQ=255)                                                                                                                 ggCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCgcgc                                                                                                                           >  1:659441/1‑250 (MQ=255)                                                                                                                      gagcAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAAt                                                                                                                     >  1:177733/4‑251 (MQ=255)                                                                                                                       agcAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAAtt                                                                                                                    >  1:59655/3‑251 (MQ=255)                                                                                                                         cAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAAtttt                                                                                                                  >  1:543249/1‑251 (MQ=255)                                                                                                                                          caaTCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTcaac                                                                                                 <  1:236200/250‑1 (MQ=255)                                                                                                                                                                            catgtgttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCtaat                                                               >  1:550749/16‑251 (MQ=255)                                                                                                                                                                               gtgttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATggg                                                                              <  1:864945/221‑1 (MQ=255)                                                                                                                                                                               gtgttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGa                                                                        <  1:857766/227‑1 (MQ=255)                                                                                                                                                                               gtgttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAg                                                            >  1:207534/13‑251 (MQ=255)                                                                                                                                                                               gtgttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAg                                                            >  1:745130/13‑251 (MQ=255)                                                                                                                                                                               gtgttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAg                                                            >  1:699008/13‑251 (MQ=255)                                                                                                                                                                                 gttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATAc                                                                                                                                                                                       >  1:857539/11‑126 (MQ=255)                                                                                                                                                                                 gttggagactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAAt                                                                    >  1:524064/11‑241 (MQ=255)                                                                                                                                                                                       gactgAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGAt                                                                                 >  1:302231/5‑222 (MQ=255)                                                                                                                                                                                               gAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATAc                                                                                                                                                                                       <  1:333726/112‑1 (MQ=255)                                                                                                                                                                                                          atcaCCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCa                                 >  1:134219/4‑251 (MQ=255)                                                                                                                                                                                                                gactATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAAGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGAc                             >  1:451873/4‑249 (MQ=255)                                                                                                                                                                                                                  ctATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATc                                  <  1:827053/241‑1 (MQ=255)                                                                                                                                                                                                                  ctATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGgcg                         >  1:337650/2‑251 (MQ=255)                                                                                                                                                                                                                  ctATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGgcg                         >  1:577265/2‑251 (MQ=255)                                                                                                                                                                                                                             aTTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGGCGCGCCATTTCAg              >  1:191091/1‑251 (MQ=255)                                                                                                                                                                                                                              ttttGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATAc                                                                                                                                                                                       <  1:45263/81‑1 (MQ=255)                                                                                                                                                                                                                              ttttGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGGCGCGCCATTTCAGc             <  1:651485/251‑1 (MQ=255)                                                                                                                                                                                                                                        gATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATAAAGACGGCGCGCCATTTCAGCTTTTCGcaca   >  1:375028/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGGCGCGCCATTTCAGCTTTTCGCacac  <  1:720862/251‑1 (MQ=255)                                                                                                                                                                                                                                         aTTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGGCGCGCCATTTCAGCTTTTCGCacac  <  1:886721/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 CAAGGTTTTTCGAACTACGCCACCGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATCTGGAACATGGCTTGCCGTGGCTTGAACAAATGGGGCGCGGTGTGAATCTGGCGGATTGGGAAATTAAAATAAAGCGGCATGTACTGGGAACGGAAGGAGTGTCAAAAATTACCGATTATCAGGCGGATTTTGAACCCGATTCGCGCAAGCTGACAATCTCAATTGACTATCAAGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAACAAAATCAATTTTGGATGGGTGGCGAAAATCTAATAAGCGATCGACTTGCGGCATGTGAACATCAAGACGGCGCGCCATTTCAGCTTTTCGCACAC  >  1/910243‑910725

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