Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,501,500	C→T	L81L (CTG→TTG)	yjgZ →	uncharacterized protein YjgZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,501,500	0	C	T	97.4%	261.3 / NA	77	L81L (CTG→TTG)	yjgZ	uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/41);  minor base G (0/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAC  >  NC_000913/4501255‑4501742                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cccTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTgacg                                                                                                                                                                                                                                                <  1:216201/251‑1 (MQ=255)     tttGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGCCGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGACCCGCAGTCTGGAGCAACATGCTTGGTTgacggacg                                                                                                                                                                                                                                            >  2:216618/1‑251 (MQ=255)     gttGTTACCTTCTGGACCTCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGCGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTgacggacg                                                                                                                                                                                                                                            <  2:98736/250‑1 (MQ=255)         ttACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCt                                                                                                                                                                                                                                        >  2:108126/1‑251 (MQ=255)           aCCTTCTGGACCGCTTCGTGTCCTTCCAGGCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGAGCGTGCATTTCAAGGCACACACAACAATAGCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGGTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACTGACGTCCTGa                                                                                                                                                                                                                                      <  2:561588/251‑1 (MQ=255)                                   gctaGTCGGAGCGAGAGTTAGCCGGGGCGGACTATCACCGACGGAAAGAGCGTGCATTTCAAGGCACACACAACAAGATCTGCGAGCGTGCTATCCGGCCGGGCGTGATGGGACGAAAGGCCTGGTTGTTCGTTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGAGAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggag                                                                                                                                                                                                               <  1:586835/247‑1 (MQ=255)                                    ccAGTCGGAGCGAGAGTTATCCGGGGCGGACTTCCACCGACGGAAAGATCGGGCATTTCAAGGCACACACAACTATAGCGGCGAGCGTGCTATCCGTCAGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCGGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggagg                                                                                                                                                                                                              <  2:498923/250‑1 (MQ=255)                                    ccAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggagga                                                                                                                                                                                                             <  1:108126/251‑1 (MQ=255)                                    ccAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTGCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTCTCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCggagga                                                                                                                                                                                                             <  1:401350/251‑1 (MQ=255)                                          ggAGCGAGAGTTAGCCGGGGCGGACTTTCACCGCCGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGCGGAGAGGtt                                                                                                                                                                                                       >  2:176541/1‑251 (MQ=255)                                           gtgCGAGAGTTAGCCGGGGCGGAATTTCTCCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTTGTTGTTCTCGGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAAGGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTg                                                                                                                                                                                                      <  1:568355/249‑1 (MQ=255)                                           gAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTg                                                                                                                                                                                                      <  2:1438/251‑1 (MQ=255)                                                  agTTAGCCGGGGCGGACATTCACCGACGGAAAGATCGTGCATGTCAAGGCACACACAAAAAGATCTGCGAGCGTGCTAGCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTTTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGt                                                                                                                                                                                               <  2:149676/251‑1 (MQ=255)                                                      aGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTAc                                                                                                                                                                                            <  1:117593/250‑1 (MQ=255)                                                           tggCGGGCTTTCAACGCCGGACAGAACGTGCATCTCAAGGAACACCCAAAAATATCCGAGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGCCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCCCACGCTTGGTTGCCGGACGTCCTGACGCGTCTGCAGGAGTGGCCGGAGGAGCGGTTGGAAGAGTTACTGCctt                                                                                                                                                                                      >  2:546424/2‑250 (MQ=255)                                                             gCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCtt                                                                                                                                                                                    <  1:352269/251‑1 (MQ=255)                                                              cGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTg                                                                                                                                                                                   <  1:393872/251‑1 (MQ=255)                                                               ggACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGa                                                                                                                                                                                  <  2:311815/251‑1 (MQ=255)                                                               ggACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGa                                                                                                                                                                                  <  2:164128/251‑1 (MQ=255)                                                                  cgttCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTTGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAggg                                                                                                                                                                               <  1:341924/249‑1 (MQ=255)                                                                     tCACCGACGGAAAGATCGTTCATTTCAAGGCACACACAACAATATCTGCGAGCGGTCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACATATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGTAGTGGCCGGAGGAGATGTTGGCTGAGTTACTGCCTCTTGAGGGCt                                                                                                                                                                             <  2:384340/250‑1 (MQ=255)                                                                           aCGGAAAGATCGTTCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGGACAGATAATGAGCCTTCTGGGAACCGGAGTCTGGAGCCACGTGCTTGGTTGACGGACGTCCTGTCGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGGTACTGCCTCTTGAGGGCGTTACctt                                                                                                                                                                      <  2:203660/251‑1 (MQ=255)                                                                           aCGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCGGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCATATGCTTGGTTGGCGGACGTCCTGACGCGTCTGCCGGAGTGGCTGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACctt                                                                                                                                                                      <  2:336989/251‑1 (MQ=255)                                                                                 aGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCgg                                                                                                                                                                <  1:329920/251‑1 (MQ=255)                                                                                  gATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGCGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCggg                                                                                                                                                               >  1:83003/1‑251 (MQ=255)                                                                                   aTCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGt                                                                                                                                                              >  1:59374/1‑251 (MQ=255)                                                                                          aTTTCAGGGCACAGACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACCGATAATGAGCCTTCTGGGAACCGAAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGCGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTCCCTTCTTCGGGTAAGTGAt                                                                                                                                                       >  2:85888/1‑251 (MQ=255)                                                                                          aTTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGAt                                                                                                                                                       >  1:441586/1‑251 (MQ=255)                                                                                            ttCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATa                                                                                                                                                      <  2:128864/250‑1 (MQ=255)                                                                                              cAAGGCAAACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAAGCATGGTTGATCGCTGGTTCACTGGTGGCCGGGAACCGAGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGCGCCACATGATTGGTTGACGGACGTCCTGACGCGTCTGCCGGCGTGGCCGGAGGAGAGGGTGGCTGAGTTACTGCCTCTTGAGGGCTTTCCCTTCTTCGGGTAAGTGATACCt                                                                                                                                                   >  2:541045/1‑251 (MQ=255)                                                                                               acggCACACACAACAATATCTGCGCGCGTGCTCTCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTCCCTTCTTCGGGTAAGTGATACCTg                                                                                                                                                  >  2:572084/3‑251 (MQ=255)                                                                                               aaGGCACACACAACCATATCTGCGAGCGTGCTATCCGTCCGGGCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACGGCCTCTTGAGGGCTTTACCTGCTTCGGGTAAGTGATACCTg                                                                                                                                                  <  1:152853/251‑1 (MQ=255)                                                                                                 ggCACACACAACCATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGcc                                                                                                                                                >  1:436201/1‑251 (MQ=255)                                                                                                 ggCACACACAACAATATCTGCGATCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAGCCGCGCGGCACAGATAATGATCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGGTGACGTACGTCCTGACGCGTCTGCCGGGGTGGCCGGATGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAATTGATATCTGcc                                                                                                                                                <  2:334147/251‑1 (MQ=255)                                                                                                      acacaACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTGGACGGACGTCCTGACGCGTCGGCCGGAGTGGCCGGAGGAGAGGTGGGCTGAGTTACTGCCTCTTGAGGGGTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAg                                                                                                                                           <  2:579847/251‑1 (MQ=255)                                                                                                        acaacaATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGgt                                                                                                                                         <  1:216618/251‑1 (MQ=255)                                                                                                            caagatCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGttt                                                                                                                                     <  1:550244/247‑1 (MQ=255)                                                                                                            caaTATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGttt                                                                                                                                     <  1:14279/251‑1 (MQ=255)                                                                                                             aaTATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTtg                                                                                                                                    <  1:257987/251‑1 (MQ=255)                                                                                                                     cGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCgg                                                                                                                            >  1:103589/1‑251 (MQ=255)                                                                                                                       aGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGc                                                                                                                          >  1:322812/1‑251 (MQ=255)                                                                                                                         cGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGcc                                                                                                                         >  1:325308/1‑250 (MQ=255)                                                                                                                                  cgtccgGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGt                                                                                                                                              >  1:540281/1‑220 (MQ=255)                                                                                                                                  cgtccgGGAGTTAGGGGACGAAAGGCCGGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGAACAGATAAGGAGCCTTCTGGGAACCGCAGTCTGGAGCCACAGGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGTTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGt                                                                                                                                              <  2:540281/220‑1 (MQ=255)                                                                                                                                    tccgGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGGGTTTGTGCACCGGGCCATAACCTGCAGt                                                                                                             >  2:478651/1‑251 (MQ=255)                                                                                                                                     ccgGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCg                                                                                                                                                                                                                           >  1:207906/1‑140 (MQ=255)                                                                                                                                          cGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGa                                                                                                        <  1:441150/250‑1 (MQ=255)                                                                                                                                          cGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGa                                                                                                        <  2:500361/250‑1 (MQ=255)                                                                                                                                           gTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAAt                                                                                                      >  1:486098/1‑251 (MQ=255)                                                                                                                                           gTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGCTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAAt                                                                                                      >  1:275346/1‑251 (MQ=255)                                                                                                                                                gggACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAAc                                                                                                 >  1:96036/1‑251 (MQ=255)                                                                                                                                                    cGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCt                                                                                             >  1:582922/1‑251 (MQ=255)                                                                                                                                                          gCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGGTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGacaac                                                                                       <  1:567535/251‑1 (MQ=255)                                                                                                                                                          gCCTGGTTGTTCGCTGGTTCACTGGGGGCCGGGACCCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGCGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCGGacaac                                                                                       >  2:83504/1‑251 (MQ=255)                                                                                                                                                          gCCTGGTTGATCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGCACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGCGTGGCCGGAGGAGAGGTTGGCTGAGGTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGacaac                                                                                       >  2:198241/1‑251 (MQ=255)                                                                                                                                                             tGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAgg                                                                                    >  2:304559/1‑251 (MQ=255)                                                                                                                                                              ggTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGa                                                                                   <  1:176541/251‑1 (MQ=255)                                                                                                                                                               gttgttCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGa                                                                                   >  1:560364/1‑250 (MQ=255)                                                                                                                                                                       cTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGAGTCTGCCGGAGTGGCCGGCGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAAAGCCTGACCACAGGAAATGAGCt                                                                           >  2:579531/1‑250 (MQ=255)                                                                                                                                                                         gggtCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTag                                                                         <  2:417982/247‑1 (MQ=255)                                                                                                                                                                             cACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCTCATGCTGGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCGCTTGATGGCTTTACCTTCTTCGGGTAATTGCTACCTGCCGTCAGGTGGTTGTGCACCGGGCCATATCCTGCAGTTGGGAATTGAACGCCTGACAACATGAAATGAGCTAGAGCa                                                                     <  2:502473/250‑1 (MQ=255)                                                                                                                                                                                      ccGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAc‑gc                                                          <  2:322812/251‑1 (MQ=255)                                                                                                                                                                                      ccGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTGGTGCACCGGGCCATAACCTGCAGTGGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAc‑gc                                                          <  1:85888/251‑1 (MQ=255)                                                                                                                                                                                       ggggAACCGCGCGGAACAGATAATGAGCCTTCTGGGAACCGCGGTCTGGAGCCACAGGCTTGGTTGACGGCCGTCCTGACGGGTCGGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACGGCCTCGTGAGGGTTTTACCTTCTTCGGGTAAGGGATACCGGCCGTCAGGTGTTGGTGCCCCGGGCCATAACCTGCAGTTGGGAATGGAACGCCTGACAACGGGAAATGAGCTAGAGCAGTAAACAc‑gcg                                                         <  2:436201/250‑1 (MQ=255)                                                                                                                                                                                            aCCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCATGACGCGTACGCCGGAGTGGCAGGAGGAGAGGTTGGCTGAGTTACCGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGACGTCAGGGGTTTGTGCACCGGGACATAACCTGCAGTTGGGAATTGAACGCCTGCCAACAGGAAATGCGCTAGACCAGTAA‑CACCGCGGGCTg                                                    >  2:497255/1‑251 (MQ=255)                                                                                                                                                                                                     cacaGATAATGAGCCTTCTGGGATCCGCAGTCTGGAGCCTCATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTGGAGGGCTTTACCTTCTTCGGGTAAGTGATACCGGCCGTCAGGGGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATggtt                                           <  2:507615/251‑3 (MQ=255)                                                                                                                                                                                                     cacaGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATgggg                                           >  2:301979/1‑250 (MQ=255)                                                                                                                                                                                                              tGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAg                                                                                                                                           >  2:260116/1‑147 (MQ=255)                                                                                                                                                                                                              tGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAg                                                                                                                                           <  1:260116/147‑1 (MQ=255)                                                                                                                                                                                                                ggctgtctgggatccgccttgtggggCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGATGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGTTAAGTGATACCTGTCGTCAGGTGTTTGTTCACCGGGCCATAACCTGCAGTTGGGAAGTGAACGCGTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATGGGTCAGAATATGAg                                <  2:239165/226‑1 (MQ=255)                                                                                                                                                                                                                aGCCTTATGGGAACCGCAGTATCGAGCCACATGCTTGGTTGACGGACGTCATGACGAGTCTGCCGGAGTGGCAGGAGGAGAGGTAGGCTGAGTTCCTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGAAACCTGCCGTCAGGTGTTTGTGCACCGGGCCATACCCTGCAGTTGGGAATTGAACGCCTGACAAAAGGAAATGAGCTAGAGCAGTAAACTC‑GCGGGCTGACAATGGGTCAGAATATg                                  >  2:106474/1‑249 (MQ=255)                                                                                                                                                                                                                    ttCTGGGAACCGCAGTCTGGAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATGGGTCAGAATATGAGtt                              >  1:583154/1‑249 (MQ=255)                                                                                                                                                                                                                         gggaCCGCATTCGGGAACCACATGCTTGGTTGAAGGACGACCTGACGCGTCTGCCGGAGTGGCAGAAGGAGAGGTTGGCTGGGTTACTGCCTCTTGATAGCTTTACCTTCTTCGGGTAATGGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGGTGGGAATTTAAAGCATGACAGCAGGAAATGAGTTAGAGCCGGCAACAC‑GGGGGCTGACAATGGGACAGAATATGAGTGAGCCCgg                       >  2:418900/4‑251 (MQ=255)                                                                                                                                                                                                                          gAACCGTAGTCTGGAGCCCCAGGCTTGGGTGACGTGCTGCTTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGGTGAGTTACTGCCTCGTGAGGGCTTTCCCCTCTGCGGGTAAGGGATACCTGCCGTCAGGTGTTTTTGCACCGGGCGAGATCCTGCAGTTGGGAATTGAAGGCCTTACAATATGAAATAAGCGAGAGCAGTAAACAC‑GCGGGCTGGGAATGGGGCAGAATATGAGtttgtctggc                      <  2:554524/251‑9 (MQ=255)                                                                                                                                                                                                                                       gAGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATGGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACTC‑GCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTgggg         <  1:572084/251‑1 (MQ=255)                                                                                                                                                                                                                                        aGCCACATGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGt        >  2:482417/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTGCTTGGTTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAc  >  1:89305/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CCCTTTTGTTACCTCCTGGACCGCTTCTTGTCCTTCCAGTCGGAGCGAGAGTTAGCCGGGGCGGACTTTCACCGACGGAAAGATCGTGCATTTCAAGGCACACACAACAATATCTGCGAGCGTGCTATCCGTCCGGTCGTTATGGGACGAAAGGCCTGGTTGTTCGCTGGTTCACTGGTGGCCGGGAACCGCGCGGCACAGATAATGAGCCTTCTGGGAACCGCAGTCTGGAGCCACATGCTTGGCTGACGGACGTCCTGACGCGTCTGCCGGAGTGGCCGGAGGAGAGGTTGGCTGAGTTACTGCCTCTTGAGGGCTTTACCTTCTTCGGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACAC‑GCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAC  >  NC_000913/4501255‑4501742

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