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%	212.2 / NA	61	K252R (AAA→AGA)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/34);  total (27/34)

AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGT  >  NC_000913/498562‑499045                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             aaTGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAAt                                                                                                                                                                                                                                           <  1:93570/251‑1 (MQ=255)    gCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATTCCGACGAAACGCTGAGAATGCt                                                                                                                                                                                                                                        <  1:130175/251‑1 (MQ=255)     cACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTc                                                                                                                                                                                                                                       <  2:100742/251‑1 (MQ=255)                  cGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCAGTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCg                                                                                                                                                                                                                          <  2:223329/251‑1 (MQ=255)                   gTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGt                                                                                                                                                                                                                         >  1:33235/1‑251 (MQ=255)                       gcgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGt                                                                                                                                                                                                                     >  1:62294/1‑251 (MQ=255)                        cgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGTCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTTCCCGCAACTTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAATTGATGATGACCGTTCAGTCTCTCTTTGGTCTGGATCCCTGGCTGATGGCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTc                                                                                                                                                                                                                    <  1:187281/251‑1 (MQ=255)                         gcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCGCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGATGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCa                                                                                                                                                                                                                   <  2:335312/251‑1 (MQ=255)                            tcttTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGCTGATTCCCCGCCACGTTGACGCACAACGACTCGTGCACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCatat                                                                                                                                                                                                                >  1:383472/1‑251 (MQ=255)                            tcttTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCGGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCat                                                                                                                                                                                                                  <  1:375606/249‑1 (MQ=255)                            tctggtGCCAAACATGGCGAACCGGATCTGCTACTGCTATCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCGCGGAGAAAAAGGCGTAGGTCata                                                                                                                                                                                                                 <  2:322642/245‑1 (MQ=255)                              ttttGCCAAACATGGCGAACCGGCTCTGCTACTCCTATCTTCTAACGGCATTCCCCAGCGTTATACAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGCTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGCTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCCTAttt                                                                                                                                                                                                              >  2:254045/1‑250 (MQ=255)                               tttGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCatat                                                                                                                                                                                                                <  1:84833/248‑1 (MQ=255)                               tttGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCAGTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCATTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCGCGGAGAAAAAGGCGTAGGTCATAtt                                                                                                                                                                                                               <  1:311641/249‑1 (MQ=255)                                         aTGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCCTTCCCCCACGTTATGCAGATGAAGGCGATGATTCCCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGACTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGc                                                                                                                                                                                                   >  1:401009/1‑251 (MQ=255)                                           ggCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATGCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCGCGGAGAAAAAGGCGTAGGTCGTATTCAGGTGATGTGccc                                                                                                                                                                                                 <  2:276323/251‑1 (MQ=255)                                             cGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCCACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCgg                                                                                                                                                                                               >  2:379933/1‑251 (MQ=255)                                                   ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCt                                                                                                                                                                                            >  2:139873/1‑248 (MQ=255)                                                   ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCt                                                                                                                                                                                            <  1:139873/248‑1 (MQ=255)                                                   ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtg                                                                                                                                                                                         >  1:143514/1‑251 (MQ=255)                                                   ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAAGGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtg                                                                                                                                                                                         >  2:299187/1‑251 (MQ=255)                                                         gCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCg                                                                                                                                                                                    <  2:143514/250‑1 (MQ=255)                                                          cTACTGCTCTCTTATCATGGCATTACCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTACGCATTGGGGATGGCACCGGACAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATCTTCAGGTGATGTGCCCGGGCTTTGCTGCgg                                                                                                                                                                                   >  2:248138/1‑250 (MQ=255)                                                           tACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAt                                                                                                                                                                                 >  1:286674/1‑251 (MQ=255)                                                            aCTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCACACGACGGCGATGACTACCCGCAACGTTACCGCACACCGACTCGCGAACTGGCTTCCGCATTGGGTATCGCACCGGCAAAAAAGCTGATGACATTTCCGTCGCGCCATGGTCGGCAACCCTGGCTGATGCCTTATACCGACGAAACCCTGAGCATGCTCGGAGACAAAGGCGTAGGCCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGaat                                                                                                                                                                                >  1:179434/1‑249 (MQ=255)                                                                    ctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGAGGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAATGCGGAGGTCAGATTGAGGTGGTGTGCCCGGGCTTGGCTGCTGATTGTCTgg                                                                                                                                                                          <  2:267832/249‑1 (MQ=255)                                                                      tATCATGGCATTCCCCAGCGTTAGGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCAGTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGgag                                                                                                                                                                        <  2:234822/249‑1 (MQ=255)                                                                          aTGGCATTCACCAGAGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTg                                                                                                                                                                  >  1:88330/1‑251 (MQ=255)                                                                               aTTCCCCAGCGTTATGAAGATGAAGGCGATGATTCCCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCAGCATTGGGGCTGGAACCGGAAAAAGCGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCCGATGCCTTATACCGACGAAGAGCTGAGAATGCTCGGCGAAAAAGGCGTAGGTCATCTTCAGGTGATGTGCCCGGGCTTCGCTGCGGATCGTCTGAAGACGCTGGAagg                                                                                                                                                             >  2:159742/1‑250 (MQ=255)                                                                                           tATGCGGATGAAAGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGAGGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATAGTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGc                                                                                                                                                  <  2:365294/250‑1 (MQ=255)                                                                                             tGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCaa                                                                                                                                                >  2:396252/1‑250 (MQ=255)                                                                                                aGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAAcc                                                                                                                                            <  1:258499/251‑1 (MQ=255)                                                                                                 gATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCg                                                                                                                                           <  2:62294/251‑1 (MQ=255)                                                                                                  aTGAAGGCGATGATTACCCGCAACGCTGCCGCACAACGCCTCGTGAAATGGCTTCCGCATTGGGGATGGCACCGGAACAAGTGCTGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGCCGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGACGAGATTGCCGAGCAAACCCGt                                                                                                                                          >  2:74102/1‑251 (MQ=255)                                                                                                          gatgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGATCCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCGGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTcat                                                                                                                                  <  2:155919/251‑3 (MQ=255)                                                                                                             gatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTGCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTAGACCGACGAAACGCTGAGAATGCTCGGAGAAAAGGGCGGAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATGGCCGAGCAAAACCGTGAGGTCTTCCt                                                                                                                               <  1:159742/251‑1 (MQ=255)                                                                                                                       aaCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggc                                                                                                                     >  2:3874/1‑251 (MQ=255)                                                                                                                                                          gCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                  <  1:379933/251‑1 (MQ=255)                                                                                                                                                           cATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCCGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAtt                                                                                                                                                         >  2:311376/1‑179 (MQ=255)                                                                                                                                                           cATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAtt                                                                                                                                                         <  1:311376/179‑1 (MQ=255)                                                                                                                                                           cATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGTGGAAAAatattactttattcctcggccttatttc                                                                                 >  2:126815/1‑227 (MQ=255)                                                                                                                                                           cATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAAGATATTCCGGCGCTTAATGcc                                                                                 <  1:97582/251‑1 (MQ=255)                                                                                                                                                           cATTGAGGATGGCACCGGAAAAAGTGGGGATGGCCTTGCAGTCGCGGTTGGGTCGGGAACCCTGGCTGATGCCTTAAACCGTAGAAGCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCGGGTGATAGGCCCGGGCGTTGCTGCGGATGGGCGGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTGTTCCTCGGTGCCGGGGGGAGAAAATACGAATGTATTCCGGCGCTTTATGcc                                                                                 <  2:432189/251‑1 (MQ=255)                                                                                                                                                              gtgggATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTGTGGTCGGGAACCCTGGCTGATGCTTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAATGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGTGGATTGTCTGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTCTTCCTCGTTTCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCAc                                                                               <  1:298103/248‑1 (MQ=255)                                                                                                                                                                gggATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGcc                                                                            <  1:349625/251‑1 (MQ=255)                                                                                                                                                                        aCCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGCAACCCTGGCTGATGCCTTATAGCGACGACACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGAGCCCGGGCTTTGCTGCGGCTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATATATTCCCGCGGTTTATTCCCCGCCGGcacac                                                                     >  1:214185/1‑249 (MQ=255)                                                                                                                                                                              aaaaGTGATGATGACCTTTACGTCGCGCTTTGGTCGGGAACCCTGGCTGATGACTTATACAGACGAAACGCTGCGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGACCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAACATGAATATATTTCGGCGGTTAATGCCCCCACGGGACATATCGaa                                                               >  2:260023/1‑250 (MQ=255)                                                                                                                                                                                    gatgatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGGCACGCCGGGACATATCGAAAAGATgg                                                        >  1:304741/1‑251 (MQ=255)                                                                                                                                                                                     atgatgaCCTTTCAGTCGCGCTTTGGACGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTCCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATGGc                                                       >  1:368030/1‑251 (MQ=255)                                                                                                                                                                                          gaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAt                                                   <  2:126770/250‑1 (MQ=255)                                                                                                                                                                                           aCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCTAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCt                                                 <  2:304741/251‑1 (MQ=255)                                                                                                                                                                                               ttCAGTCGCGCTGTGGTCGGGAGCCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTGTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCGTCCTCGGTGCCGTCGGGAAAAAATATGAATATATTCCGGCGCTGAATGCCACGCCTGAACATATCGAAATGATGGTTAATCttgt                                              <  2:214185/250‑1 (MQ=255)                                                                                                                                                                                                tCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCGTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgtt                                             <  2:380481/250‑1 (MQ=255)                                                                                                                                                                                                     cgcgCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGGCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAATCCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATAGGAATATATTCCGGCGCGTAAGGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCcgc                                        <  1:192623/250‑1 (MQ=255)                                                                                                                                                                                                     cgcgCTTTGCCCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAAGCGTAGGTCATATTCAGGTGATGGGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGTAAGAGATTGCCGAGCAAAACCGTAAGGTCTTCCTCGGTGCTGGCGGGAAAAAATATGAATATAGTCCGGCGCTTGATGCCACGCCGGATCATATCGAAATGATGGCTAACCTCGTTGCcgcg                                       <  2:399458/251‑1 (MQ=255)                                                                                                                                                                                                              gTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAAAATATCGAAATGATGGCTAATCTTTTTGCCGCGGATCGCTa                               >  1:96230/1‑250 (MQ=255)                                                                                                                                                                                                                          ggCTGATGCCTTATACCGACGAAACGCTGAGCATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGACAAAATATGAAAATATTTCGGCCGTTTATGCCACGCCGGGCCATATTGAAAAGATAGGTAATCTTGTTTCcgcggactgctaaagatgcgcgggg                  >  2:153479/1‑230 (MQ=255)                                                                                                                                                                                                                          ggCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGCGAACAAGGCGTCCGTCATATTCCGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGGGCggg                   >  1:29696/1‑249 (MQ=255)                                                                                                                                                                                                                                     tataCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGTAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAAGGATGGCGAATCTTGTTGCGGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGa         <  2:202893/249‑1 (MQ=255)                                                                                                                                                                                                                                        aCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgcgc     <  2:286674/250‑1 (MQ=255)                                                                                                                                                                                                                                          cgacgaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCGGCGCTTAATGCCACCCCCGGACATATCGAAATGATGGCTAATCTTGTTTCCGCGTATCGCTAAAGCTGAGGGGTAAAGAAcctatcggcgc  >  2:257560/1‑241 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGT  >  NC_000913/498562‑499045

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