Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	1,499,480	2 bp→AC	intergenic (+57/‑106)	trpS → / → zapE	Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,499,480	0	C	A	100.0%	111.3 / NA	32	intergenic (+57/‑107)	trpS/zapE	Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/14);  total (18/14)
*	AE015451	1,499,481	0	A	C	100.0%	118.3 / NA	32	intergenic (+58/‑106)	trpS/zapE	Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/14);  total (18/14)

CCGCTGATGGCGAGCAACTGCTGTTGTCGCGCAGCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCC‑GGTGTT‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCA  >  AE015451/1499188‑1499764                                                                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                                                               gggCTGATGGCGGGCAACTGCTGTTGTCGCGCAGCTTCGCGGACGGGAGGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                                                                                                                                                                         <  2:58804/299‑1 (MQ=255)           ggtctaagtgctgTTGTCGCGCAGCTTGGCCGACGGCAAAGGTGCGGGCGCCGTGGGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGttc                                                                                                                                                                                                                                                                              <  2:11876/293‑2 (MQ=255)              gCAACTGCTGTTGTCGCGCAGCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGttcgcg                                                                                                                                                                                                                                                                           >  3:67977/1‑297 (MQ=255)                aaCTGCTGTTGTCGCGCAGCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGttcagatg                                                                                                                                                                                                                                                                         >  4:11876/1‑295 (MQ=255)                            cacgccgctcttccgatctcAAGGGTGCTGGCGCCGTGAGCCAGCAGCTGCAGCAAGGTGACGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGATCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTcc                                                                                                                                                                                                                                                             <  1:436634/282‑1 (MQ=255)                                          gACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGc                                                                                                                                                                                                                                               <  4:204674/301‑1 (MQ=255)                                               cAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTccagatcggaagagcacacgg                                                                                                                                                                                                                                          >  3:436636/1‑282 (MQ=255)                                                aaGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTaa                                                                                                                                                                                                                                         >  3:111992/1‑301 (MQ=255)                                                           ggagttcagacgtgtgctcttccgatctGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                                                                                                                                              <  4:121215/274‑1 (MQ=255)                                                                 gAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                                                                                                                                        >  3:226597/1‑301 (MQ=255)                                                                                       ggCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttgg                                                                                                                                                                                                  >  3:62135/1‑300 (MQ=255)                                                                                       ggCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGGTGTC‑GCGGGCACGCCCGCTTCCACCGACTGATTGCCATTACGCGGGCCTaacagatcagaagagagtcgtagaagggaag                                                                                                                                                                                                  >  2:121190/1‑271 (MQ=255)                                                                                          gAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                                                                                                                                                                               <  1:62133/301‑1 (MQ=255)                                                                                                   gTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACCGTGACGGCCCGGTTTTTTTTTCCGCGCTAACGaa                                                                                                                                                                                      >  2:459341/1‑301 (MQ=255)                                                                                                              aGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                                                                                                                                           <  4:459470/301‑1 (MQ=255)                                                                                                                                                   cccggtccagggCCACAGTGCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGACGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct                                                                                                                                      <  2:531411/292‑1 (MQ=255)                                                                                                                                                       ccgatcggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGTCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt                                                                                                                                  <  1:415192/295‑1 (MQ=255)                                                                                                                                                           cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGAAGATCTGAAACGTCCGGActagatcgga                                                                                                                              >  4:531561/1‑293 (MQ=255)                                                                                                                                                              gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGTCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtcttagatcgg                                                                                                                           >  3:415194/1‑294 (MQ=255)                                                                                                                                                                                       cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTTTTTTGCCCCGCTAACGAATCATGACTCCCTTTAAACGCTATCAATCAGATCTGAAACGTCCCGACTTCTTCCATTACGCGGCGCAAGAAACTGACGTGCGTc                                                                                                  >  1:114820/1‑301 (MQ=255)                                                                                                                                                                                           aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAAGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAACCGTCCCGACTTCTTCCATGAAGCGGCACAGGAACCTGCGGTGCGTccccc                                                                                              >  1:173481/1‑297 (MQ=255)                                                                                                                                                                                            ggCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTc                                                                                             <  1:226599/300‑1 (MQ=255)                                                                                                                                                                                              gggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCa                                                                                           <  1:111989/300‑1 (MQ=255)                                                                                                                                                                                                                   aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGt                                                                      <  3:114823/301‑1 (MQ=255)                                                                                                                                                                                                                                      ggACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAAGCGGGGCCGGAAACTGAGGTGCGTCACCTACCGCGCCTGTCCGCAGACCTGGGGCACGCGCAGaacaaaaca                                                   >  1:365866/1‑297 (MQ=255)                                                                                                                                                                                                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAACCTGCGGTGCGTCCCCCCCAGCGCCTGTAAGACGACCTGGCGCACGCGCCGCACAACAAGCCGGGccag                                          >  2:694019/1‑298 (MQ=255)                                                                                                                                                                                                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGGGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGGTCgg                                     >  4:57049/1‑301 (MQ=255)                                                                                                                                                                                                                                                             gAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGtt                            <  3:173483/301‑1 (MQ=255)                                                                                                                                                                                                                                                                cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCAGGTTCGGCAAGCTGTTCgg                         >  4:448945/1‑301 (MQ=255)                                                                                                                                                                                                                                                                              aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCACGAAGGACCAGa           >  4:166176/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                       ccgcGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGCTTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGCTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGAACTGGTGCAAGCGCAGAACAACAAGCCTGGCAAGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGCCTCCGGtaa  >  1:163419/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                       ccgcGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCa  <  3:365867/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                                                               CCGCTGATGGCGAGCAACTGCTGTTGTCGCGCAGCTTCGCCGACGGCAAGAGTGCTGGCGCCGTGAGCAAGCAGCTGCAGCAAGGTGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCC‑GGTGTT‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACAAGCCGGGCATGTTCGGCAAGCTGTTCGGCAAGAAGGAGCAGACTCCGGTCA  >  AE015451/1499188‑1499764

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