Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	6,013,885	T→C	intergenic (+90/+20)	xpt → / ← PP_5266	Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	6,013,885	0	T	C	100.0%	146.5 / NA	40	intergenic (+90/+20)	xpt/PP_5266	Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/20);  total (20/20)

GATGACTTCCTGGCCAACGGCAAAGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCCGCGCATGCTGTGC  >  AE015451/6013592‑6014168                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                              gttgcattaccggcgaaaggCAAAGCCTCCACAGCGGTGAGCTCGATCATCAAGCAGCAAGGTGCAACCGTGGCGGGCCTGGGTATAGTGATCGAGAAGTCGGTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATg                                                                                                                                                                                                                                                                                        <  2:14642/283‑1 (MQ=255)   tGACTTACTGGCCAACGGCAAAGCCTCCCAAGCGCGGATCTCGATCATCAAGCAGGCAGGTGCCACCGGGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGGCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCGTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGc                                                                                                                                                                                                                                                                                      <  1:10115/301‑1 (MQ=255)         gagttcagacgtgtgctcttccgatctgcTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAAc                                                                                                                                                                                                                                                                                <  3:605326/274‑1 (MQ=255)                      aaaGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGGAACCCTTTTATCTggg                                                                                                                                                                                                                                                                   >  1:287988/1‑300 (MQ=255)                                    gcTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAAcagatcggaagagcgtcgtgtagggaaa                                                                                                                                                                                                                                                     >  1:605324/1‑274 (MQ=255)                                                   ctggaggtcagacgtgtgcttttccgattcgggTAACGTGATCGAGAAGTCGTTCCAGGGCGGCCGGGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCt                                                                                                                                                                                                                                      <  4:33836/271‑1 (MQ=255)                                                    aaGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                                                                                                                                                                                     <  3:287987/301‑1 (MQ=255)                                                        ccgatcttGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTCACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaa                                                                                                                                                                                                                                 <  4:49388/294‑1 (MQ=255)                                                               tGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTCACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGTAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGTTTGCCGCCAGCCCCGCCACCACCAAGCGCTGGTaaacatcgg                                                                                                                                                                                                                          >  2:49371/1‑295 (MQ=255)                                                               tGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGAGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttctt                                                                                                                                                                                                                          >  4:771346/1‑301 (MQ=255)                                                                  cACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTcttctt                                                                                                                                                                                                                       >  3:512350/1‑301 (MQ=255)                                                                   cccGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                                                                                                                                                                                                                      <  4:174345/300‑1 (MQ=255)                                                                          cGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGcccc                                                                                                                                                                                                               >  4:753101/1‑301 (MQ=255)                                                                                               agaAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATc                                                                                                                                                                                        >  4:301076/1‑301 (MQ=255)                                                                                                 aaGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGGTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCg                                                                                                                                                                                      >  3:81413/1‑301 (MQ=255)                                                                                                                ggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTc                                                                                                                                                                       <  1:512348/301‑1 (MQ=255)                                                                                                                             tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGc                                                                                                                                                          >  4:449548/1‑301 (MQ=255)                                                                                                                                   gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCCGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCGTCCA‑‑TGCATCCGCCCCAGGTCGCTTTCATACCTCTCCGGCCCCGGt                                                                                                                                                    >  1:419003/1‑301 (MQ=255)                                                                                                                                             ccGCGTTGAATCGCTGGGCCGAGTGACTTCGCTGGAAGGTTGTGTGTTCAGTGTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACg                                                                                                                                          <  1:81411/301‑1 (MQ=255)                                                                                                                                                                    tatcttCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCAcc                                                                                                                   <  1:30776/296‑1 (MQ=255)                                                                                                                                                                         tCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCAccagatc                                                                                                              >  3:30776/1‑296 (MQ=255)                                                                                                                                                                         ccGCTGGCAGGTGGTGTGATCACCTTCATCGAAGACAAGCGAGCCCTGATTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCt                                                                                                              <  2:449419/300‑1 (MQ=255)                                                                                                                                                                                    gggggcgggaaGCTTCATCGAATAAACGCGACACCTGCTGCGCGGGAACGCCCGCTCCAACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACt                                                                                                   <  2:771066/291‑1 (MQ=255)                                                                                                                                                                                      ggttggtCAGATTCATCGAATAACCGCGACACCGGCTTCGCGGGCACGCCCGTTCCGACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTGTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTg                                                                                                 <  2:752822/298‑1 (MQ=255)                                                                                                                                                                                             cAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGTGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCa                                                                                          <  3:419005/301‑1 (MQ=255)                                                                                                                                                                                                            cctacacgaccctctcccgacctaCGCCCGCCCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGccc                                                                           <  2:93267/278‑1 (MQ=255)                                                                                                                                                                                                                               gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCAGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGctcctc                                                        >  1:235359/1‑301 (MQ=255)                                                                                                                                                                                                                                ggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCg                                                       <  4:591679/301‑1 (MQ=255)                                                                                                                                                                                                                                   aCGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGcccagatcggaagagcacacgtctga                                                    >  4:93288/1‑278 (MQ=255)                                                                                                                                                                                                                                     gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCgc                                                  >  2:194233/1‑300 (MQ=255)                                                                                                                                                                                                                                      cccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTa                                                 >  1:423832/1‑301 (MQ=255)                                                                                                                                                                                                                                        cGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCAGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTaaa                                               <  3:235357/301‑1 (MQ=255)                                                                                                                                                                                                                                                             cgTCTCGAACGCTGTGCCGTACCTGTGTAAGCTTGCTTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTc                          <  3:145642/300‑1 (MQ=255)                                                                                                                                                                                                                                                               ccccgacgctattccggttgcggtgGGAGCGGGCGTGCCCGCGAATTGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCg                        <  2:253515/279‑1 (MQ=255)                                                                                                                                                                                                                                                                 ccGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCggg                      <  3:423834/300‑1 (MQ=255)                                                                                                                                                                                                                                                                        cTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACTCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCCCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGGAAAGTCAGTGCCCATCGCACACTCCGGGAACAGcc               >  2:44122/1‑301 (MQ=255)                                                                                                                                                                                                                                                                         tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCcg              >  4:266688/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                tACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCCGCGCATGc       >  2:287914/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                   cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCgagatcggaagagcacacgtc    >  4:253573/1‑281 (MQ=255)                                                                                                                                                                                                                                                                                     gtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATTCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCCGCGCATGCTGTGc  >  4:113222/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                              GATGACTTCCTGGCCAACGGCAAAGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCC‑‑CCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCCGCGCATGCTGTGC  >  AE015451/6013592‑6014168

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