Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,809	A→G	K252R (AAA→AGA)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,809	0	A	G	100.0%	180.9 / NA	53	K252R (AAA→AGA)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/32);  total (21/32)

AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCG  >  NC_000913/498562‑499049                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 aaTGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAAt                                                                                                                                                                                                                                               >  1:158652/1‑251 (MQ=255)  aTGCACCGGCGATCAGCGAACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCAACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGATGCCGCACACCGACTCGTGAACTGGCTTCCGCATTGGGGCTGGCACCGGAAAAAGTGATGATGCCCTTTCAGCCGCGATTTGGTCGGGACCCCTGGCTGATGCCTTATACCGACGAACCGATGAGAATg                                                                                                                                                                                                                                              >  2:419189/1‑251 (MQ=255)    gCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCt                                                                                                                                                                                                                                            >  1:311554/1‑251 (MQ=255)             aaCAGCGTACGCGCTTCTTCTGCCAAACATGGAGAAGCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGCTGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGCCGAAACGCTGAGAATGCTCGgagcaa                                                                                                                                                                                                                                    >  2:57443/1‑247 (MQ=255)                aGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGTGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAgg                                                                                                                                                                                                                                <  2:196452/251‑1 (MQ=255)                aGCGTACGCGCTTCTTTGGCCAAACAGGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAgg                                                                                                                                                                                                                                <  2:128544/251‑1 (MQ=255)                                           ggCGAACCGGATCTGCTACTGCTCTCTTATCATGGCCTTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCACCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGCTGGCACCGGAAAAAGCGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGACAAAGGCGTAGGTCATATTCAGGTGATGTGccc                                                                                                                                                                                                     >  2:220653/1‑251 (MQ=255)                                               aaCCGGATCTGCTAATGCTCTATTATAATGGAATTCCCCAGCGTTATGAAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGACGGGAACCCTGGCTGATGCCTTATACAGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGc                                                                                                                                                                                                 >  1:161263/1‑251 (MQ=255)                                                aCCGGATCTGCTACTGCTCTCTTCTCATGGCATTCCCCAGCGTTATGCAGAAGCAGGCGATGATTACCCGCGCCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGAATTGGGGCTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGCGAATGCTCGGAGAAAAAGGCGTAGGGCATATTCAGGTGATGTGCCCGGGCt                                                                                                                                                                                                >  2:451415/1‑251 (MQ=255)                                                  cGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATGACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGGCTAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCttt                                                                                                                                                                                              <  1:220653/251‑1 (MQ=255)                                                   ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtg                                                                                                                                                                                             >  2:144591/1‑251 (MQ=255)                                                       cTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGCGCCCGGGCTTtgctgc                                                                                                                                                                                         >  1:46194/1‑251 (MQ=255)                                                             cTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTg                                                                                                                                                                                   >  2:309423/1‑251 (MQ=255)                                                               gCTCTCTTAGCATGGCAGTCCCCAGCGTTAGGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTGTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGACTTGGCTGCGGATTGt                                                                                                                                                                                  <  2:232385/250‑1 (MQ=255)                                                                  ctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGCTTACCCGCAACGTTGCCGCACAACGCCTCGTGAACTGGCTTCCGCCTTGGGGATGGCACCGGAAATAGTGATGATGCCCTTTCAGTCGCGCTTTGGTCGGGAACACTGGCTGATGCCTTATACCGACGAAACGCTGAGACTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTgg                                                                                                                                                                              >  1:427722/1‑251 (MQ=255)                                                                       aTCATGGCATTCCCCAGCGTTAGGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGGTCGGAGAAAAAGGCGTAGGTCATATTCAGGTTATGTGCCCTGGCTTGGCTGCGGAGGGTCTGGAGAc                                                                                                                                                                          <  2:292500/250‑1 (MQ=255)                                                                       aTCATGGCATTCCCCAGCGTGATGCAGATGAAGGCGATGATTACCCGTAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTGTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGGTCGGAGAAAAAGGCGGAGGTCATAGTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGAc                                                                                                                                                                          <  1:401596/250‑1 (MQ=255)                                                                           tGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTgg                                                                                                                                                                     >  1:153580/1‑251 (MQ=255)                                                                                 tCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGCGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAagag                                                                                                                                                                <  2:46194/250‑1 (MQ=255)                                                                                   cccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCCCAACGACTCGTTAACTGGCTTCCGCATTGGGGATTGCACCTGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGGTTTGCTGTGGATTGTCGGGAGACGCTGGAAGAGAtt                                                                                                                                                             <  1:176339/251‑1 (MQ=255)                                                                                   cccAGCGTTATGCAGAGGAAGGCGATGAGTACCCGCCACGTTGCCGCACAACGACTCGTGAACTGGCTGCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTGTCAGGCGCGCTTTGGTCGGGAACCCTGGCGGATGCCTGAGACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGGCATATTCAGGTGATGGGCCCGGGCTTTTCTGCGGATTGTCTGGAGACGCGGGAAGAGAtt                                                                                                                                                             <  2:427722/251‑1 (MQ=255)                                                                                            aTGCAGAGGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCaa                                                                                                                                                    <  2:311554/251‑1 (MQ=255)                                                                                                    gAAGGCGATGATTACCCGCAATGTTGCCGCACAACGACTCGTGAACTGGCGTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGTGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTg                                                                                                                                             <  2:99231/250‑1 (MQ=255)                                                                                                          gatgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCtt                                                                                                                                      <  1:421995/251‑1 (MQ=255)                                                                                                                tACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCAGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCg                                                                                                                                 <  1:309540/250‑1 (MQ=255)                                                                                                                     gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTTGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcg                                                                                                                           <  1:266732/251‑1 (MQ=255)                                                                                                                      cAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcgg                                                                                                                          <  2:96010/251‑1 (MQ=255)                                                                                                                            tGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCgggg                                                                                                                     >  1:293548/1‑249 (MQ=255)                                                                                                                            tGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGa                                                                                                                     <  2:132357/250‑1 (MQ=255)                                                                                                                                 cacaACGACGCGGGGACTGGCTTCCGCATTGGGGATGGCACCTGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATTGCTTATACCTACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCGCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGATTAGAGTGCCGAGCAAAACCGATAGGGCTGCCTCGGTGCCGGCGGGAAAAAat                                                                                                               <  2:295480/251‑1 (MQ=255)                                                                                                                                  acaACGACCCGTGAACTGGCTTCCGCCTTGGGGATGGCCCGGGAAAAAGCGATGATGACCTTTCAGTCGAGCTTTGGTCGGGAACCCTGGCTGATGACTTATACCGACGAAACGCTGAGACTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGCTTGCCGCGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCgggggataaa                                                                                                               >  2:430551/1‑243 (MQ=255)                                                                                                                                       gACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGGAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTAGACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGTTGAGGTGCCCTGGCTTTGCTGCGGAGTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAat                                                                                                         <  1:297884/251‑1 (MQ=255)                                                                                                                                              gAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATAGTCAGGTGATGTGCCCGGGCTGTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggc                                                                                                                         <  2:293094/228‑1 (MQ=255)                                                                                                                                              gAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAAGGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTcc                                                                                                  <  1:430551/251‑1 (MQ=255)                                                                                                                                              gAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGACGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGATGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggc                                                                                                                         >  1:293094/1‑228 (MQ=255)                                                                                                                                                    gCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGATTTGGTCGGGAACCCAGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCCAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAGCAAATATGAATATAtttcgctgct                                                                                            >  2:96525/1‑243 (MQ=255)                                                                                                                                                          gCATTGGGGATGGCACCGGAAAAAGTGATGATGCCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGCTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAACACATGCATAAATTACGGCTCCTAATGc                                                                                      >  2:382917/1‑251 (MQ=255)                                                                                                                                                          gCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                      <  1:389098/251‑1 (MQ=255)                                                                                                                                                              tGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCAc                                                                                   <  1:27782/250‑1 (MQ=255)                                                                                                                                                                       cACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACat                                                                         <  1:378244/251‑1 (MQ=255)                                                                                                                                                                                    gatgatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCGCTTAATGCCACGCCGGAACatatagcaaagatgg                                                            >  1:431152/1‑240 (MQ=255)                                                                                                                                                                                           aCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCGCGTTATGCCCCGCCGGGACATATCGAAATGATgggttaact                                                     >  2:174320/1‑244 (MQ=255)                                                                                                                                                                                                  aggcgcgCTGTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAGTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGGGCCGGCGGTAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGTAACATATCTAAATGATGGCTAATCTTGTTGcc                                              <  2:452277/248‑1 (MQ=255)                                                                                                                                                                                                           ttGGTCTGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATc                                       <  1:261083/249‑1 (MQ=255)                                                                                                                                                                                                            tGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCt                                    >  1:377558/1‑251 (MQ=255)                                                                                                                                                                                                            tGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGc                                     <  1:255570/250‑1 (MQ=255)                                                                                                                                                                                                                 gggAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAGATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                               <  1:216427/251‑1 (MQ=255)                                                                                                                                                                                                                    aaCCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGCTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGACGGCGGGGAAAAATATGAAAATATTCAGGGGCTTAATGGCCCGCCGGAGAATAACGAAATGAAGGCTCAACttgtttgcgcggttcggcaaaacttg                            >  2:33256/1‑230 (MQ=255)                                                                                                                                                                                                                       ccTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCg                         <  1:257679/251‑1 (MQ=255)                                                                                                                                                                                                                                  ccTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTg              <  1:273136/251‑1 (MQ=255)                                                                                                                                                                                                                                      ataCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgc           <  2:153580/250‑1 (MQ=255)                                                                                                                                                                                                                                              gaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgcg  <  1:78553/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGa                                                                                                           <  1:446632/144‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCG  >  NC_000913/498562‑499049

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