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

TCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCT  >  NC_000913/389360‑389845                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      tCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCGGCGTGCAGATGCGCGACATACGCGCCCCCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATTCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGTCTTTGTAGTTGTCAAtt                                                                                                                                                                                                                                                <  2:261555/249‑1 (MQ=255)                   gACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                                   >  1:93345/1‑244 (MQ=255)                   gACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                                   >  1:417992/1‑244 (MQ=255)                    aCATAACGTTCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGc                                                                                                                                                                                                                          <  1:318639/251‑1 (MQ=255)                    aCATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAGCGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGGGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGc                                                                                                                                                                                                                          >  2:17346/1‑251 (MQ=255)                        aaCGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAaca                                                                                                                                                                                                                      <  1:207205/251‑1 (MQ=255)                            aTCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAg                                                                                                                                                                                                                  <  1:364265/251‑1 (MQ=255)                               aTTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCCGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTc                                                                                                                                                                                                               >  2:190925/1‑251 (MQ=255)                                tttCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCAGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAAAGGCTTTGTAGTTGTCAATTTCTTCTTCT‑CAAACGATCGGCAACACCAggacg                                                                                                                                                                                                              >  2:318581/1‑248 (MQ=255)                                 ttCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCg                                                                                                                                                                                                              <  1:166378/250‑1 (MQ=255)                                 ttCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCg                                                                                                                                                                                                              <  1:161100/250‑1 (MQ=255)                                     ggCAGGGTCTGCTTGCAGTTGCACGACATACGCGCCTCCATGCCATCTTCCCGCCAGGCATCTCTGCCGGTTTCATCGGTATTGTGAGAGATGTCAAAATTCAGCACGGAACGAATACCGGCGTGGGCGATTCGTTCAATTTAGCGTGCCAGATTTTTCTCTTGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAg                                                                                                                                                                                                         <  2:260340/251‑1 (MQ=255)                                       cACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGATTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGc                                                                                                                                                                                                       >  1:282486/1‑251 (MQ=255)                                                cTTGCAGATGCGCGACATACGCGCCACCAGTCAATCTTCCCGCCAGACATCGCTGCCGGTTTCAACGGTAAGGTGAGAGATCCCAAAAGTCATCCCGGAACGAAAACCGGCGTTGGCGATGCGTTCACTTTCGCGTGCCAGATGTTTCTCTGGAATGCGACTCACGCCTGGCATGGCTTCCCCGCCTTTGTAGTTGTTCATTTCTTCTTCAACAAA‑GGTCGGCAAAACCAGGTCGTTAAGGCTacatgttg                                                                                                                                                                                              >  2:375585/1‑244 (MQ=255)                                                                  aCGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcg                                                                                                                                                                             <  1:446145/250‑1 (MQ=255)                                                                     cgcCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCCCGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgcgcc                                                                                                                                                                         >  2:85408/1‑250 (MQ=255)                                                                       ccaccaGTCCATCTGCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCGCTTCAAACATAGCgcgcagc                                                                                                                                                                       <  1:343131/251‑1 (MQ=255)                                                                        caccaGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCgcgcagcg                                                                                                                                                                      <  1:71717/251‑1 (MQ=255)                                                                            aGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTNTTGTCTCTTCAAACATAGCGCGCAGCGCAgg                                                                                                                                                                  >  2:208934/1‑251 (MQ=255)                                                                             gTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTATTCAACAAA‑GATCGGCAACACCAGGTCGTTACGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAggg                                                                                                                                                                 >  2:319877/1‑250 (MQ=255)                                                                                            aGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGCTTTGCGCAGGcg                                                                                                                                                   >  1:292341/1‑250 (MQ=255)                                                                                                   gcggCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcgc                                                                                                                                                        <  2:461401/235‑1 (MQ=255)                                                                                                   gCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcgc                                                                                                                                                        >  1:461401/1‑238 (MQ=255)                                                                                                      gCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAgt                                                                                                                                                                                                  >  1:389476/1‑193 (MQ=255)                                                                                                      gCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAgt                                                                                                                                                                                                  <  2:389476/193‑1 (MQ=255)                                                                                                          gTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGGTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTGCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGAt                                                                                                                                    <  2:27/251‑1 (MQ=255)                                                                                                            ttCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGAtt                                                                                                                                   >  1:363882/1‑250 (MQ=255)                                                                                                              cATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTa                                                                                                                                  >  1:432150/1‑249 (MQ=255)                                                                                                                 cGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCt                                                                                                                             >  1:289135/1‑251 (MQ=255)                                                                                                                      tGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCAt                                                                                                                        >  1:31721/1‑251 (MQ=255)                                                                                                                        gTGAGAGATGCCAAAAGTCATCACGGAACGAATTCCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCATGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATgg                                                                                                                      <  2:306287/251‑1 (MQ=255)                                                                                                                           agagaTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCt                                                                                                                   <  1:219849/251‑1 (MQ=255)                                                                                                                             agaTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCACGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAAAGCGCGCAGCGCAGGAGATTTGCGCAGGCAAAGAGGGCGTTGGATTAAGTCTGTCATGGTCTGc                                                                                                                 >  1:433477/1‑251 (MQ=255)                                                                                                                              gaTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGGCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCGCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGcc                                                                                                                <  2:405375/251‑1 (MQ=255)                                                                                                                                    aaaaGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGt                                                                                                          <  1:167162/251‑1 (MQ=255)                                                                                                                                                  aaCGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGTCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAgg                                                                                            <  2:63857/251‑1 (MQ=255)                                                                                                                                                    cGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAgggg                                                                                          <  1:385964/251‑1 (MQ=255)                                                                                                                                                          ccGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGt                                                                                    >  1:290446/1‑251 (MQ=255)                                                                                                                                                          ccGGCGTTGGCGATGCGGTCAAGGTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGt                                                                                    <  2:309754/251‑1 (MQ=255)                                                                                                                                                                 tGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAGAGCGCGCAGCGCAGGAGATTTGCGCAGGGGACTAGGGCGTTGGAGTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCt                                                                              <  1:48351/250‑1 (MQ=255)                                                                                                                                                                  ggCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTg                                                                             <  1:85408/250‑1 (MQ=255)                                                                                                                                                                              aaTTTCGCGTGCCAGATGCTTATCTGGAATGCGCATAACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGTCGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGCATCGAGTGGCATAGTATACCTGAAGCAGGGTAggg                                                                >  2:173562/1‑251 (MQ=255)                                                                                                                                                                                   cgcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGATTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAAAACCCGGTCGTTAACGCTAAGTGTTGTCTCTTCAAACATAGAGCGCAGCGCAGGCGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCTTGGTCTGCCTGATGTTTGTGGAATAGTGGGGCATAGTATACCTGAAGCAGGGTAGGGATGtt                                                           >  1:439365/1‑251 (MQ=255)                                                                                                                                                                                     cgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGtttt                                                         <  1:395239/251‑1 (MQ=255)                                                                                                                                                                                           aGATGTTTCTCTGGAATGCGCATCCCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                                                   <  2:432150/251‑1 (MQ=255)                                                                                                                                                                                                        gAATGCTCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAAGCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGt                                      <  1:56374/251‑1 (MQ=255)                                                                                                                                                                                                               cATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaag                               <  1:17346/251‑1 (MQ=255)                                                                                                                                                                                                                   ccGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAagaaga                           <  1:208934/250‑1 (MQ=255)                                                                                                                                                                                                                         ggCATGGCTTCCACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                                                   >  2:295923/1‑221 (MQ=255)                                                                                                                                                                                                                         ggCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                                                   <  1:295923/221‑1 (MQ=255)                                                                                                                                                                                                                            aTGGCTTCAACGGCGTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAag                  <  2:93345/251‑1 (MQ=255)                                                                                                                                                                                                                               gCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAagaga               <  1:190925/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGcccgg        >  1:195976/1‑247 (MQ=255)                                                                                                                                                                                                                                        gCTTTGTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAc      <  1:205654/251‑1 (MQ=255)                                                                                                                                                                                                                                             gTAGTTGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCt  <  1:169486/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      TCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAA‑GATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCT  >  NC_000913/389360‑389845

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