Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	353,067	A→G	intergenic (+34/+84)	hisF → / ← cbcV	Imidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	353,067	0	A	G	100.0%	129.9 / NA	39	intergenic (+34/+84)	hisF/cbcV	Imidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/20);  total (19/20)

AGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTG‑CCCATTTCTGCATGTCC  >  AE015451/352794‑353358                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                     aGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCGGCTTTGCAGCCCTTCTCGGGCACGCCCGCCCCCACATGTTTGCTg                                                                                                                                                                                                                                                                           >  2:314967/1‑301 (MQ=255)      gtacacgggcGGCAGGAACAACGGCTTCGCCCTGGTCGTTACCCGGGCTATCAGTGATGCGCTTGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                                                                                                                                                                                      <  2:536091/292‑1 (MQ=255)               cGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGCTCTCtgttggg                                                                                                                                                                                                                                                             >  2:193480/1‑297 (MQ=255)               cGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgggg                                                                                                                                                                                                                                                             >  1:263182/1‑300 (MQ=255)                     gaagaaCGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTTCAGCCCTTCGCGGGCAAGCCCGACCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                                                                                                                                       >  1:521112/1‑301 (MQ=255)                        gaaCGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTg                                                                                                                                                                                                                                                    <  1:133481/301‑1 (MQ=255)                                                      tATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                                                                                                                                                                                      >  1:99368/1‑301 (MQ=255)                                                           gttatgggcTGGGGATTCCGGTTATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                                                                                                                                                                                 <  1:315073/294‑1 (MQ=255)                                                             gaggcgcttggcatgccGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGAATTCTGGAGGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGTTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaa                                                                                                                                                                                                               <  1:20636/286‑1 (MQ=255)                                                             gaggcgcTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCCGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaa                                                                                                                                                                                                               <  1:287107/298‑1 (MQ=255)                                                                 cgcTGGGGATTCCTGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGAACGTGGCAGACGGCATTCTGGAAGGGCATGCAAGGGCGGGGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCGGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCAGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                                                                                                                                           <  1:656052/301‑1 (MQ=255)                                                                                            ggggggTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACACGTTTGCTGAGATCTCTGTGGTAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATGCCTTCAACCGTGGTTCTTCCCGACCAGCGCGTGGtcaa                                                                                                                                                                                >  1:244234/1‑298 (MQ=255)                                                                                            ggggggTGGGTAACCTGAAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGTTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGGGCGGGATGGTCGTTCGCTGATCCCTTGACCGGGGCTGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGACCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                                                                                                                                <  1:55213/301‑1 (MQ=255)                                                                                                ggtcggtcaCCTGCAGCACCTGGCGGACGGCATTGTCGAAGGGCAGGCGAGCGCGGTGCTGGCGGCCAGCATCTTCCACTCTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACGTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGATGCTTTGCAGCCCTTCGCGGGCACGCCCGCCCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                                                                                                                                                                            <  1:181555/293‑1 (MQ=255)                                                                                                                           cGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTTGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGTGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTACTCCCGAGTATCCCCACAAcccc                                                                                                                                                 >  1:838941/1‑300 (MQ=255)                                                                                                                                ttCTGGAGGGGCATGCCAGCGCGGGGCTGGCGGCCAGCAGCTTCCACTTTGGCGAGTACCCAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGt                                                                                                                                            <  2:76219/301‑1 (MQ=255)                                                                                                                                     gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCCCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                                                                                                                                       >  1:406691/1‑301 (MQ=255)                                                                                                                                                   ccgggTGCTGGCGGCCAGCATCTTCCACTGTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt                                                                                                                         <  1:329281/298‑1 (MQ=255)                                                                                                                                                    gcgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttt                                                                                                                        >  1:571698/1‑300 (MQ=255)                                                                                                                                                           gggcggcGAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCgg                                                                                                                 <  2:99338/300‑1 (MQ=255)                                                                                                                                                                     cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACg                                                                                                       >  2:613156/1‑301 (MQ=255)                                                                                                                                                                         taacACTTTGACGAGTGCACGGTGCCGGAGGTCAAGGCCTAGATGGGTTACCGGGGCAGCGACGTTCGCTGATCCCTTGACCGGGACTGCTATGCAGCTCTTAGCGGGCCCGCCCGCTCCCACCGGTTTGCTGAGTTCTCTGTGGGAGCGAGCGTGCACGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATc                                                                                                   <  1:134044/298‑1 (MQ=255)                                                                                                                                                                              cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGTATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAg                                                                                              <  1:822415/301‑1 (MQ=255)                                                                                                                                                                                  ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTGGCCGAGCAGCGCGTGGTAAAGCTCGGTTTCCCCGAGTATCCCCACCACCCTGTCGTTACCCTGCCGCACCATCTTGTTTCCGGGTTGACTACGAATcctcaattcc                                                                                          >  2:587920/1‑292 (MQ=255)                                                                                                                                                                                      gagagacggTGCCGGGGGCAAAGGACTACATGGGCTGGCGAGGGATGGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCAc                                                                                      <  2:571500/293‑1 (MQ=255)                                                                                                                                                                                        tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGTTTTGATAACGCATCTGCAGTGCCTCACGc                                                                                    >  2:679695/1‑301 (MQ=255)                                                                                                                                                                                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTc                                                                         <  1:193549/301‑1 (MQ=255)                                                                                                                                                                                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCCTGTTGCCGGTTTGATAACGAATCTGCAGTGCCCCACGCATGCCGCTGTc                                                                         >  1:91196/1‑301 (MQ=255)                                                                                                                                                                                                          aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGCTCTCTGCGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTGGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACAACCCTGTCGTTAGCCTGCAGCACCAGCTCGTTTCAGGTTCGATAACCAATATGCACTGCTTCCCGCATGCCGATGTCGGCgcgt                                                                  >  2:807078/1‑298 (MQ=255)                                                                                                                                                                                                              ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTTCAGCACCAGCTTGTTGTCGGTTTGATAACGAATCTGCAGTGCCCCACGCATGCCGATGTCGGCGTCcaccc                                                              >  1:449835/1‑300 (MQ=255)                                                                                                                                                                                                                 acgagcgttcgcgcgtgagcgtagggcGCTGATCCAATGACCGGGGCTTGTCTGCAGCCCTTCGCGGGCACACCCGCGGCCACAGTAGTGCGGAGATCTCTGTGGGAGCGGGCGTGCCAGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACgg                                                           <  1:234016/275‑1 (MQ=255)                                                                                                                                                                                                                      gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGACGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCCGTTTGATAACGAATCTGCCGTTCCTCACGCATCCCGATGTCGGCGTGCACCACGGTCGGc                                                      >  2:167011/1‑301 (MQ=255)                                                                                                                                                                                                                            ggcgGTATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTgcgc                                                <  2:263092/300‑1 (MQ=255)                                                                                                                                                                                                                                        gtccGCTGATCCATTGTGCGTGGCTGGTTTGCAACGCTGCGCGGGGACGCCCGCTGCCAAAGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTcc                                    <  1:868002/298‑1 (MQ=255)                                                                                                                                                                                                                                               gATCCATTGACCGGGGCAGCTTTGCAGCCCTTCGCGGGCAGGCCCGCTACCACAGGTTTGCTGAGGTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCtt                             <  1:331629/301‑1 (MQ=255)                                                                                                                                                                                                                                                         ccGGGGCTGCTTTGCGGAACTTCGCGGGGACGCGCGCTCCCACAGGTTTGCGGAGATCTCTGTGGGAGCGGGCGTGCCCACGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTG‑c                  <  1:290587/301‑1 (MQ=255)                                                                                                                                                                                                                                                                    ttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGCTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCACCGAGTATCGCCACCACCCGGTCGTTATCCTGCAGCACAAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGACGATGTCGGCGTGCACCACGGTCGGCCTGTGCTCCAGCAGATCACCATCTTGCCCCGGTGACCCAttcagcc        >  2:89549/1‑296 (MQ=255)                                                                                                                                                                                                                                                                       cAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAACGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTGGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTG‑CCCATTTCCGCaggt    >  2:806331/1‑298 (MQ=255)                                                                                                                                                                                                                                                                          cccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGCGGGAGCGGGCGTGCCCGCGCAGCGCCGCGAAGCCGCCCCAATTGCTTCAATCGTGGCTCTTGCCGAGCAGCGCGTGGCACAGCTCGGTATCCCCGAGTACTCCCACCACCCTGTCGTCATCCTGCAGCAACAGCCTCTTTGCGGTTTGATTACTAATCTGCAGCGCCTCACCCATGCGGCTGCCGGCCTCCACCCTGGTCTGCCTCCGCTCCAGCCGCTCCACTTCGTGccccagagccctattcttcatgcccc  >  2:558038/1‑279 (MQ=255)                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                     AGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTG‑CCCATTTCTGCATGTCC  >  AE015451/352794‑353358

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