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

AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAA  >  NC_000913/498562‑499051                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    aaggcgcTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAGCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCTCTTTGGTCGGGAACCCTGGCTTG‑TGCCTTATACCTACGAAACGCTGAGAAt                                                                                                                                                                                                                                                  <  2:609475/245‑1 (MQ=255) aaTGCGCTGGCGAACAGCGTACGCGCTTCTTTTGACAAACATGGCGAACCCGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGATGGCGATTATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAAt                                                                                                                                                                                                                                                  <  2:432798/251‑1 (MQ=255) aaTGCACTGGCGAACAGCCTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTTATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAAt                                                                                                                                                                                                                                                  <  1:510255/251‑1 (MQ=255)   tGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGc                                                                                                                                                                                                                                                >  1:181977/1‑251 (MQ=255)         ggCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTTGCTTCCGCATTGTTGATGGCACCGGAAAATGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGTAACCCTGTCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGgag                                                                                                                                                                                                                                          <  1:71633/251‑1 (MQ=255)         ggCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCgg                                                                                                                                                                                                                                            <  1:542968/249‑1 (MQ=255)              accgctgaCGCGCTTCGTTAGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGAGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGaaaaa                                                                                                                                                                                                                                     <  2:148121/244‑1 (MQ=255)               cAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATTCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAg                                                                                                                                                                                                                                    <  1:172684/251‑1 (MQ=255)                    tACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGt                                                                                                                                                                                                                                >  1:593896/1‑250 (MQ=255)                     aCGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGCGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGt                                                                                                                                                                                                                                >  2:143062/1‑249 (MQ=255)                       gcgcTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGCATGCTCGGAGAAAAAGGCGTAGGt                                                                                                                                                                                                                            >  1:495514/1‑251 (MQ=255)                        cgcTTCTTTTGCCAAACATGGCGAACCGGATCGGCTACTGCTCCATTATCATGGCATTCCCAAGCGTTATGCAGATGAAGGCTAGGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGGATTGGTGATTGCACCGGAAAAAGTGATGATGACCTTTCAGGCGCGCTGTGGTCGGGAACCCTGGCTGA‑TGCCGTATGCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTc                                                                                                                                                                                                                           <  1:261550/251‑1 (MQ=255)                            tcttTTGCCAACCATGGCGCACAGGATCTGCTACTGCTCTCTTATCATGGCATTACCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGACGCACAACGCCTCGTGAACTGGCTTCCGCATTGGGGATGGCACGGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAGCGCTGAGAATGCTCGGGGAAAAAGGAGTAGGTCctat                                                                                                                                                                                                                       >  2:242730/1‑251 (MQ=255)                                   ccacaCATGGCGAACCGGATCTGCTACTGCTCGCTGATCATGGCCTTCCCCAGCGTTCTGCAGATGAAGGCGATGATTACCCGCAACGTCGCCGCCCAACGAGCCGTGAACTGGCTTCCGCCTTGGGGATGGCAGCGGAAAAAGCGCGGATGATCTTGCAGTCGCGTTTTGGTCGGGAACCCGGGTTGA‑TGCCTTATATCGACGACACGCTTAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTg                                                                                                                                                                                                                <  1:348327/247‑1 (MQ=255)                                     aaaCATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGCTTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGAt                                                                                                                                                                                                              >  2:261490/1‑251 (MQ=255)                                          gggCGAACCGAATCTGCTACTCCTCCCTTATCATGGCAGTCCCCAGCGTTATGCAGATGCAGGCGATTATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGTCCCCGGACAAAGTGATGATGACCTTGCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGcc                                                                                                                                                                                                         <  2:324545/250‑1 (MQ=255)                                           ggCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGccc                                                                                                                                                                                                        <  2:129738/251‑1 (MQ=255)                                                aCCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCt                                                                                                                                                                                                   >  2:299676/1‑251 (MQ=255)                                                    gATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGCCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGATTTGGTCGGCAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGCGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtgc                                                                                                                                                                                               >  1:352025/1‑251 (MQ=255)                                                         gCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTtgct                                                                                                                                                                                              <  2:550409/247‑1 (MQ=255)                                                           tACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGCAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGACAAAGTGATGATGACATTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGCCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAt                                                                                                                                                                                        >  2:75150/1‑251 (MQ=255)                                                              tGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGCTTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGt                                                                                                                                                                                     >  2:517277/1‑251 (MQ=255)                                                               gCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGAGGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGt                                                                                                                                                                                     <  1:63341/250‑1 (MQ=255)                                                                       aTCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACTCGCAACGTTGCCGCACAACGCCTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGCATGCTCGGAGATACAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACg                                                                                                                                                                            >  1:567858/1‑251 (MQ=255)                                                                          aTGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCTTAGGTCATATTCAGGTGAtgtg                                                                                                                                                                                                           <  1:162998/217‑1 (MQ=255)                                                                          aTGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGAtgtg                                                                                                                                                                                                           >  2:162998/1‑217 (MQ=255)                                                                            ggCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGa                                                                                                                                                                       <  1:333459/251‑1 (MQ=255)                                                                              cATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAAGGCTCGGAGAAAAAGGCTTAGGTCATATTCAGGTGATTTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAag                                                                                                                                                                     <  1:317648/251‑1 (MQ=255)                                                                                    ccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTg                                                                                                                                                               >  1:94554/1‑251 (MQ=255)                                                                                      aGCGTTATGCAGAGGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTGCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGcc                                                                                                                                                             <  2:318541/251‑1 (MQ=255)                                                                                       gCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCg                                                                                                                                                            <  2:525740/251‑1 (MQ=255)                                                                                                         cGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCt                                                                                                                                          <  1:135023/251‑1 (MQ=255)                                                                                                          gatgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCCTTGGGGATGGAACCGGAAAAAGCGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCtt                                                                                                                                         >  1:237931/1‑251 (MQ=255)                                                                                                          gatgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCtt                                                                                                                                         >  1:297699/1‑251 (MQ=255)                                                                                                           atgatTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTc                                                                                                                                        <  1:200967/251‑1 (MQ=255)                                                                                                                    cGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGcc                                                                                                                               >  1:524432/1‑251 (MQ=255)                                                                                                                    cGCAACGCTGCCGCACAGCGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCTGAAAAAGTGATGATGACCTTTCAGTCGGGCTTTGGTCGGGAACCCAGAC‑GACTGCCTTATACCGACGAATCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTAATATTCAGGTGATGTTCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGCGAGCACTTCCTCGATGcc                                                                                                                               <  2:24684/251‑1 (MQ=255)                                                                                                                     gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcg                                                                                                                              >  1:274566/1‑251 (MQ=255)                                                                                                                       gaCGTTGCCGCACACCGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATTACCTTTCAGTCGCGCTTTGGTCGGGAAACCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGACTGCTCGGAGAAAAAGGCGTAGGCCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGACGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggg                                                                                                                            >  2:418689/2‑250 (MQ=255)                                                                                                                       aaCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCTTGGCTGA‑TGCCTTATACCGATGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggc                                                                                                                            <  2:500740/251‑1 (MQ=255)                                                                                                                        aCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGTCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTGTTCCTCGGTGCcggcg                                                                                                                           <  2:277744/251‑1 (MQ=255)                                                                                                                        aCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcggcg                                                                                                                           <  1:491191/251‑1 (MQ=255)                                                                                                                          gTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCggg                                                                                                                         >  1:174595/1‑251 (MQ=255)                                                                                                                             gCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaa                                                                                                                      <  1:299676/251‑1 (MQ=255)                                                                                                                                   caACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTTATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTTCCCGGGCGTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTTCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAatat                                                                                                                <  1:55413/251‑1 (MQ=255)                                                                                                                                   caACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGCGATGATGACCTTTCAGTAGCGCTTTGGTCGGGCACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGCAAAAGGCGTAGGACATATTCAGGCGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGGGCCGGCGGGaaaaaaaat                                                                                                                >  1:581501/1‑248 (MQ=255)                                                                                                                                      cGACTCGTGAACTGGCTGCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTGTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCGCGGTGCCGGCGGGAAAAAATATGaa                                                                                                             <  1:505402/251‑1 (MQ=255)                                                                                                                                               aaCTGCCTGCCGCATGGGTGATGGCACCGCAAAAAGGTATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGGATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCGTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGACTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCCGTGAAAAATACGAATATATTCCg                                                                                                    <  2:321458/251‑1 (MQ=255)                                                                                                                                                   ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCACAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAAGATATTTCggggg                                                                                                >  2:443582/1‑248 (MQ=255)                                                                                                                                                      ttCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTCCGGCGCtt                                                                                              >  1:322045/1‑250 (MQ=255)                                                                                                                                                        ccGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAGTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTGCCTCGGTGCCGGCGGGAAAAAATATGAATATATGCCGGCGCTTAAt                                                                                           <  2:274566/251‑1 (MQ=255)                                                                                                                                                              tGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTCTACCGACGAAACGCTGAGAATGCTCGGAGAAACCGGCGTAGGTCATCTTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTCAT‑Gacacc                                                                                    >  1:348251/1‑250 (MQ=255)                                                                                                                                                                        aCCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACata                                                                          <  1:371763/251‑1 (MQ=255)                                                                                                                                                                           ggAAAAAGTGATGATGACCGTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTGAAT‑GCCACGCCGGAACATATCg                                                                       <  2:265936/251‑1 (MQ=255)                                                                                                                                                                                  gTGATGATGGCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTAGACCTACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCATGTGATGTGCCCGGGCTTTGCGGCCGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAGATTCCGGCGCGTAAT‑GCCTCGCCGGAACATATCGAAatgat                                                                <  2:567858/251‑1 (MQ=255)                                                                                                                                                                                   ggaagatgaCCTTGCAGGCGCGCTTTGGGCGGGAACCCTGGCTGA‑GGCCTTAGACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCTGAGGTCAGATTCAGGTGATGGGCCCGGGGTTTGCTGCGGATTGTCTGGAGACGTTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAA‑GCCACGCCGGACCATATCGAAatgatg                                                               <  2:352025/247‑1 (MQ=255)                                                                                                                                                                                     atgatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAAAGATggg                                                             >  1:170194/1‑250 (MQ=255)                                                                                                                                                                                        atgaCCTTTCAGTCGCGCTTTGGTCGTGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCGGCGTATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAACCCGTGAGTTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCTGAACATATCGAAATGATGGCTa                                                           <  1:242730/250‑1 (MQ=255)                                                                                                                                                                                              tttCAGTCGAGCCTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATGTTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTATGGAGACGCTGGACGAGATTGCCGAGCAACATCGTGAGGTCTTCCTCGGTGCCGGGGGGAAGAAATATGGATATATTCCGGGGCTCGAT‑TCCACGGCGGAACAAATCGCAatgatagcgcattttgt                                                    >  2:350058/1‑239 (MQ=255)                                                                                                                                                                                                 cAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTAATCttgttg                                                  <  1:79425/250‑1 (MQ=255)                                                                                                                                                                                                 cAGTCGCGCTTAGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGCTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGCGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAATTATTTCGGCGCCTAATTGCC‑CGCCGGAGCTTATCGaaaagtacggtaagcttgtttt                                                 >  2:106069/1‑232 (MQ=255)                                                                                                                                                                                                                cGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCGTAAT‑GCCACGCCAGAACATAACGAAATGAGGGCTAATCTTGTTGCCGCGTATTGCTAAAg                                  <  1:72205/251‑1 (MQ=255)                                                                                                                                                                                                                 gggAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTTATCTTGTTGCCGCGTATCGCTAAAGc                                 >  1:75829/1‑251 (MQ=255)                                                                                                                                                                                                                        cTGGCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGGGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAGGAACATATTTCGGCGCTTCAT‑GCCCCGCCGGGCAATATATAAATGATGGCTAATATTTTTGcgcagtatcgcgcaagccgcgcgg                          >  2:470936/1‑228 (MQ=255)                                                                                                                                                                                                                           gCTGA‑TGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAAAATATTTCGGCGCTTATT‑GCCCAGCCGGAACATATCGCACAGATGGCTAATATTGTTGGcgcggatagctcacgctgtgcgggt                        >  1:26835/1‑229 (MQ=255)                                                                                                                                                                                                                                  gCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATGCCGGCGCTTAAT‑GCCTCGCCGGAACATATCGAAATGATGGCTAGTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACt                 <  1:37343/251‑1 (MQ=255)                                                                                                                                                                                                                                            gacgaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTc       <  2:420730/251‑1 (MQ=255)                                                                                                                                                                                                                                            gacgaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTc       <  2:181977/251‑1 (MQ=255)                                                                                                                                                                                                                                             acgaAACGCTGAGAATGCTCGGCGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAAGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTTCGGCGCTTTAA‑GCCCCGCCGGGACATATCGAAAAGATGGCTAATCTTTTTTCCGCGGATCGCTAAAGCTGAGCGGGAAAGAACAGAGCGCCggcg      >  1:476519/1‑248 (MQ=255)                                                                                                                                                                                                                                                 aaCGCTGAGAATGCTCGGAGAAAAAGTCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCGTAAG‑GCCACGCCGGAACATATCGAAATTATGGCTAATCTTGTTGCCGCGTATCGCTAAATCGGAGCGGTAAAGAACTGAGCGCGGCCGCGaa  <  2:567228/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aaCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGaa  <  2:451376/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGA‑TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAT‑GCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAA  >  NC_000913/498562‑499051

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