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%	128.8 / NA	38	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 (17/21);  total (17/21)

ATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGA  >  AE015451/2935521‑2936070                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                        aTTTCTGCCTGCGCATGGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGGGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCg                                                                                                                                                                                                                                                            <  1:741761/301‑1 (MQ=255)      ggcctgcGCATCGGGCCGTGCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTAt                                                                                                                                                                                                                                                       <  1:316601/300‑1 (MQ=255)       cagacgtgtgctcttccgatctcgaGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATg                                                                                                                                                                                                                                                      <  1:743501/279‑1 (MQ=255)             gcATCGAGCCGTTCGACGAGGCCGCGTTTACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCAt                                                                                                                                                                                                                                                >  1:673517/1‑301 (MQ=255)                gcgagccgtgcgagcgggCCGGGTTTACGGCCTATCTGGAAGCCGCGGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTGTGGTGCCGAAGGGGAGGGGCTTTCGCTGAACTGCTACGGCCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGaa                                                                                                                                                                                                                                             <  2:325774/285‑1 (MQ=255)                      ggTTCGACGAGGCCGCGTTGACGGCCTAGCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTtgc                                                                                                                                                                                                                                       <  1:348650/300‑1 (MQ=255)                        gtctagcgagcagcgTTGACGGGCTACATGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGGGAA‑ACGTTATGGTGCCGAAGGGGAGGGGCGTTCGCTGTACTGCTACGACCTGGATGGTAATCAGATCGAGTGGAAGGGGCCGGTAGGGCCTGGCGGGTGAAACCGCCCCCGCCAGAATCGTCAGCCAGCCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTtgctg                                                                                                                                                                                                                                     <  2:409747/289‑1 (MQ=255)                             cgaGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTTCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATgagatcggaagagcggcgggtag                                                                                                                                                                                                                                >  2:743187/1‑279 (MQ=255)                               gcggcgcggcggcGGCCTATATGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCAGAGGGGAGGGGCTTGCGCTGTACTGCTACGACCCGGATGGTAATCAGGTGGAGTTGAAGGGGCGGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCg                                                                                                                                                                                                                              <  2:229514/289‑1 (MQ=255)                                             gCCTATCTGGTAGCCGCCGGGGTAAAGGTGGGGCCTAGCGAGGA‑GCGTGAGTGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATTGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGGGAACCCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtc                                                                                                                                                                                                                <  2:60256/301‑1 (MQ=255)                                              ccTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcg                                                                                                                                                                                                               <  1:699609/301‑1 (MQ=255)                                                ggagttcagacgtgtgctcttccgatctaGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                                                                                             <  1:137962/273‑1 (MQ=255)                                                                 ggTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCCATCCGAGGCGGCCCGCCTCGTCGGGCTTCGCTCAAACaa                                                                                                                                                                                            >  1:482356/1‑301 (MQ=255)                                                                   taatggTGGAGCCTGACG‑GAACGCGTTATGGTGCCGAAGGGGAGGAGCTTTAGCTGTACTGCTCCGACCCGGATGGTGATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaag                                                                                                                                                                                          <  2:515757/297‑1 (MQ=255)                                                                            aGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCccgtcagatcggaaaggcgtcggggagggaaag                                                                                                                                                                                 >  2:137903/1‑273 (MQ=255)                                                                                  ccGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTTATGAACCTTGCTGCGATCCGAGGCGTCCCGCCTCGTCGTGATTGCCGCAACCAAAGCAAACGCCacgcccc                                                                                                                                                                           >  2:743291/1‑295 (MQ=255)                                                                                               tgggTGCCGAAGGGGGGGGGCTTTCGCTGTACGCCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGt                                                                                                                                                               <  2:136675/299‑1 (MQ=255)                                                                                                               ggggCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCTAAGCGTCCCGCCTCGTCGTGCTTCCCGCAACCCAAGCAACCGCCAGGGCCTGGGTGCGCCGGTGAAGT‑GGTTGTGAggg                                                                                                                                              >  2:68656/1‑300 (MQ=255)                                                                                                                          tGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATTCGAGGCGGCCCGCCTCGTCGTCCTTCGCCCAACCAAAGCAAGCGCCAGGCCCTGGGTGCGACGGTGACGT‑GGTTGTGAGGTGCTCCGGCCTg                                                                                                                                   >  2:111961/1‑301 (MQ=255)                                                                                                                               tGCTACGACCCGGATGGGAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTTGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGa                                                                                                                              <  1:388029/301‑1 (MQ=255)                                                                                                                                     gACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGATGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACAt                                                                                                                        >  1:678592/1‑301 (MQ=255)                                                                                                                                           gATGGTAATCAGGTCGAGTTGAAGAGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAAACAAAGCAAGCGCCAGGGCCTGGGTGCGTCGGTGAAGT‑GGGTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGa                                                                                                                  >  1:387196/1‑301 (MQ=255)                                                                                                                                                tAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCCGCGGTTGAGT‑GGTTGTGAGGTTCGCCGGCCTGGTGGACAACATGACGGAAGCaa                                                                                                             >  2:745017/1‑301 (MQ=255)                                                                                                                                                 aaTCAGTTCGAGTTGAGGGGGCCGGTGGGGCCTGGCGGGTGAACCCCCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAc                                                                                                            <  2:295011/301‑1 (MQ=255)                                                                                                                                                    cAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGCGGCGGCCCGCCTCGTCGTGCTTCGCCCAACCAACGCAAGCGCCAGGGCCTGGGTGCGACGGTGACGTGGGTTG‑GAGGTGCTCGGGCGCGTCTGACAACCTGACgcggggaccccc                                                                                                         >  2:327267/1‑290 (MQ=255)                                                                                                                                                     gggTCGAGTTGAGGGGGCCGGGAGTGCCTAGCTGGTGCAACCGCCACCGGCAAAATCGTCAGCCAACCCGCTAGCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATAGCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGCTCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCAt                                                                                                        <  2:240442/300‑1 (MQ=255)                                                                                                                                                                              gCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACg                                                                               >  1:278661/1‑301 (MQ=255)                                                                                                                                                                              gCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACg                                                                               <  1:474372/301‑1 (MQ=255)                                                                                                                                                                               ccTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGGGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACgg                                                                              <  1:327394/301‑1 (MQ=255)                                                                                                                                                                                cTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCaacgg                                                                             >  2:135238/1‑297 (MQ=255)                                                                                                                                                                                              cgcacccGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCTCCCTCTCGGGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAAc                                                               <  1:735654/297‑1 (MQ=255)                                                                                                                                                                                                                  gCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGgaggag                                           >  1:430340/1‑301 (MQ=255)                                                                                                                                                                                                                        ccGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCCGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAg                                     >  1:257279/1‑301 (MQ=255)                                                                                                                                                                                                                                       aCCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGccccc                      >  2:332705/1‑298 (MQ=255)                                                                                                                                                                                                                                         gcgggTGCTGCGTGTCGGCCCGTCCTCACGCTCGTGCAGGACCGTGCGGGGTGAAATCCAGCCAATCCCGATTATGCGTCATGAACCTGTCTGCGATCCGCGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTc                    <  2:730221/298‑1 (MQ=255)                                                                                                                                                                                                                                               gcttcgtttcgtcCAGTCCGCCCTCTCGTGCAGGACCCATGGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCTCCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTggg              <  2:734349/293‑1 (MQ=255)                                                                                                                                                                                                                                                  gcGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCAc           <  1:135297/301‑1 (MQ=255)                                                                                                                                                                                                                                                           ccAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCgg  >  1:754242/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                        ATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAA‑GCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGT‑GGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGA  >  AE015451/2935521‑2936070

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