Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	2,935,792	G→A	intergenic (+87/‑175)	PP_2568 → / → PP_2569	Ring‑cleaving dioxygenase/Metabolite MFS transporter, MHS family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	2,935,792	0	G	A	100.0%	149.3 / NA	43	intergenic (+87/‑175)	PP_2568/PP_2569	Ring‑cleaving dioxygenase/Metabolite MFS transporter, MHS family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/22);  total (21/22)

GGGCGTGGAAGGGCGCAACCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGT  >  AE015451/2935495‑2936072                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                         gggggTGGAAGGGCGCAACCTGAACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCGATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGAccc                                                                                                                                                                                                                                                                                       <  1:626361/297‑1 (MQ=255)                  cccTGCACCATTTCTGCCTGCGCATAGGGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAg                                                                                                                                                                                                                                                                      <  1:251722/300‑1 (MQ=255)                        cccATTCCTGCCTGCGCATCGAGCCGTTGGACGAGGCCGCGGTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATc                                                                                                                                                                                                                                                                <  2:96054/300‑1 (MQ=255)                                 ccctgcgctacgggccggtcgaccgggCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTGCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATg                                                                                                                                                                                                                                                       <  2:474997/276‑1 (MQ=255)                                    ccttccgatcacccgtggacgaGGCCGCGTTGACGTCCTATCTGGAGGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCtt                                                                                                                                                                                                                                                    <  1:567146/284‑1 (MQ=255)                                       gcATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCAt                                                                                                                                                                                                                                                 >  3:151654/1‑301 (MQ=255)                                        ggctcttccgatctacgaGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGATGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATg                                                                                                                                                                                                                                                <  4:411801/287‑1 (MQ=255)                                               ccGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCttagatcggaaga                                                                                                                                                                                                                                         >  3:567148/1‑290 (MQ=255)                                                     gacgaGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGCTTATGCTTCATTAACCTttctgacc                                                                                                                                                                                                                                   >  2:549270/1‑298 (MQ=255)                                                      acgaGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGATGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGTTGGATATCCAGCCAATCCCGATTAttccccatgaagtccgaaaagtg                                                                                                                                                                                                                                  >  2:411692/1‑279 (MQ=255)                                                                          tatctaggaggccccgGGGGAAAGGTGGTGCCTGCGTAGAAGCGTTATGGGGCCGCAGGGGGGGGGCTCTCGCTGTACTGCTCCGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAACCCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                                                                                              <  2:595127/288‑1 (MQ=255)                                                                             cTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGCTATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCCATCCGCTGCGGCCCGACTCCTCGTg                                                                                                                                                                                                           >  2:749408/1‑301 (MQ=255)                                                                             cTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTg                                                                                                                                                                                                           >  3:92794/1‑301 (MQ=255)                                                                               ggAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGCTTATGCTTCATGAACCTTTCTTCCATCCGAGGCCGCCCCCCtcggccttcc                                                                                                                                                                                                         >  2:3440/1‑294 (MQ=255)                                                                                gAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCtt                                                                                                                                                                                                        <  1:151652/301‑1 (MQ=255)                                                                                 aaGCCGCCGGGGTAAAGGTGGAGCCTGCCGTGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTc                                                                                                                                                                                                       <  1:92791/301‑1 (MQ=255)                                                                                          gggTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAAccc                                                                                                                                                                                              >  3:459904/1‑300 (MQ=255)                                                                                                      ttttcccatcacgacgcctttccgtactcggggagggcCTTTCGCGGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCa                                                                                                                                                                                  <  2:33549/272‑1 (MQ=255)                                                                                                            ccGAGAAGCGTTATGGGGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCt                                                                                                                                                                            <  4:3440/301‑1 (MQ=255)                                                                                                                gaAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTgggg                                                                                                                                                                        >  4:512710/1‑300 (MQ=255)                                                                                                                gaAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGt                                                                                                                                                                        <  4:549427/301‑1 (MQ=255)                                                                                                                 aaGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCCCGCAACCAAACCAAGCCCCCGGGcccggtgg                                                                                                                                                                       >  2:31242/1‑295 (MQ=255)                                                                                                                          ggTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTgg                                                                                                                                                              >  1:563879/1‑300 (MQ=255)                                                                                                                           gTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCCATCCGAGGCGGTCCGCCTCGTCGTGCTTCTCGCAACCCAAGCCACGGCCAGGGCCTGGGTGCGACGGGGaa                                                                                                                                                             >  1:560129/1‑301 (MQ=255)                                                                                                                                gAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAACGCAATCGCCCGGGCCTGGGTGCGACggggagcggct                                                                                                                                                        >  2:70835/1‑292 (MQ=255)                                                                                                                                  aggggaggggCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCCCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCacgatcggcagagcccccgtctgaactcc                                                                                                                                                      >  4:33560/1‑273 (MQ=255)                                                                                                                                       gggggCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAgg                                                                                                                                                 <  1:459902/300‑1 (MQ=255)                                                                                                                                                 gCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGcc                                                                                                                                       <  3:563881/301‑1 (MQ=255)                                                                                                                                                    gTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGt                                                                                                                                    <  4:31252/301‑1 (MQ=255)                                                                                                                                                                                    aGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGt                                                                                                    <  4:70854/301‑1 (MQ=255)                                                                                                                                                                                              cccggtgtgcccTGGCGGGGCAAACCGCCCCGTCCAAAACCGTCGGCCCACCCGCAATCGCCCGGCACCCCGGGGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGc                                                                                          <  2:512560/291‑1 (MQ=255)                                                                                                                                                                                                             ctgcttcccacGCCCCCGCCAACATGGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAgc                                                                           <  2:450723/291‑1 (MQ=255)                                                                                                                                                                                                                     aCCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCa                                                                   >  4:597405/1‑301 (MQ=255)                                                                                                                                                                                                                       cctacacgatccccttccgttcgaccccGCTATCCCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAAc                                                                 <  2:745905/276‑1 (MQ=255)                                                                                                                                                                                                                          ccccGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAgcagatcggaagagc                                                              >  4:450852/1‑288 (MQ=255)                                                                                                                                                                                                                               ccAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATc                                                         >  1:494056/1‑301 (MQ=255)                                                                                                                                                                                                                                       gTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGGCCATCCACGGGCAGCCCGCCCCAACGCGGTATcccccaggg                                                 >  2:28735/1‑293 (MQ=255)                                                                                                                                                                                                                                              aaCCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAAcagatcggaagagcacacgtctga                                          >  4:746184/1‑278 (MQ=255)                                                                                                                                                                                                                                                  cGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAgg                                      <  4:749687/301‑1 (MQ=255)                                                                                                                                                                                                                                                    cTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCCAGGGCAGCCCGCCCCAACCCGGTATCACCAaggaggcacgccaggtg                                    >  2:781491/1‑290 (MQ=255)                                                                                                                                                                                                                                                                                 cgtcCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTAtt       <  4:28745/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                  ggccAGGCCGCCCCCCCGGGCAGGACCCGGAGGGTGGATATCCTGCCAACCCCGATTATGCTCCATGAACCTGGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTc      <  2:597208/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                      aGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAgg  >  3:299509/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                         GGGCGTGGAAGGGCGCAACCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGT  >  AE015451/2935495‑2936072

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: ‑