Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,116	C→T	P21L (CCT→CTT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,116	0	C	T	100.0%	201.2 / NA	55	P21L (CCT→CTT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/25);  total (30/25)

CGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGG  >  NC_000913/497887‑498349                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          cGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCt                                                                                                                                                                                                                      <  2:389654/251‑1 (MQ=255)       aGCAAAATTATCAGCAATATCAGGCCGGATATTCATTCGGCCTTTTACAAAAAACTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAAACTGGGTACGCCCGATGGCCCCACACTTGAAGCGGGAAAACGCTATCTGAAAc                                                                                                                                                                                                                 >  2:168935/1‑250 (MQ=255)        gCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTCATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACa                                                                                                                                                                                                                >  2:195823/1‑250 (MQ=255)         cAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAt                                                                                                                                                                                                              >  2:527239/1‑251 (MQ=255)          aaaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAtt                                                                                                                                                                                                             >  1:409294/1‑251 (MQ=255)           aaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAttt                                                                                                                                                                                                            <  1:137644/251‑1 (MQ=255)           aaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAttt                                                                                                                                                                                                            >  2:109859/1‑251 (MQ=255)                    gCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAc                                                                                                                                                                                                   <  1:437685/251‑1 (MQ=255)                    ccaatccaaGGCCGGATATTCATCCGGCCTTT‑ACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAACACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACa                                                                                                                                                                                                  <  2:35879/243‑1 (MQ=255)                           cAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCGGCGCTATGCCTTTTTCTGATTTGACCGCTCACAGCAATTAGTTCTTCTTCCTCGCTTTTCCGCTACAATAATCAACAAGTTGAATCGATAAGAGTCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcg                                                                                                                                                                                              <  2:179548/249‑1 (MQ=255)                            aGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcgcgg                                                                                                                                                                                           >  1:26387/1‑250 (MQ=255)                            aGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGt                                                                                                                                                                                           <  1:414985/251‑1 (MQ=255)                                       tcattcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCt                                                                                                                                                                                >  1:146441/1‑251 (MQ=255)                                         attcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTcc                                                                                                                                                                              >  2:50767/1‑250 (MQ=255)                                           tcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCAc                                                                                                                                                                             <  1:487189/250‑1 (MQ=255)                                              gCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgtt                                                                                                                                                                         <  2:407274/251‑1 (MQ=255)                                              gCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgtt                                                                                                                                                                         >  1:546066/1‑251 (MQ=255)                                                       aaaaaCCTCAATCTGCGCTAGTCCTTCTTCTGATTTGGCCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAATCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCTCACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggtg                                                                                                                                                                <  1:309691/251‑1 (MQ=255)                                                        aaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggtgg                                                                                                                                                               >  1:544002/1‑251 (MQ=255)                                                        aaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggtgg                                                                                                                                                               >  1:17716/1‑251 (MQ=255)                                                           aCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTCAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGcc                                                                                                                                                             >  2:386415/1‑250 (MQ=255)                                                                aaTCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtgct                                                                                                                                                       <  2:26387/251‑1 (MQ=255)                                                                aaTCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtgct                                                                                                                                                       >  1:254071/1‑251 (MQ=255)                                                                          aTTCCTTTGTCAGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACGCTAGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGGTGATACCTCACGGTTGTTATGGTGGCCATtgctgctaggcgt                                                                                                                                              <  2:71424/250‑8 (MQ=255)                                                                                   tCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGcc                                                                                                                                     <  1:204235/250‑1 (MQ=255)                                                                                      gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCt                                                                                                                                  <  2:146441/250‑1 (MQ=255)                                                                                         ttGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgcgc                                                                                                                              >  2:55339/1‑251 (MQ=255)                                                                                          tGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgcgc                                                                                                                              <  1:291326/250‑1 (MQ=255)                                                                                           gACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTTCCGCTGCGCTc                                                                                                                            >  2:198103/1‑251 (MQ=255)                                                                                           gACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTc                                                                                                                            <  1:354271/251‑1 (MQ=255)                                                                                                                   nttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGCGGCGGTAATGCGTAAGACTAAAACCGGTAAACTGCTGGCAAAACTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGGTACCTCACGGTTGTTATGGTGGCCATTTCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGccttatgc                                                                                                    >  2:255532/2‑244 (MQ=255)                                                                                                                      ccTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCctc                                                                                                 <  1:55339/251‑1 (MQ=255)                                                                                                                      ccTCAATTTTCCGCTACAATTATCAACCAGTTGAATCGATAAGAGCCGGTACTGTGTCAGACTAAAACCGGTATCCTGCTGGCCAACCTGGGCACGCTCGATGCCCCCACACTTGAAGCCGCAAAACGCTATCTGAAACAATTTTTAAGCGATAGACGCGTCGGTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGATTTTGACGCCGCGCTCGCCGCGTGTGGCGAAGcgggttgccct                                                                                                 >  2:338093/1‑241 (MQ=255)                                                                                                                         cACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACCGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCtctgt                                                                                              >  1:302695/1‑251 (MQ=255)                                                                                                                          aCTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCtctgtc                                                                                             <  1:180973/251‑1 (MQ=255)                                                                                                                                 cGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGGACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGGCtggatgg                                                                                      >  2:216059/1‑251 (MQ=255)                                                                                                                                       aaTTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                <  1:50767/251‑1 (MQ=255)                                                                                                                                        aTTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGTAggggg                                                                               >  2:157743/1‑248 (MQ=255)                                                                                                                                              aacaaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCTCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgcc                                                                         <  2:560802/251‑1 (MQ=255)                                                                                                                                              aacaaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgcn                                                                         <  2:254071/251‑2 (MQ=255)                                                                                                                                                     tGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATg                                                                  >  2:326381/1‑251 (MQ=255)                                                                                                                                                           gATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGAGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                            >  1:498678/1‑251 (MQ=255)                                                                                                                                                            aTAAGAGGCGGAAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCCCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                            >  2:483203/1‑250 (MQ=255)                                                                                                                                                              aaGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgc                                                         >  1:190493/1‑251 (MQ=255)                                                                                                                                                                gagGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccg                                                       <  1:198103/251‑1 (MQ=255)                                                                                                                                                                    cGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgcc                                                     >  1:203939/1‑249 (MQ=255)                                                                                                                                                                        aaTGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAc                                               <  1:533216/251‑1 (MQ=255)                                                                                                                                                                         aTGCGTCAGACTAAAACCGGTATCCTGCTGGCAACCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGGGTGGCGAAGATGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGccac                                               >  2:368963/1‑247 (MQ=255)                                                                                                                                                                                     aaaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTg                                    <  1:67240/249‑1 (MQ=255)                                                                                                                                                                                         ccGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcaca                               >  1:10059/1‑250 (MQ=255)                                                                                                                                                                                           ggTATCCTGCTGGTAAACCTGGGTTCGCCCGATGCCCCCTCACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCTTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCGGATGGTTTACAGCCGCCAGCATCAGCAGTCGCTGGCACAACg                            <  2:469264/251‑1 (MQ=255)                                                                                                                                                                                                   gctgGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCgg                    <  1:323111/251‑1 (MQ=255)                                                                                                                                                                                                         aaaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGATCGCCGAGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACATCAGGAGATGGCACAACGTTTACCGgagaagc              >  2:162406/1‑248 (MQ=255)                                                                                                                                                                                                            ccTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGAATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCg           >  2:3688/1‑251 (MQ=255)                                                                                                                                                                                                                     gCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTgg  >  1:379411/1‑251 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          CGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGG  >  NC_000913/497887‑498349

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