Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,405,925	C→T	I413I (ATC→ATT)	purA →	adenylosuccinate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,405,925	0	C	T	100.0%	227.4 / NA	68	I413I (ATC→ATT)	purA	adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/45);  total (23/45)

CGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATCATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCG  >  NC_000913/4405688‑4406169                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     cGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTc                                                                                                                                                                                                                                         <  2:1511/251‑1 (MQ=255)  gTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGCTGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTcc                                                                                                                                                                                                                                        >  1:352892/1‑251 (MQ=255)  gTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTcc                                                                                                                                                                                                                                        >  1:426279/1‑251 (MQ=255)    tCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCATGCGGCGCTGAACTAAATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATAGTATATATACCGGTCCgg                                                                                                                                                                                                                                      >  2:185499/1‑251 (MQ=255)       ggATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCgg                                                                                                                                                                                                                                      <  1:433756/248‑1 (MQ=255)       ggATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGAt                                                                                                                                                                                                                                    <  1:262866/250‑1 (MQ=255)       ggATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATc                                                                                                                                                                                                                                   <  1:648643/251‑1 (MQ=255)       ggATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATGTACGAAACCAGGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGAt                                                                                                                                                                                                                                    <  2:259268/250‑1 (MQ=255)        gATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCTTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCg                                                                                                                                                                                                                                  <  1:228215/251‑1 (MQ=255)        gATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCg                                                                                                                                                                                                                                  <  2:362263/251‑1 (MQ=255)          tGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGt                                                                                                                                                                                                                                 >  1:281625/1‑250 (MQ=255)            gtcTGATAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACCACCACTCCGCTGGCAGCTGGCGACTGGAAAGGTGTAGAGCCGATTTACGAACCCATGCCGGGCTGGTCTTAATCCCCCTTCGGCGTGAAAGATTGTAGCTGCCGGCGGCAGTCGGCGCTGAACTTTATCACGCGTATTGAAGAGCTGGCTGGTGTGTCGATCGAGATTATCTCTAcctctcccgctcttact                                                                                                                                                                                                                              <  2:376515/249‑16 (MQ=255)                gAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGaa                                                                                                                                                                                                                           >  1:375811/1‑250 (MQ=255)                    gagGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGACAGGTGTAGAGCCGATTTACGAAACCCTGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGAATGGTGTGCCGATCGATATTATATCTACCGGTCCGGATCGTACTGAAACCAt                                                                                                                                                                                                                      >  1:254507/1‑251 (MQ=255)                         tAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTc                                                                                                                                                                                                                 >  1:243487/1‑251 (MQ=255)                          aaaCTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCCCTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTcc                                                                                                                                                                                                                                        >  2:314454/1‑227 (MQ=255)                          aaaCTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTcc                                                                                                                                                                                                                                        <  1:314454/227‑1 (MQ=255)                           aaCTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTg                                                                                                                                                                                                               <  1:417679/251‑1 (MQ=255)                               ctGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGcgcg                                                                                                                                                                                                           >  1:46512/1‑251 (MQ=255)                                  cGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATAAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGAcc                                                                                                                                                                                                        >  2:556865/1‑251 (MQ=255)                                       cTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTc                                                                                                                                                                                                   >  2:115812/1‑251 (MQ=255)                                       cTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTc                                                                                                                                                                                                   >  1:205978/1‑251 (MQ=255)                                         tACCGGATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTc                                                                                                                                                                                                   <  2:573616/249‑1 (MQ=255)                                                   cGGATGGTCGCGAATTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTc                                                                                                                                                                                       <  1:556865/251‑1 (MQ=255)                                                    ggATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCACGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGCTCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCt                                                                                                                                                                                      >  2:18402/1‑251 (MQ=255)                                                               aaGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTgg                                                                                                                                                                           >  1:502886/1‑251 (MQ=255)                                                                aGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGc                                                                                                                                                                          <  1:180174/251‑1 (MQ=255)                                                                    actacaccTCCGCTGGCAGCTGACGACTCGAAAGGTGTAGCGGCGATTTACGAAGCCATGCCGGGCTGGACTGAATCCTCCTTCGGCGTGAAAGATCTTAGCGGCCTGCCGCAGGCGGCTCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTTATTCCATGATTCTGCGCGACCCGTTCGACTCGTAATTCTGGTACGCCTGGCAGat                                                                                                                                                                      <  2:1227/244‑1 (MQ=255)                                                                             ccGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCCCCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTtgc                                                                                                                                                               <  2:637508/249‑1 (MQ=255)                                                                                  ggCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGTTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTtgcctgcc                                                                                                                                                          <  1:464308/249‑1 (MQ=255)                                                                                  ggCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCgg                                                                                                                                                        <  2:223272/251‑1 (MQ=255)                                                                                          acgacTGGAAAGGTGTAGAGCCGATTTGCGAAACCATGCCGGGCTGGTCTGAATCCACCGTCGGCGTGAAAGATCGTAGCGGGCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAg                                                                                                                                                <  1:309495/251‑1 (MQ=255)                                                                                          acgacTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAg                                                                                                                                                <  2:399333/251‑1 (MQ=255)                                                                                             acTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAgtgt                                                                                                                                             <  1:115812/251‑1 (MQ=255)                                                                                                          agagCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGTCCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCtgt                                                                                                                                <  2:126176/251‑1 (MQ=255)                                                                                                            agCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGCCGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTgtgt                                                                                                                              >  2:501825/1‑251 (MQ=255)                                                                                                              ccGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTTCTGTGTCg                                                                                                                            >  2:500151/1‑251 (MQ=255)                                                                                                               cGATTTACGAATCCATGCTGGGCTGGTCTGCAGCCACCTACGGCGTGAAATATCGTAGCGGCCTGCAGCAGGCGGCGCTGAACTATATCAAGCGTATTGCAGAGCTGACTGGTGTTCCGATCGATAGTATCTCTTCCGGGCCGGATCGTACTGAAACCATGAGCCTGCGCGACCGGTTCGACGCGTAATTCTGGTACTCCTGGCAGGTAGTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCgg                                                                                                                           <  2:651338/251‑1 (MQ=255)                                                                                                                   ttACGAAACCATGCAGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTa                                                                                                                        >  2:446429/1‑250 (MQ=255)                                                                                                                    tacaaaaCCATGCCGGGCTGGTCTTAATCCCCCTTCGGCGTTAAAGATCGTAGCGGCCTGCCGCAGGTGGTGCTGAACTATATCAAGCGTATTGAAGAGCTGTCTGGTGTGCCGATCTATATTAGCTCTGCCGGTCCTGATCCTATTGCAACCATGATTCTTCTCGTCGCGCTCGGCGCGTAATTCTGGTTCGCTTGTTCGATATTTTTCCTGCCGGGCTAAGACTGTTATACATTGccgggtcttttaa                                                                                                                       <  2:278222/246‑13 (MQ=255)                                                                                                                    tACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTaa                                                                                                                       <  1:150346/250‑1 (MQ=255)                                                                                                                          cgcATGCCGGGCTGGTTTGAGTCCACCTCCGGTGTGAAAGATCGTAGTGGCCGGCCGCAGGCGGCGCTGAATTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGCGCATTATTGTCGCTACCGGTCCGGATCGTACTTAAACCAGGATTCTGCGCGTCCCGTTCGGGGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGTGTACAGTGTTATACATTGCTTGGTCGGGTAAGCCATTa                                                                                                                <  2:537807/249‑1 (MQ=255)                                                                                                                            cATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAATCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACg                                                                                                              <  2:91826/251‑1 (MQ=255)                                                                                                                               gCCGTGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTa                                                                                                           <  1:116667/251‑1 (MQ=255)                                                                                                                               gCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCt                                                                                                            <  1:472678/250‑1 (MQ=255)                                                                                                                                       ggTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTGCCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCGGTGTCGGGTAAGCCATTACGCTATCCGac                                                                                                     <  1:202356/249‑1 (MQ=255)                                                                                                                                             aaTCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGtt                                                                                             <  1:535418/251‑1 (MQ=255)                                                                                                                                             aaTCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGtt                                                                                             <  2:281625/251‑1 (MQ=255)                                                                                                                                                   ccTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATcc                                                                                       <  2:243487/251‑1 (MQ=255)                                                                                                                                                           gTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCtttt                                                                               <  2:193311/251‑1 (MQ=255)                                                                                                                                                               aaGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTttcc                                                                           <  1:399215/251‑1 (MQ=255)                                                                                                                                                                 gATCGTAGCGGACTGCCGCATGCGGGGCTGAACTATTCCAAGCGTATGGAAGAGCTGACTGGTGTGCCGATCGATAATATCTCTACCGGTCCGGATAGAACTCAACCCATGATTCTGCGCGACCCGTTCTACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGTCGGGCGAACAGTGGCATGCATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACTCAGTGGTAAAGCCTCGCTTTTttcctt                                                                         <  2:346510/251‑1 (MQ=255)                                                                                                                                                                                      gcggcgCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAAtt                                                      >  2:118363/1‑249 (MQ=255)                                                                                                                                                                                             tGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTcttttacgccat                                                                                                                  <  1:363878/182‑12 (MQ=255)                                                                                                                                                                                             tGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCAt                                                                                                                  >  2:363878/1‑182 (MQ=255)                                                                                                                                                                                              gAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGAcaca                                             <  1:185499/250‑1 (MQ=255)                                                                                                                                                                                                   atatCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATAGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATGCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCtt                                         <  2:254507/249‑1 (MQ=255)                                                                                                                                                                                                    tatCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTCCTTCTGTTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACCGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACCCAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGCAATAAATTCGCTCCCCAGCTtgtg                                      >  1:99377/1‑251 (MQ=255)                                                                                                                                                                                                     atCAGGCGTATTGGAGAGCTGACTGTTGTGCCGAACGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGTCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACTTTGCTGTCTCTGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACGGAAATAAAGTAGGGACACAGCATGTgg                                     <  2:420635/251‑1 (MQ=255)                                                                                                                                                                                                     atCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTtgtg                                      <  1:600368/250‑1 (MQ=255)                                                                                                                                                                                                           cGTATTGAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGtt                               <  1:477424/251‑1 (MQ=255)                                                                                                                                                                                                              aTTGAAGAGCTGACTGGGGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTCAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACTCAGCTTGTGGCTGGTTTat                            >  1:181231/1‑251 (MQ=255)                                                                                                                                                                                                                 gAAGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTatca                          <  2:46512/250‑1 (MQ=255)                                                                                                                                                                                                                  aaGAGCTGACTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACCGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGttttttctc                        >  2:94391/1‑245 (MQ=255)                                                                                                                                                                                                                          aCTGGTGTGCCGATCGATATTATCTCTACCGGTCCGGATCGGACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTACCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTGATCATAAAGATaaa                <  2:99377/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCCGATCGATATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGGACTGAAATAAATTAGCGACCCAGCTTGTGGCTGGTTTATTATTAATATAAATGTAtttt        >  2:461881/1‑251 (MQ=255)                                                                                                                                                                                                                                     aatcgatATTATCTCTCCCGGTCCGGATCGTACTGAAACCATGATTCTGGGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCGGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCGTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTc     <  2:549069/250‑1 (MQ=255)                                                                                                                                                                                                                                        cgatATTATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCg  <  2:115016/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     CGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAAAGATCGTAGCGGCCTGCCGCAGGCGGCGCTGAACTATATCAAGCGTATTGAAGAGCTGACTGGTGTGCCGATCGATATCATCTCTACCGGTCCGGATCGTACTGAAACCATGATTCTGCGCGACCCGTTCGACGCGTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCG  >  NC_000913/4405688‑4406169

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