Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,553,678	T→C	C152R (TGC→CGC)	hemF →	coproporphyrinogen III oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,553,678	0	T	C	100.0%	123.0 / NA	36	C152R (TGC→CGC)	hemF	coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/17);  total (19/17)

GATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATG‑GCCTACGGCGAGCGCGAGCGCAATTTC  >  NC_000913/2553431‑2553920                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     gATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGcgcc                                                                                                                                                                                                                                                   >  2:149485/1‑251 (MQ=255)     ccTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGcgcctgcc                                                                                                                                                                                                                                               >  1:130002/1‑251 (MQ=255)                     cTCATAGCCAGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGCAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGtt                                                                                                                                                                                                                               >  1:299437/1‑251 (MQ=255)                          cGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATcc                                                                                                                                                                                                                          <  1:229600/251‑1 (MQ=255)                           gCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATccc                                                                                                                                                                                                                         <  1:14304/251‑1 (MQ=255)                            cccGGAATGTGCCGTGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCCCCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGTTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATcccc                                                                                                                                                                                                                        <  2:87295/251‑1 (MQ=255)                              cGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGt                                                                                                                                                                                                                      >  2:113332/1‑251 (MQ=255)                                 aaCTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTAc                                                                                                                                                                                                                   <  1:190101/251‑1 (MQ=255)                                                    cGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAAt                                                                                                                                                                                                <  1:301969/251‑1 (MQ=255)                                                          gATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCt                                                                                                                                                                                          >  2:281665/1‑251 (MQ=255)                                                               gCGTTTCACTGGTAGTTCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCGTCGACTTAACCCCATTCTATGGTGTTTAAGAAGATGCTATTCACTGGCATCTCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGGCGTTTATCTCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTc                                                                                                                                                                                     <  1:53297/251‑1 (MQ=255)                                                                    tCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCTCCCGGGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAAc                                                                                                                                                                                 <  2:187360/250‑1 (MQ=255)                                                                     cACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACAt                                                                                                                                                                               <  2:165928/251‑1 (MQ=255)                                                                       cTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCg                                                                                                                                                                             >  1:55842/1‑251 (MQ=255)                                                                                        tAACCCGTCTGTTCCCACCAGCCACGCGAATGCGCGGATCTTTATTGCCGAAAAACCGGATGCCGATCCCGTCTGGTGGCTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGt                                                                                                                                                             >  2:27732/1‑250 (MQ=255)                                                                                          aCCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTAtt                                                                                                                                                          <  1:81431/251‑1 (MQ=255)                                                                                            ccGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTgg                                                                                                                                                        <  1:316529/251‑1 (MQ=255)                                                                                                aTGTTCCTACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGCGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCggg                                                                                                                                                    >  2:56193/1‑251 (MQ=255)                                                                                                        accaGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACCTCGCAACGAACAGCGCGGGATTGGCGGGCTGTTCtt                                                                                                                                            >  2:45662/1‑251 (MQ=255)                                                                                                         ccaGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCttt                                                                                                                                           <  2:63489/251‑1 (MQ=255)                                                                                                            accACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGGTGTTCtttttt                                                                                                                                        >  2:212001/2‑248 (MQ=255)                                                                                                                               ttttATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTTATGACCTGAAAACGccccatt                                                                                                                     >  2:97547/1‑246 (MQ=255)                                                                                                                                 ttATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGAACTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGCGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGTCTGTTCTTTGATGACCTGCACACGCCAGAtttt                                                                                                                   >  1:209680/1‑250 (MQ=255)                                                                                                                                              cccGGGTGCCGATCCCGTCTTGTGGTTTGTCGGTGGCTTCTACTTACCCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGCCCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATTCTTCTACCTCAAACATCGCACCGAACAGCGCGGTATTGTCGGTCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGtttt                                                                                                       <  1:348093/249‑1 (MQ=255)                                                                                                                                                  ggTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCtt                                                                                                  >  1:56317/1‑251 (MQ=255)                                                                                                                                                           cccGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCCGTTCTTTGATGACCTGAACACGCCCGATTTCGACCGCTGTTTTTCCTTTATGCAGGc                                                                                         >  1:41776/1‑251 (MQ=255)                                                                                                                                                                          ggCGGTGGCTTCGACTTATCCCCATTCTATGGTTTTGAAGAAGAGTCTGTGCATTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCt                                                                          <  1:307908/250‑1 (MQ=255)                                                                                                                                                                            cGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTTTGCAGGCGGGAGG‑GAAAGGccac                                                                        >  1:37609/1‑247 (MQ=255)                                                                                                                                                                                 gCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAGCGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGGAGG‑CAAAGGCTACACCGAc                                                                  >  1:214730/1‑251 (MQ=255)                                                                                                                                                                                           aaCCCCATTCTATGGTTTTGAAGAAGAAGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTTATGACCTGAACACGCCAGAATTCGACCGCTGTTTTTCCTTTATGCCGGCGGTAGGGCAAAgctaaacccacggtttttttt                                                         >  1:278435/1‑231 (MQ=255)                                                                                                                                                                                             ccccATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTTATTACCTGAACAAGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGGAGG‑GAAAAGCCACACCCACCCCTATTTAccc                                                      >  1:205750/1‑250 (MQ=255)                                                                                                                                                                                             ccccATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCTACACCGACGCTTATTTACCa                                                      <  1:35812/251‑1 (MQ=255)                                                                                                                                                                                                                  aagaTGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCCGTTCTTTGATGACCTGAAAACGCCAGATTTCGACCGCTTTTTTGCCTTTATGCCGGGGGTAGG‑CAAAAGCTACACCGCCGCCTTTTTACCCATTGTCTAGCGGCGGAAAGCg                                 >  1:198897/1‑251 (MQ=255)                                                                                                                                                                                                                     aTGCTATTCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATg                              <  1:143563/251‑1 (MQ=255)                                                                                                                                                                                                                            tCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAGCGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATG‑GCCTAc                       <  2:214730/250‑1 (MQ=255)                                                                                                                                                                                                                            tCACTGGCATCGCACCGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATG‑GCCTAc                       <  2:37609/250‑1 (MQ=255)                                                                                                                                                                                                                                           cGCCCGTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCCGTTCTTTTATTACCTTAACACCCCAGAATTCGACCGCTGTTTTTCCTTTTTGCATGCGGTGGG‑CAAAAGCTACAACCACCCCTTTTTATCAATTGTCGCACGCCGGAAAGCGATG‑GCCTACGgcgcgcgccaacgcc       >  2:5106/1‑244 (MQ=255)                                                                                                                                                                                                                                                gTGACCTGCGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGGTAACCCGACGCTTTTTTACCAATTTTTGCGCGACGGAAAGCG‑TGTGCCTACGGCGCGCGGGAGGGCAATTTc  >  1:137163/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGG‑CAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATG‑GCCTACGGCGAGCGCGAGCGCAATTTC  >  NC_000913/2553431‑2553920

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