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%	184.1 / NA	51	D43N (GAC→AAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/27);  total (24/27)

CACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAC  >  NC_000913/389361‑389841                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 cacaGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttct                                                                                                                                                                                                                                        >  1:352822/1‑251 (MQ=255)       aGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttcttc                                                                                                                                                                                                                                   >  1:192642/1‑250 (MQ=255)          aGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaac                                                                                                                                                                                                                                >  1:155516/1‑250 (MQ=255)              gtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                              >  1:58463/1‑248 (MQ=255)              gtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAg                                                                                                                                                                                                                            <  1:186892/250‑1 (MQ=255)                  gACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATc                                                                                                                                                                                                                         <  1:85705/249‑1 (MQ=255)                          gATCATTTCTGGCAAGGTATGCTTGCAGATGCGCGACATACGCGCCACCAGGCCAACTTCCCGCCAGGCATCGCTGCCGGTTTCATAGGTATGATGAGCGATGCCATAAGTCATCACGGACCGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGCATATCTGGAATGCGAATCACGCCTGGCATGGATTCAACGGCTTTGTAGTTGTCTATTTATTCTTCATCAAAGCTCGGcaccactc                                                                                                                                                                                                               >  2:108996/1‑249 (MQ=255)                             cATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGt                                                                                                                                                                                                            <  1:385845/251‑1 (MQ=255)                             cATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGt                                                                                                                                                                                                            >  1:175479/1‑251 (MQ=255)                                atatGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCAGCTTACCGCGAGGCATCGCTGCCGGTTTCATCGGTATGGAGAGAGATGCCAACAGTCATCACGGAACGAATACCGGCGTTGGCGCTGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGACTGCGCAACACGCCTGGCTTGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGt                                                                                                                                                                                                         >  2:462/4‑251 (MQ=255)                                      cACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGAATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACTGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGc                                                                                                                                                                                                   <  1:93513/251‑1 (MQ=255)                                        cGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGCTGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGGTGTCACTTTCTTCTTCAACAAAGATCGGCAACACCAGGACGTTAAGGCt                                                                                                                                                                                                  >  2:298301/1‑250 (MQ=255)                                        cGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTa                                                                                                                                                                                                 <  1:427599/251‑1 (MQ=255)                                         ggTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTaa                                                                                                                                                                                                <  1:11774/251‑1 (MQ=255)                                                 tGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGtgttgt                                                                                                                                                                                         <  1:352765/250‑1 (MQ=255)                                                    aGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTctct                                                                                                                                                                                     >  1:118864/1‑251 (MQ=255)                                                    aGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTctct                                                                                                                                                                                     <  2:192642/251‑1 (MQ=255)                                                        gcgcGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCaa                                                                                                                                                                                 <  1:37683/251‑1 (MQ=255)                                                                          cagtcactCTTCCCGACAAGCATCGCTGCCGGTTTAATCTGCATGGCTAGAGATGCCACAAGTCAACACGGAACGAATACCAGCGTTGGCGATACGTTCAATTTCGCGTGCCAGATGCTTCTATGGAATGCGCATCACGCCTGGCATGGCATCAACGGCTTTGCAGTTGTCAATTTCTTCTTCAACAAAGATCGGCACCACCAGGTCGTTAAGGCTAAACGTTGTATCTTCAAACATAGCgcgcagcgc                                                                                                                                                                 >  2:121798/8‑249 (MQ=255)                                                                                  cctCCCGCCAGGCATCGCTGCCGGTTTAGTCGGTATGGGGAGAGAGGCCAAAAGTCAGCACGGAACGAATACCGGCGTTGGCGATGCGGTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCTCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAAGTTCGTCTTCAACAAAGATCGGCATCACCAGGTCGTTAAGGCTAAGGTTTGTCTCGTCAAACTTAGCGAGCAGCGCAGGAGAttt                                                                                                                                                        <  2:356488/248‑1 (MQ=255)                                                                                    tCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTg                                                                                                                                                       <  1:161524/249‑1 (MQ=255)                                                                                           aGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcga                                                                                                                                              <  1:413100/251‑1 (MQ=255)                                                                                                 cGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCAGCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAgg                                                                                                                                         <  1:298301/250‑1 (MQ=255)                                                                                                  gCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGc                                                                                                                                       <  1:412556/251‑1 (MQ=255)                                                                                                                  gTATGGTGAGAGATGCCAAAAGGCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGt                                                                                                                       <  2:227878/251‑1 (MQ=255)                                                                                                                          agagaTGCCAAAAGTCATCACGGAACGAATACAGGCGTTGGCGATGCGTTCAATTTAGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGACGTTACGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCCGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTCGGATTAAGACTGTTATGGtat                                                                                                               >  1:1221/1‑249 (MQ=255)                                                                                                                              aTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCt                                                                                                           >  1:219075/1‑251 (MQ=255)                                                                                                                                               ggAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCg                                                                                           >  1:207037/1‑250 (MQ=255)                                                                                                                                                gAACAAATCCCGGCGTTGGCGATGCGTTCAATTTCTCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTTAATTTCTTCTTCAACAAAGATCGTCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCATGCGCCGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGTAATCGAg                                                                                         <  2:329171/251‑1 (MQ=255)                                                                                                                                                                           tCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTa                                                               <  2:117156/250‑1 (MQ=255)                                                                                                                                                                                  cgcgTGCCAGATGGTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTGAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCTTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                                        <  1:122886/250‑1 (MQ=255)                                                                                                                                                                                       gCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCCTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTAc                                                   <  2:213580/250‑1 (MQ=255)                                                                                                                                                                                        ccAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                                                  >  2:92473/1‑250 (MQ=255)                                                                                                                                                                                          aGATGTTTCTCTGGAATGCGCATCACGCCTTGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                                               <  2:296120/251‑1 (MQ=255)                                                                                                                                                                                          aGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAGCGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGTAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                                               <  2:73689/251‑1 (MQ=255)                                                                                                                                                                                           gctGTTTCTCTGGAATGGGAATCAAGCCTGGCCTGGCTTCAACGCCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGCGATTTGCGCAGGCGACGAGGGAGTTGGATTAAGTCTGTCATGGTCTGCCTGCTGTTTGTGGAATCGAGGGGCATAGTATATCTGCAGCAGGCTAGGGATGTTTTACGaa                                                >  2:147210/3‑249 (MQ=255)                                                                                                                                                                                           gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaa                                                >  1:406225/1‑249 (MQ=255)                                                                                                                                                                                           gATGTTTCGCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACGCCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCGGATGTTTTTGTAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGGTTTACGAAAg                                              <  2:134869/251‑1 (MQ=255)                                                                                                                                                                                                  ctctGGAATGCGCATCAAGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgc                                       >  2:407179/1‑251 (MQ=255)                                                                                                                                                                                                              cATCACGCCTGGCATGGCTACAACGGCCTTGTAGTTGTCAATTTCTTCCTCAACAAAGATCGGCAACACCCGGTCGTTCAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGAAGAGGGCGTTGGATTCAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGCTGTTGACa                             >  2:318586/1‑249 (MQ=255)                                                                                                                                                                                                               aTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGCGATTTGCGCAGGCGACGTGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGATGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACa                             >  2:180178/1‑248 (MQ=255)                                                                                                                                                                                                                           aTGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGTTGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGa                >  2:40097/1‑249 (MQ=255)                                                                                                                                                                                                                               cTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAagag            >  1:209090/1‑249 (MQ=255)                                                                                                                                                                                                                               cTTCAACGGCTTTGTAGTTGTCAATTTCTACTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAagag            >  1:395300/1‑249 (MQ=255)                                                                                                                                                                                                                                     cGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAgg    <  2:266931/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAgg    >  1:192783/1‑250 (MQ=255)                                                                                                                                                                                                                                      ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGa   <  1:367113/251‑1 (MQ=255)                                                                                                                                                                                                                                       gCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAc  <  1:360859/251‑1 (MQ=255)                                                                                                                                                                                                                                       gCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAc  <  1:81565/251‑1 (MQ=255)                                                                                                                                                                                                                                          ttGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGt                                             >  1:444235/1‑205 (MQ=255)                                                                                                                                                                                                                                          ttGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGt                                             <  2:444235/205‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 CACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAC  >  NC_000913/389361‑389841

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