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

CACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCG  >  NC_000913/498566‑499049                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 cACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCTAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGTCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTATACCGACGAAACGCTGAGAATGCTc                                                                                                                                                                                                                                           <  2:447515/251‑1 (MQ=255)  aCTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCg                                                                                                                                                                                                                                          >  1:129289/1‑251 (MQ=255)         aaCAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGaaaa                                                                                                                                                                                                                                   <  1:114801/251‑1 (MQ=255)            aGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTCTACCGACGAAACGCTGAGAATGCTCGGAGAAAAAgg                                                                                                                                                                                                                                >  2:269570/1‑251 (MQ=255)               gTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCgg                                                                                                                                                                                                                             >  2:188580/1‑250 (MQ=255)                   gcgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGt                                                                                                                                                                                                                         >  2:352542/1‑251 (MQ=255)                   gcgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGt                                                                                                                                                                                                                         <  1:188580/251‑1 (MQ=255)                            ttGCCAAACATGGCGCACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTc                                                                                                                                                                                                                  >  1:400653/1‑249 (MQ=255)                              gCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGt                                                                                                                                                                                                              <  1:186233/251‑1 (MQ=255)                                   aCATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGAtgt                                                                                                                                                                                                         >  1:353911/1‑251 (MQ=255)                                        gCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGAAGGCACCGGAAAAAGAGATGATGAACTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGCAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCg                                                                                                                                                                                                    >  1:278928/1‑251 (MQ=255)                                         cGAACCGCATCTGCTACTGCTCTCTTATCATGGCATCCCCCCGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGCCTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGCAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGCAGGCCATATTCAGGGGATGTGccccg                                                                                                                                                                                                   >  2:30917/1‑249 (MQ=255)                                          gAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGTTGAGAATGCTCGGAGAAAACGGCGTAGGTCATATTCAGGTGATGTGCCCggg                                                                                                                                                                                                  >  2:406108/1‑251 (MQ=255)                                              cGGATCTGCTACTGCTCTCTTATCCTGGCAGGCCACGGCGTTATGCAGATGAAGGCTATGATTACCCGCAACGTTGCAGCACAACGACTCGTGAACTTGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATTGCCGTACGGTCGCGCTTTGGTCGTGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCttt                                                                                                                                                                                              <  2:20127/251‑1 (MQ=255)                                              cGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCttt                                                                                                                                                                                              <  1:33176/251‑1 (MQ=255)                                               ggATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACACTGGCTGATGCCTTATACCGCCGAAACGCTGCGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtg                                                                                                                                                                                             >  2:258607/1‑251 (MQ=255)                                                gATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtgc                                                                                                                                                                                            >  2:222220/1‑251 (MQ=255)                                                       tACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAt                                                                                                                                                                                     >  1:123775/1‑251 (MQ=255)                                                       tACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGATTCCGCATTGGGGATGGCACCGGAAAAAGTGATAATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTCTACCGACGCAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTAAGGTGATGTGCCCGGGCTTTGCTGCGGAt                                                                                                                                                                                     >  1:307989/1‑251 (MQ=255)                                                               tctTATCACGCCATTCCCCAGCCTCATGCAGATGTCCGCGATGATTCCCCACAACGTTGCCGCACAACCACCCGCGAACTGGCTTCCCCATTGGGGATGCCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGCGCATGCTCGGAGAAAAAGGCGTAGGTCTTATTAAGGTGATGTGCCAGGGCTTTGCTGCGGATTGTCTgg                                                                                                                                                                              >  2:51494/1‑250 (MQ=255)                                                                ctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTgg                                                                                                                                                                              <  2:353911/249‑1 (MQ=255)                                                                  tATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGgag                                                                                                                                                                            <  1:106040/249‑1 (MQ=255)                                                                   aacATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAGTGTCTGGAGAc                                                                                                                                                                          <  2:38717/248‑1 (MQ=255)                                                                    tCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGc                                                                                                                                                                        >  2:126321/1‑251 (MQ=255)                                                                          cATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGAGGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTTAGAATGCTCGGAGAAAAAGGCGTAGGTCAGATTCAGGTGATGTGCCCGGGCTTTTCTGCGGTTTGTCTGGAGATGCTGGAag                                                                                                                                                                  <  2:58945/251‑1 (MQ=255)                                                                               cccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAtt                                                                                                                                                             <  1:26802/251‑1 (MQ=255)                                                                               cccAGCGGTATGCATATGAAGGCGATTATTACCCACAGCGTTGCCGCACAACGACTCGTGAACTCGCTTCCGCATTGGGGATGGCACCGTAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTTGCTGATGCCTTATTCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCTTAGGTCATATTCAGGTTATGTGCCTGGGCTTTGCTGCGGATTGTCTGGAGACGCTTGAAGAGAtt                                                                                                                                                             <  1:157653/251‑1 (MQ=255)                                                                                                 aaGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGAc                                                                                                                                                                          <  2:51636/220‑1 (MQ=255)                                                                                                 aaGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGAc                                                                                                                                                                          >  1:51636/1‑220 (MQ=255)                                                                                                     cGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGCACTGGCTTCCGCATTGGGGTTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGtgt                                                                                                                                       >  2:148217/1‑249 (MQ=255)                                                                                                      gatgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCtt                                                                                                                                      >  1:468459/1‑251 (MQ=255)                                                                                                       atgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTc                                                                                                                                     <  2:478669/251‑1 (MQ=255)                                                                                                       atgatTACCCGCAACGTTGACGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTATGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTTCCGAGCAAAACCGTGAGGTCt                                                                                                                                       <  2:449004/249‑1 (MQ=255)                                                                                                          atTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTc                                                                                                                                  <  2:96076/251‑1 (MQ=255)                                                                                                                   aaCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggc                                                                                                                         <  1:102381/251‑1 (MQ=255)                                                                                                                    aCGTTGCCGCACAACGACTCGTGAACTGGCTTCAGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggcg                                                                                                                        >  2:222976/1‑251 (MQ=255)                                                                                                                     cGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggcgg                                                                                                                       <  1:396347/251‑1 (MQ=255)                                                                                                                      gTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCggg                                                                                                                      <  2:228899/251‑1 (MQ=255)                                                                                                                        tGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGa                                                                                                                     <  1:15239/250‑1 (MQ=255)                                                                                                                        tGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGa                                                                                                                     <  1:134021/250‑1 (MQ=255)                                                                                                                         gCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGACTTATACCGACGAAACGCTGAGAATGCTCGGAGAACAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGACGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCggggg                                                                                                                    >  2:85455/1‑248 (MQ=255)                                                                                                                               caACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGATTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAatat                                                                                                             >  1:326146/1‑251 (MQ=255)                                                                                                                                aaCGACTCGGGAACTGGCTTCCGCATTGGGGATGGCCCCGGAAAAAGTGATGATGACCTGTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAGTTCCGAGCAAAACCGTGAGGTGTTCCTCGGTGCCGGCGGGAAAAAATATg                                                                                                            <  2:455561/251‑1 (MQ=255)                                                                                                                                   gACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAACAAGTGATGATGCCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAACAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCggggaggagtaagga                                                                                                          >  2:29991/1‑238 (MQ=255)                                                                                                                                              tGGCTTCCGCATTGGGGATGGCACCGGCCAAAGTGATGATGCCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAACCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCCGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAATATTTCGgcg                                                                                              >  2:477059/1‑251 (MQ=255)                                                                                                                                               ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGACAAAATAAGAatatttttcggggg                                                                                             >  2:245291/1‑241 (MQ=255)                                                                                                                                                gCTTCCGCTTTGGGGATGGCACCGGAACAAGCGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGACACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGCGACGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGAGCCGGCGGGTAACAATATGAATATAtttagccgct                                                                                            >  2:297036/1‑243 (MQ=255)                                                                                                                                                    ccGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGTTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAt                                                                                        <  1:406108/251‑1 (MQ=255)                                                                                                                                                            gggATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCTCTTAATGCCACgac                                                                                <  2:123930/251‑3 (MQ=255)                                                                                                                                                               aTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAACCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCggg                                                                             >  1:305885/1‑250 (MQ=255)                                                                                                                                                                            aaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGCTGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGGCAAAAATTAATATAATCCTCCGCTTAATTCCCCGCCGGGACatatacaaat                                                                 >  1:286719/1‑244 (MQ=255)                                                                                                                                                                                gatgatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGCGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATATATTTCGGGGCTTTATGCCACCCCgtacaatatcgacaagatgg                                                            >  2:207005/1‑232 (MQ=255)                                                                                                                                                                                 atgctgtccTTTCAGTCGCGCTTTGTTCGGGAACCCTGGCTGATGCCTTATACCTACGAAACGCTGAGAATGCTCTGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGGCCGTGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAATCCGTGAGGTCTGCCTCGGTGCCGGCGGGAAAAAATATGAGTATATTCCGGCGCTTAAGGCCCCGCCTGAATATATCGAAATGATGGc                                                           <  2:28183/244‑1 (MQ=255)                                                                                                                                                                                   gatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGTAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCt                                                          <  1:183482/250‑1 (MQ=255)                                                                                                                                                                                   gatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATAGTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCt                                                          <  2:39128/250‑1 (MQ=255)                                                                                                                                                                                          tttCAGTCGGGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAGAAAATATGAATATATTTCGGCGCGTAATGGCAACCCGGGACATAACGAAATGATggatcatcctgt                                                  >  2:31876/1‑241 (MQ=255)                                                                                                                                                                                                       ttGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATc                                       <  1:245291/249‑1 (MQ=255)                                                                                                                                                                                                          gTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCGAATCTTGTTGCCGCGTATCGCTaa                                  <  1:178019/251‑1 (MQ=255)                                                                                                                                                                                                             gggAACCCTTGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTGGTTGCCGCGTATCTCTAAAGc                               <  2:200515/251‑1 (MQ=255)                                                                                                                                                                                                             gggAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                               <  1:31876/251‑1 (MQ=255)                                                                                                                                                                                                                   ccTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGACAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGCGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGCAAATATTTCCGCGCCTTATTCCACGCCGGGAAACATCGAAAAAATGGCTAATCTTGTTGTCGCGTATCGCTCACGCTGAGCg                         >  2:175481/1‑251 (MQ=255)                                                                                                                                                                                                                      ggCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTTCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAGTGATGGGTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGt                       <  1:51494/250‑1 (MQ=255)                                                                                                                                                                                                                          gATGCCTTATACCGACGAAACGCTGAGCATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCGGGGCTTTATGCCACGCCGGGAAATATCGCAATGAAGGCTAATCTTGTTGGCGCGTATTGCTAAAGCTGAGCgggaa                     >  1:168411/1‑245 (MQ=255)                                                                                                                                                                                                                           gtGCCTTATACCGCCGAAACGCTGAGAATGCGCGGAGGAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATGGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAGTATATTCCGGCGCTTAATGCCACGCCGGAACGTATCGAAAGGATGGCTAAGCGTGTTGCCGCGTATCGCTAAAGCTTAGCGGTAAAGaa                 <  2:93804/250‑1 (MQ=255)                                                                                                                                                                                                                                    aCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGcc        <  1:5273/251‑1 (MQ=255)                                                                                                                                                                                                                                     ccGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCg       <  2:123775/251‑1 (MQ=255)                                                                                                                                                                                                                                       gacgaAACGCTGAGAATGCTCGGATAAAAAGGCGTAGGTCATATTCAGGTGATGTGCGCGGGCTGTCCTGCGGATTGGCTGGAGACGCTGGGAGAGATTACCGAACAAAACCGTGATGTCTTCCGCGGTGCCGGCGGTAAAAAATATGAATATATTCCGTCGCGTAATGCCACGCCGTAACATATCGAAATGATGGCTAAGCTTGGTGCCGCGTATCTCTAAAGGTGAGCGGTTAGGAACTTATCGCCGTc     <  2:200516/251‑1 (MQ=255)                                                                                                                                                                                                                                       gacgaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAGGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTc     <  1:477059/251‑1 (MQ=255)                                                                                                                                                                                                                                          gaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCGCTTAATTCCACGCCGGGACATAACGAAATGATGGTTAATCTTGTTGGCGCGGATCGCTAAAGCTGAGCGGGAAAGAACTGAGGGCCGTcgcg  >  2:288968/1‑251 (MQ=255)                                                                                                                                                                                                                                              cGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATAAATTTCGGGGCTTAATGCCACGCCGGGACATATCGAAATGAAGGCTTATCTTTTTGCCGCGTATCGCTaa                                  >  2:470853/1‑215 (MQ=255)                                                                                                                                                                                                                                              cGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTaa                                  <  1:470853/215‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 CACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCG  >  NC_000913/498566‑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: ‑