Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	389,601	C→T	D43N (GAC→AAC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	389,601	0	C	T	100.0%	234.2 / NA	67	D43N (GAC→AAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/39);  total (28/39)

TGTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAC  >  NC_000913/389355‑389841                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 ggtATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCa                                                                                                                                                                                                                                              <  1:47957/250‑1 (MQ=255)  gtATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCCCGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCaa                                                                                                                                                                                                                                             <  2:57373/251‑1 (MQ=255)  gtATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCAGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGATGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGCTGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACAGCGTTGAAGTTGTCa                                                                                                                                                                                                                                              >  1:231411/1‑250 (MQ=255)      tCACAGAATGAGGTGTCTGACATAACGATCATTTCTGGCACGGTCGGCTTGCAGGTGCGCCACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCGGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGCCTTTGTAGTTGTCAAttt                                                                                                                                                                                                                                          <  2:490556/250‑1 (MQ=255)       cacaGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttc                                                                                                                                                                                                                                         <  1:529862/250‑1 (MQ=255)        acaGAAGCAGGTGGCGGACATAACGATCATTTCTGGCACGGGCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttctt                                                                                                                                                                                                                                       <  1:240235/251‑1 (MQ=255)         caGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGAATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttc                                                                                                                                                                                                                                      >  2:506020/1‑251 (MQ=255)            aaGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttcttc                                                                                                                                                                                                                                   >  1:50661/1‑251 (MQ=255)            aaGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttcttc                                                                                                                                                                                                                                   >  1:55615/1‑251 (MQ=255)             aGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACAGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCCACGGCTTTGTAGTTGTCAATttcttcttc                                                                                                                                                                                                                                   >  2:144297/1‑250 (MQ=255)                  gtgtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                              >  1:102421/1‑250 (MQ=255)                   ggtCTGACATAACGAACATTTCTGGCACGGTCTGCTGGCAGATGCGCGACAT‑CCGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTGCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATCCCGGCGTTGGCGATGCGTTCAATGTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCACCGGCTTTGTAGTTGTCAATTTCTGCTTCAACAAAg                                                                                                                                                                                                                            <  2:198320/250‑1 (MQ=255)                    gtcggccATAACGAGCATTTCTGTCACGGGTTGCTTGCAGAGTCGCGACAT‑ATGCGCGACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTAGTCGGTATGGTGAGAGATGCCAAAAGGTATCACGGAACGAATTCCGGCGTTGGCTATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTCGTAGTTGTCCATTTCTTCTTCAACaaaa                                                                                                                                                                                                                            <  2:140448/244‑2 (MQ=255)                              aCGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGACACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGATTCATCGGTATGGTGAGAGATGCCAAACGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCACTTTCTTCTTCAACAAAGATCGGCAacac                                                                                                                                                                                                                 >  2:475203/1‑251 (MQ=255)                                   cATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAgg                                                                                                                                                                                                             <  1:164194/250‑1 (MQ=255)                                       tCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGtt                                                                                                                                                                                                        >  2:58525/1‑251 (MQ=255)                                           gCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAgg                                                                                                                                                                                                    <  1:77166/251‑1 (MQ=255)                                            ccgggtcggtttGCCGATGCGCG‑CATTACGCGCCACCAGTCCATCGTCCCGCCACGCATCGCTGCCGGTTGTATCGGGATGGTTAGAGATGCCAAAAGTCATCCCGGAACGATTACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCATCGGCTTTGTAGTTGTCACTTTCTTCTTCAACCAAGATCGTCAACACCAGGTCGTTAAgtc                                                                                                                                                                                                   <  2:402096/241‑3 (MQ=255)                                             gcGGTCTGCTTGCAGATGAGCGGCAT‑TCGCGCAGCCAGTTCCTCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATTCCAAAAGTCATCACGGAACGAAGACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCt                                                                                                                                                                                                  <  2:208151/250‑1 (MQ=255)                                                gTCTGCTTGCAGCTGAGCGACAT‑ACGCGCCACCAGTCCCTCTTCCAGCCAGACATCGCTGACGGTTTCATCGGCATGGTTAGAGATGCCAACAGTCATCGCGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCCTGGCTTCAACGGCTTTGTAGGTGTCAATTTCTTATTCAACAAAGATCGGCAACACCAGGTCGtttcggctacg                                                                                                                                                                                               >  2:573919/1‑242 (MQ=255)                                                gTCTGCTTGCAGATGCGCGACAT‑AAGCGCCACCAGTCCATCTTCCCGCCAGGAATCGCTGCCGGATTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCAAGATGTTTCTCTGGAATGCGCATCACGCAGGGCATGGCTTCACCGGCATTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGACAACACCAGGACGTTACGGCTa                                                                                                                                                                                                 >  1:205116/1‑249 (MQ=255)                                                              gcgcGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCa                                                                                                                                                                                  <  1:100740/250‑1 (MQ=255)                                                                  gACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAt                                                                                                                                                                             <  1:27779/251‑1 (MQ=255)                                                                  cgaccgccgcgCCACCAGTCCATCTACCCGCCAGGCATCGCTGCCGTTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGATCGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTTGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCaacca                                                                                                                                                                              <  2:410169/244‑4 (MQ=255)                                                                               accaGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCgcgcagcgc                                                                                                                                                                 <  1:58525/251‑1 (MQ=255)                                                                                   gTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGATAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAgg                                                                                                                                                              >  1:414191/1‑250 (MQ=255)                                                                                     ccATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGgaga                                                                                                                                                           <  1:517922/251‑1 (MQ=255)                                                                                            cccGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcg                                                                                                                                                     <  1:154348/250‑1 (MQ=255)                                                                                                  aGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcga                                                                                                                                              <  2:415304/251‑1 (MQ=255)                                                                                                         gCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGc                                                                                                                                       <  1:130808/251‑1 (MQ=255)                                                                                                            gCCGGGTTCATCGGTATGGTGAGAGATGCCAAAAGGCATCACGGAACGAATACCGGCGGTGGCGATGCGTTCAATTTCGCGTGCCAGATGGTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTGTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCTCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGtt                                                                                                                                    <  2:505214/251‑1 (MQ=255)                                                                                                                  ttCAGCGGTATGGTGAGAGATGCCAAAAGTCAGCATGGAACGAATACCGGCGTTGGCGATGCGTTTAATTTCTCGTGCCAGATTTTTCTCTGGAATGCGCATCACGCCTTGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCATCTTCAATAAAGATCGGCAACACCAGGTCGTTAATGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGGTGGAtt                                                                                                                               <  1:538731/250‑1 (MQ=255)                                                                                                                   tCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTa                                                                                                                              <  1:327814/250‑1 (MQ=255)                                                                                                                          tATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTc                                                                                                                      >  1:374125/1‑251 (MQ=255)                                                                                                                              gTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATgg                                                                                                                  >  1:534445/1‑251 (MQ=255)                                                                                                                                 ggagaTTCCAGAAGGCATCACTGGACGACTACCGGCGTTGGCGATGCGTTCAATTTCGTGTGCCAGATTTTTCTCGGTAATGCGCATCACGCCTGGCATGGCTTCTACGGCTTGGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGGGTTGTCGCTTCAACCATAGCGCGCAGCGCAGGAGATTTGCGGAGGCGGCGAGGGCGTTTGAGGAAGTCTGTCATTGTCt                                                                                                               <  2:335138/250‑1 (MQ=255)                                                                                                                                         cAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATg                                                                                                       >  1:412688/1‑251 (MQ=255)                                                                                                                                          aaaaGTCATCACGGAACGAATAACGGCGTTGGCGATGAGTTCAATTTCTAGTGCAAGATGTTTCTCTGGCATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTAAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAACGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGt                                                                                                      >  2:106665/1‑251 (MQ=255)                                                                                                                                           aaaGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGtt                                                                                                     <  1:480031/251‑1 (MQ=255)                                                                                                                                            aaGTCATCACGGCACGAATCCCGGCGTTGGCGTTGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGTTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGttt                                                                                                    <  2:13657/251‑1 (MQ=255)                                                                                                                                             aGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCACTTTCTTCTTCAACAAAGATCGGCAACACCCGGTCGTTAAGGCTAGGTGTTGTCTCTTCAAACATAGCGCGCCGCGCAGGAGATTTGCGCAGGCGTCGGTGGCGTTGGATTAAGTCTGTCCTGGTCTGCCTGATGTTtg                                                                                                   >  2:16079/1‑251 (MQ=255)                                                                                                                                                  tcaCGAAACGAATAACGGCGTTGGCGATGCGGTCAATTTCGCGAGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGAAGGAGATTTGCGCCGGCGACGCGGGCGTTTGATTAAGTCTGTCCTGGTCTGCCCGCTGTTTGTggcc                                                                                              >  2:86945/1‑249 (MQ=255)                                                                                                                                                    aCGGAACGAATTCCGGCGTTGGCGATGCGTTCAATGTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGGCATGGTCTGCCTGATGTTTGTGGAATc                                                                                            <  1:438986/251‑1 (MQ=255)                                                                                                                                                     cGGAACGAATGCCGTCGTTGGCGATGCGCTCTAGTTCGCGTGCCAGATGTTTCTCTGGCATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATAGGCAACACCCGGTCGTTAAGGCTAAGTGTTTTCTCTTCAAACATAGCGCGGAGCGCAGGAGATTTGCGTAGGCGACGAGGGCGTTGGATTAAGTCTGTCGTGGTCTGCCTGATGTTTGTGGAATc                                                                                            <  2:179552/250‑1 (MQ=255)                                                                                                                                                             aTACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                                   >  1:404805/1‑251 (MQ=255)                                                                                                                                                               aCCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAg                                                                                 <  2:55615/251‑1 (MQ=255)                                                                                                                                                                 cGGCGTTGGCGATTCGGTCAATTTCGCGGGCCAGATGTATCGCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTATGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAGCAGAGCGCGCAGCGCATGAGATTTGCGCAGGCGACGAGGGCGTTGGATGAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGta                                                                               <  2:8713/251‑1 (MQ=255)                                                                                                                                                                  ggCGTTTGCGATGCGTTCATGTGCGCGTGCGAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATTGCTTCAAGGGCGTTGTAGTTGTCAATGTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGAGCAGCGCAGGAGATTTGCGCGGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGGCTGATGTTTTTGGAATCGAGGGGCATAGtag                                                                              <  2:2040/251‑2 (MQ=255)                                                                                                                                                                   gCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGtata                                                                             <  1:399142/251‑1 (MQ=255)                                                                                                                                                                           gATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGTCGAGGGCTTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGGATAGTATACCTGAAGc                                                                     >  2:179034/1‑251 (MQ=255)                                                                                                                                                                                gTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCCATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTc                                                                                                                      >  2:27403/1‑197 (MQ=255)                                                                                                                                                                                gTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTc                                                                                                                      <  1:27403/197‑1 (MQ=255)                                                                                                                                                                                    aaTTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAggg                                                            <  1:475203/251‑1 (MQ=255)                                                                                                                                                                                        tCGCGTGCCAGATGTTTCTCTGGAATGCGCACCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATg                                                         <  1:399254/250‑1 (MQ=255)                                                                                                                                                                                        gcgcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCGGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGGTTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGAGGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                                        <  2:567050/250‑1 (MQ=255)                                                                                                                                                                                              gCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                                                  >  1:132904/1‑251 (MQ=255)                                                                                                                                                                                                 aGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaa                                                >  1:414489/1‑250 (MQ=255)                                                                                                                                                                                                 aGATGTTTCTATGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTTCGCAGGCGACGCGGGCGTTGGATTAAGCCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATAACTGAAGACGGGTAGGGATGTTTTACGaa                                                >  2:225935/1‑250 (MQ=255)                                                                                                                                                                                                       ttCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTc                                         >  1:483721/1‑251 (MQ=255)                                                                                                                                                                                                               aaTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGtt                                 <  1:449074/251‑1 (MQ=255)                                                                                                                                                                                                                 tGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTANGGATGTTTTACGAAAGTTGTCGCGATGTTg                                >  2:252758/1‑250 (MQ=255)                                                                                                                                                                                                                  gcgcATCACCCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAATATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGAc                              <  2:534445/251‑1 (MQ=255)                                                                                                                                                                                                                           gCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGGCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGCATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGACAGTTGTCGCGATGTTGACacgaag                        >  2:165539/1‑248 (MQ=255)                                                                                                                                                                                                                                gCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATgg                 >  1:518246/1‑250 (MQ=255)                                                                                                                                                                                                                                      cTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAagag            >  2:212916/1‑249 (MQ=255)                                                                                                                                                                                                                                         cAACGGCTTTGTAGTTGGCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGATAGGcc       <  2:102421/251‑1 (MQ=255)                                                                                                                                                                                                                                              ccTTTGTAGTTGTCAATTTCTTCGTCAACAAAGATCGGCAACACCAGGTCGTTAAGGTTACGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGTCGAGGGCGTTGGATTAAGTCTTTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAc  <  1:29305/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 TGTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACAT‑ACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAC  >  NC_000913/389355‑389841

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