Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	1,126,647	+C	intergenic (‑101/+48)	tdcG‑II ← / ← gcvP‑I	L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,126,645	1	.	C	100.0%	103.0 / NA	30	intergenic (‑99/+50)	tdcG‑II/gcvP‑I	L‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/14);  total (16/14)

GTGGCCTTTGCCTGTGGCGCCCAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGTCCTCCACCGCGCGGATT  >  AE015451/1126364‑1126938                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                        gTGGCTTTTGCCTGTGGCGCCCAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGcc                                                                                                                                                                                                                                                                                       <  1:735443/301‑1 (MQ=255)        ttGCCTGTGGCGCCCAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGCAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAA‑ccccctccgacccacagccggccacc                                                                                                                                                                                                                                                                               >  2:318783/1‑277 (MQ=255)                  cgcCCAGTGAGCCATACAGCTCGGCTGTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGGAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGt                                                                                                                                                                                                                                                                      <  2:78712/301‑1 (MQ=255)                                            ttgaccatGGCGGTGCGGGCCAGCAGGCCGTCCCGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg                                                                                                                                                                                                                                            <  1:60755/294‑1 (MQ=255)                                                  cTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGct                                                                                                                                                                                                                                      >  1:862678/1‑301 (MQ=255)                                                   gtggggTGCGGGCCAGGAGGCCGTCACGGCGCGGCCCTTCGGCAAACCGCGCCGCGGCGCGCATCGGGAGGGCCGGGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                                                                                                                                                                                     <  2:8894/297‑1 (MQ=255)                                                               ccAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTATACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGGCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAGCAGGCGTCCTGATAGCTCTCGATCGACggg                                                                                                                                                                                                                         >  1:275112/1‑301 (MQ=255)                                                                                               aaCCGCGCAGCGGCGAGCATCGGTCCGACCGTGTGGGAGCTGGAGGGGGTGAAGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCAACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac                                                                                                                                                                                         <  1:875746/301‑1 (MQ=255)                                                                                                   gcgcAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                                                                                                                                                                                     <  1:323799/301‑1 (MQ=255)                                                                                                      cAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCGGATAGCTCTCGATCGACGGGCAGGCGCACACAAGGGTGCGGTCGCCGAACACggtggcg                                                                                                                                                                                  >  1:19696/1‑295 (MQ=255)                                                                                                                 ctCGGGCCGACCGGGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc                                                                                                                                                                       <  2:356323/300‑1 (MQ=255)                                                                                                                            ggtgtGGGAGCTGGAGTGGCCGATGCCGATCTTGAACAGGTCGAAGACCCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGt                                                                                                                                                            <  1:17105/300‑1 (MQ=255)                                                                                                                                          aGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGTCGGGGAGGGCTATCAGGCGTCCTGATAGCCCTCGATCGACGGGCAGACGCAAACCGGGCTGCGGTCGCCGAACACTTTGTCGACCCGGCCGCCCGGCGTCCATCACTTGCCTTgcacc                                                                                                                                              >  2:681107/1‑301 (MQ=255)                                                                                                                                            gggCCGATGCCGATCTTGAACAGGTCGAAGGCGCTCAGTGACATGTGGGCCTCCTGCCAGGTTGTTGGGTGGACATTGTGGCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAg                                                                                                                                            <  2:19694/301‑1 (MQ=255)                                                                                                                                                                   gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGGTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCt                                                                                                                     <  2:862370/301‑1 (MQ=255)                                                                                                                                                                      gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCATGCACGCCCGGGGGTCTACCGCCTGcc                                                                                                                  >  2:480799/1‑300 (MQ=255)                                                                                                                                                                              tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCAACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGAGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCCGCCAGTCCTTGCCTTCCACCAGGCTCGCCCGGGCGTATACCGCCTGCtaccgcctg                                                                                                          >  2:574044/1‑293 (MQ=255)                                                                                                                                                                                       tgtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTAGACCCGGCCGAACGGCGGCCAGTACTTGCCTTCCACCCGGCTCGCACGGGGGTCTGCCGCCTGCTCGCGTCCGTACCCATgg                                                                                                 >  2:223137/1‑301 (MQ=255)                                                                                                                                                                                            gcctcctGCCAGGTTGTTGGGTGGACAGTGAGACCTGTAGGCGCGGGTGTAACCGCGAATGCGGAGGTATAGCAACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCa                                                                                            <  1:712242/300‑1 (MQ=255)                                                                                                                                                                                                cctaccAGGTTGGTGGTTGGACATTCAGCCCTGTTGGAGCGGGTTTACCCGCGGATGCGAAGGTATTACCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCg                                                                                        <  1:93443/297‑1 (MQ=255)                                                                                                                                                                                                                      aTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCCCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CCGCCGCAgtgg                                                                 >  2:844404/1‑300 (MQ=255)                                                                                                                                                                                                                               tGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGCTAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAAGCTCGCCACGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAATCCGCT‑CAGCCGCCgtgccccgtgcgt                                                        >  1:727463/1‑291 (MQ=255)                                                                                                                                                                                                                                                   cgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCGGCTCGCGGCTGTAGCCATTGCTCCATTCGCCAACCAGCT‑CAGCCGCCGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTcc                                   >  2:247958/1‑301 (MQ=255)                                                                                                                                                                                                                                                      aTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGTTATACCGCCCGCCCGCGGCTTTAACCATGGGCCCATTCGCCAACCAGCTACA‑CCGCAGCGTGCGGTCGCG‑CTTTCAGCGGGTTTTTGTCCttt                                >  2:446793/1‑300 (MQ=255)                                                                                                                                                                                                                                                          ggcggTACAGCCACCATCCCGTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCa                            <  1:139976/298‑1 (MQ=255)                                                                                                                                                                                                                                                           gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAg                           <  1:104520/301‑1 (MQ=255)                                                                                                                                                                                                                                                                      accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGtcct                >  1:439341/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                  gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGTCCTCCACCGCGCggg    >  1:546980/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                   cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGCGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGCAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGCCGACCGTCCTCCACCGCTCGGAt   >  1:494719/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                    gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCC‑CAGCCGCAGTGTGCGGT‑GCGTTTTTGAGGGGGTTGTTGTCCTTGTCCAGGCTCCCGTCCTCCCCCGCGCGGAtt  >  2:813936/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                    gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCC‑CAGCCGCAGGGTGCGGT‑GCGTTTTTCAGCGGGTCGTCCTCATTGTCCAGGCTACCGTCCTCCACcgctccgatt  >  1:347482/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                        GTGGCCTTTGCCTGTGGCGCCCAGTGAGCCATACAGCTCGGCTTTGACGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCCATTCGCCAACCAGCT‑CAGCCGCAGTGTGCGGT‑GCGTTTTTCAGCGGGTTGTCGTCCTTGTCCAGGCTACCGTCCTCCACCGCGCGGATT  >  AE015451/1126364‑1126938

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