Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	910,413	G→T	R74R (CGG→CGT)	unknown →	phage‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	910,413	0	G	T	100.0%	164.0 / NA	61	R74R (CGG→CGT)	unknown	phage‑related protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (32/29);  total (32/29)

CGGAACAACAGGAACACCACAATGACAACAAGAGTAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACCGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATCTGGAACATGGCTTGCCGTGGCTTGAACAAATGGGGCGCGGTGTGAATCTGGCGGATTGGGAAATTAAAATAAAGCGGCATGTACTGGGAACGGAAGGAGTGTCAAAAATTACCGATTATCAGGCGGATTTTGAACCCGATTCGCGCAAGCTGACAATCTCAATTGACTATCAAGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAAC  >  1/910171‑910630                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                          cggAACAACAGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgt                                                                                                                                                                                                                    <  1:391640/250‑3 (MQ=255)       aacaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagac                                                                                                                                                                                                             <  1:538397/251‑10 (MQ=255)       aacaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagac                                                                                                                                                                                                             <  1:528246/251‑10 (MQ=255)        acaGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggaga                                                                                                                                                                                                              >  1:301100/1‑241 (MQ=255)         caGGAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggaga                                                                                                                                                                                                              >  1:771209/1‑240 (MQ=255)            gAACACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagactgaag                                                                                                                                                                                                        >  1:880160/1‑237 (MQ=255)             aaCACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagactgaagg                                                                                                                                                                                                       <  1:484478/251‑16 (MQ=255)             aaCACCACAATGACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGtgttggagactgaagg                                                                                                                                                                                                       <  1:381610/251‑16 (MQ=255)                        gACAACAAGAATAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCaaaaa                                                                                                                                                                                            >  1:400948/1‑251 (MQ=255)                                gaatAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGa                                                                                                                                                                                     >  1:636485/4‑250 (MQ=255)                                gaatAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGa                                                                                                                                                                                     >  1:707404/4‑250 (MQ=255)                                gaatAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTAGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGATCGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGTCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGa                                                                                                                                                                                     >  1:545062/4‑250 (MQ=255)                                                   aaaTCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAAc                                                                                                                                                                 <  1:186166/251‑1 (MQ=255)                                                    aaTCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAAcc                                                                                                                                                                <  1:334215/251‑1 (MQ=255)                                                                         cAAGGTTTTTCGAACTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACaa                                                                                                                                            <  1:653799/250‑1 (MQ=255)                                                                                     aaCTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTCTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGatt                                                                                                                                <  1:292221/250‑3 (MQ=255)                                                                                       cTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAAc                                                                                                                                                                 <  1:850259/215‑1 (MQ=255)                                                                                       cTACGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTAt                                                                                                                              >  1:828515/1‑250 (MQ=255)                                                                                         aCGCCACTGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGaa                                                                                                                                                                  <  1:642525/212‑1 (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:342967/1‑251 (MQ=255)                                                                                                           cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGcaacaa                                                                                                         >  1:30810/1‑251 (MQ=255)                                                                                                           cTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGcaacaa                                                                                                         >  1:466254/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)                                                                                                                 cacaGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAAc                                                                                                                                                                 >  1:655820/1‑189 (MQ=255)                                                                                                                 cacaGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTATCAGGCTGATTTTGAACCCGATTCGCGCAAACTGACAATATCCATTGATTATCAGGACATTTACGGGCAACAACAGAc                                                                                                    <  1:654092/250‑1 (MQ=255)                                                                                                                     gAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATTTAGAACATGGCTTGCCATGGCTTGAGCAAATGGGGCGCGGTGTCAATCTGGCAGATTGGGAAATCAAAATTAAGCGTCATGTGTTGGAGACTGAAGGAGTGTCAAAAATCACCGACTAt                                                                                                                                                                                 <  1:491850/169‑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)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                          CGGAACAACAGGAACACCACAATGACAACAAGAGTAAGACGGCTGGATAAAAATCACGATTGGACGTTTGGGCAAGGTTTTTCGAACTACGCCACCGAATCTGAAGCTATCGCACAGAATGTACAAACCCGCCTTTGGTCGTTTGCCAATGACTGGTTTCTAGATCTGGAACATGGCTTGCCGTGGCTTGAACAAATGGGGCGCGGTGTGAATCTGGCGGATTGGGAAATTAAAATAAAGCGGCATGTACTGGGAACGGAAGGAGTGTCAAAAATTACCGATTATCAGGCGGATTTTGAACCCGATTCGCGCAAGCTGACAATCTCAATTGACTATCAAGACATTTACGGGCAACAACAGACCGCACGCTATACAGCTTAAAGTGCGGTCAAATCGACCGCATTTTTTAATGAACCGATAAATTAATTTCGCGCCCTAATTTTTTAAAGGCTTTTTCAAC  >  1/910171‑910630

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