Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,740,811	(T)7→6	intergenic (+48/+705)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,740,805	0	T	.	100.0%	197.5 / NA	43	intergenic (+42/+711)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (22/21);  total (22/21)

ACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTT  >  AE015451/4740531‑4741099                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                       ccTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                                                                                                                                                                             <  1:403213/300‑1 (MQ=255)  cctCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCATGTGCCGGCCTCtt                                                                                                                                                                                                                                                                            <  3:55763/299‑1 (MQ=255)                cgctGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCTTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGTCCTTGcccc                                                                                                                                                                                                                                                              >  2:706309/1‑300 (MQ=255)                 gctGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                                                                                                                                                             >  4:701767/1‑301 (MQ=255)                    gTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                                                                                                                                          >  1:289979/1‑301 (MQ=255)                     tgctcttccgatctcTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                                                                                                                                                                         <  4:423987/287‑1 (MQ=255)                       cccGACCGCCGACTTCTACGC‑GGGAACATTCCTCAAGGCGACCGGCATCCCGACCCGCATGTTCACCGTGCTCTTCCCCCTGGCACGTACCGTGGGCTGGATCTCGAACTGAAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                                                                                                                                       <  1:214793/301‑1 (MQ=255)                       cccGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGCTCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCCGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGccctcgcagcg                                                                                                                                                                                                                                                       >  2:556560/1‑293 (MQ=255)                       cccGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGGGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCTGCCTCTTCGCCGCCTTcccccccaaagg                                                                                                                                                                                                                                                       >  2:80856/1‑293 (MQ=255)                                   cTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGCCTTGCCCGCGAagagatcggaagagcg                                                                                                                                                                                                                                           >  2:423878/1‑289 (MQ=255)                                                       cacAAGGCGATCGGCAGTCCGACAAGCATGTTCACCGTGATCTGCGCCCCGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtttt                                                                                                                                                                                                                        <  4:556730/299‑1 (MQ=255)                                                                   ggCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTGATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCGCAAAGGCCCGCCCAGTCCCCACACAACCCCAGTTTTtccccccacccc                                                                                                                                                                                                            >  2:468386/1‑291 (MQ=255)                                                                          cGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGTCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTCTCTCCGGCCCCCCACGCAg                                                                                                                                                                                                     >  2:677089/1‑301 (MQ=255)                                                                            cccAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                                                                                                                                                                                   <  3:289978/300‑1 (MQ=255)                                                                             ccAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGGGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                                                                                                                                                                                  <  4:80875/301‑1 (MQ=255)                                                                             ccAGCATGTTCACCGTGATCTTCGCCCTGGCACCTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGTGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                                                                                                                                                                                  <  4:706554/301‑1 (MQ=255)                                                                                      gacgccgttccgatgtccTGGCACGTACCGGGGGCTGGATCTCGAACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGt                                                                                                                                                                                         <  2:605685/285‑1 (MQ=255)                                                                                                cTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCAAAGCCACCACAAAACCTCAGTTTTTCTCCGCCCGCTAACGCAGCCCTTTCAGGGTATTGAACGgc                                                                                                                                                                               >  3:485874/1‑301 (MQ=255)                                                                                                      ccTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTagatcggaagagcaca                                                                                                                                                                         >  4:605882/1‑286 (MQ=255)                                                                                                       tcggCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                                                                                                                                        <  4:651469/299‑1 (MQ=255)                                                                                                          gCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGACCGGCCCATcgcaga                                                                                                                                                                     >  2:651258/1‑296 (MQ=255)                                                                                                                         gggATCTCGCACTGGAAAGAAATGCTCTCCCGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAAAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca                                                                                                                                                      <  4:468515/300‑1 (MQ=255)                                                                                                                          ggaacacGCACGGGAAAGAAATGCACTCCCGCCCGTACAAGATCGACCGTCCG‑CGCCAGCTGTAAACCGGCGAACAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAACTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccac                                                                                                                                                     <  1:485872/295‑1 (MQ=255)                                                                                                                             tctcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCAt                                                                                                                                                  <  4:677311/301‑1 (MQ=255)                                                                                                                                  cagtggaaagaaatctccccccgCCCGTACCACATCGGCCGTCCGACG‑CAGCTGTACACCGGCGAGCAAGAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCg                                                                                                                                             <  2:701522/279‑1 (MQ=255)                                                                                                                                    cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCGGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGATTTTCTCCGCCAGCGCACGCAGCCCTTCCAGGGTATAGAACGTCGAATACACCACAAACTTGTGGGGCACCATCGCCggg                                                                                                                                           >  2:43347/1‑300 (MQ=255)                                                                                                                                         aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCt                                                                                                                                      <  4:43364/301‑1 (MQ=255)                                                                                                                                                               cAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAcc                                                                                                                >  4:230518/1‑301 (MQ=255)                                                                                                                                                                      ggCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAgg                                                                                                         >  3:127148/1‑301 (MQ=255)                                                                                                                                                                                        gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCa                                                                                        <  1:127146/301‑1 (MQ=255)                                                                                                                                                                                                                   cacgacgctcttccgctctaaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCg                                                             <  2:26278/282‑1 (MQ=255)                                                                                                                                                                                                                            gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACtttt                                                    >  4:438773/1‑301 (MQ=255)                                                                                                                                                                                                                             aaataCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTg                                                   <  2:230464/297‑1 (MQ=255)                                                                                                                                                                                                                               aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTTAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCCCCCTTGGGCACCAGCTTCCAGAAGCCGTCGCCGTTTGcc                                                 >  1:367061/1‑301 (MQ=255)                                                                                                                                                                                                                                      aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTTGGCCCCAGCTTCCAGAAGCCgagatcggaagagcacacgt                                          >  4:26288/1‑282 (MQ=255)                                                                                                                                                                                                                                        ttccgatctagatgCCGAAGGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCAGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACa                                        <  2:189321/288‑1 (MQ=255)                                                                                                                                                                                                                                                 aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCAACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGACCaagatcggaa                               >  4:189367/1‑293 (MQ=255)                                                                                                                                                                                                                                                  ggagttcagacgtgtgctcttccgatctggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTtct                              <  4:720769/273‑1 (MQ=255)                                                                                                                                                                                                                                                     gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCtt                           >  4:604396/1‑301 (MQ=255)                                                                                                                                                                                                                                                      ccgatctcAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTc                          <  4:318582/294‑1 (MQ=255)                                                                                                                                                                                                                                                      ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTACGTTCTTTCACCCTTGGGCACCAGCTTCCAGACGCCGTCGCCTTTTGCCCCCCGATCCAAAACCTTCTATTc                          >  2:561859/1‑301 (MQ=255)                                                                                                                                                                                                                                                             cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCCCCCGAACAAACCCCTTCTCTTcccaccgc                   >  2:318509/1‑295 (MQ=255)                                                                                                                                                                                                                                                                              ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTtctagatcggaagagcgtggtgtcggggagg  >  2:720524/1‑273 (MQ=255)                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                       ACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCA‑TCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCG‑CGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTT  >  AE015451/4740531‑4741099

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