Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,592	C→A	P294P (CCG→CCT)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,592	0	C	A	100.0%	233.2 / NA	68	P294P (CCG→CCT)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/34);  total (34/34)

GGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCC  >  NC_000913/2446346‑2446837                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      ggCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTA‑CAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGtt                                                                                                                                                                                                                                                   <  1:334778/250‑1 (MQ=255)           gCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCat                                                                                                                                                                                                                                         >  1:201715/1‑251 (MQ=255)                cagcGAAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATggg                                                                                                                                                                                                                                    >  1:116964/1‑251 (MQ=255)                   cgcAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCaa                                                                                                                                                                                                                                 <  1:89251/251‑1 (MQ=255)                         cgcgGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGAtt                                                                                                                                                                                                                            <  2:189927/250‑1 (MQ=255)                          gcgGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTtg                                                                                                                                                                                                                          <  2:35425/251‑1 (MQ=255)                           cgGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTtgc                                                                                                                                                                                                                         <  1:523516/251‑1 (MQ=255)                             gTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCTACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTtgctg                                                                                                                                                                                                                       >  1:215202/1‑251 (MQ=255)                                      aTTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCgctgct                                                                                                                                                                                                              >  1:39157/1‑251 (MQ=255)                                        tttttACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGctgctg                                                                                                                                                                                                             >  1:295834/1‑250 (MQ=255)                                             cccAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccg                                                                                                                                                                                                       <  2:116964/250‑1 (MQ=255)                                               cAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccgcc                                                                                                                                                                                                     <  1:171188/251‑1 (MQ=255)                                                aGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGTAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATCccgccg                                                                                                                                                                                                    >  2:74307/1‑251 (MQ=255)                                                  cGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCgag                                                                                                                                                                                                  <  1:220121/251‑1 (MQ=255)                                                  cGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGCTGATCCATTAAAACCATCGCCAGCATCGCTTCTGCGATCGGCACTGCCAGGATCCAGACACAGGGATCGTGACGGCATTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACACGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccgccg                                                                                                                                                                                                    >  1:521993/1‑249 (MQ=255)                                                    gggAATATCAGTCTTCACATCGGCATTTGGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGaa                                                                                                                                                                                                <  1:553019/250‑1 (MQ=255)                                                     ggaatttcatccttcactgcGGCATTTTGCTCCCGTTGCCGTACCAGGGGATCCATTAAAACGAGCGCCGGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCGTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                                                                                                                                               <  2:423314/232‑1 (MQ=255)                                                             aGTCTTCCCATCGGCATTTTGCGCCCGTTGCCGAAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCGGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCTCGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGAGTTGCTGCCCGCTGCTGATACCGCCTAGAATGccgcccg                                                                                                                                                                                       <  2:372973/251‑1 (MQ=255)                                                                      aTCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGc                                                                                                                                                                              >  1:13116/1‑251 (MQ=255)                                                                          gCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTg                                                                                                                                                                          >  2:150945/1‑251 (MQ=255)                                                                                 gcgcCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCg                                                                                                                                                                   >  1:117558/1‑251 (MQ=255)                                                                                  cgcCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGt                                                                                                                                                                  >  1:370056/1‑251 (MQ=255)                                                                                   ccccGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCGGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTc                                                                                                                                                                 <  1:39263/250‑1 (MQ=255)                                                                                      cGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCttt                                                                                                                                                              >  2:178238/1‑251 (MQ=255)                                                                                      cGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCttt                                                                                                                                                              >  2:124288/1‑251 (MQ=255)                                                                                          gCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAAGGATTTGCTGCCCGCTGCTGATTCCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGt                                                                                                                                                           <  1:202554/250‑1 (MQ=255)                                                                                            cGTAACAGGTGATCCATTAAAACGATCGCCGGCATCGCTTCTGCGATCGGCACGGCGCGGACCCCGACAAAGGGATCGTGCCGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGATACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAt                                                                                                                                                        <  2:407452/251‑1 (MQ=255)                                                                                            cGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCACCTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAt                                                                                                                                                        >  1:372465/1‑251 (MQ=255)                                                                                              tAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAttt                                                                                                                                                      >  1:251431/1‑251 (MQ=255)                                                                                              tAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAttt                                                                                                                                                      >  1:484916/1‑251 (MQ=255)                                                                                               aaCAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACTGTCTTTGGTGAttt                                                                                                                                                      <  1:116863/250‑1 (MQ=255)                                                                                                    gTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATcgc                                                                                                                                                >  1:310965/1‑251 (MQ=255)                                                                                                           aTTAACACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGGGGCTCCCGACACAGGGCTCGTTACGGCCTTTGGTGATCATCTCAACTTCTCCGCCAAAGGGGTTAATGGTCCGACCCGGCACGGTAATGCTGGAGGTAGGTGTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTg                                                                                                                                         <  2:501080/251‑1 (MQ=255)                                                                                                                                 cTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACg                                                                                                                   >  1:327370/1‑251 (MQ=255)                                                                                                                                        gATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGc                                                                                                            >  2:133853/1‑251 (MQ=255)                                                                                                                                           cGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGt                                                                                                         <  2:509065/251‑1 (MQ=255)                                                                                                                                               aCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGcc                                                                                                     >  1:445503/1‑251 (MQ=255)                                                                                                                                                cTGCGCGGATCCCGACACAGGGATCGTGACGGCATTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTCCGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGcc                                                                                                     >  2:296194/1‑250 (MQ=255)                                                                                                                                                  gcgcGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAGGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCaa                                                                                                   <  2:85746/250‑1 (MQ=255)                                                                                                                                                               cacaGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGGCGCCAATTTCCACGCCTTTc                                                                                     <  1:496105/251‑1 (MQ=255)                                                                                                                                                                 caGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGGAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCTCGCCTTTCAc                                                                                   <  2:357118/251‑1 (MQ=255)                                                                                                                                                                  aGGGATCGGGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGAGTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCAcc                                                                                  <  1:94264/251‑1 (MQ=255)                                                                                                                                                                    ggATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGTTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACcgc                                                                                <  1:140691/251‑1 (MQ=255)                                                                                                                                                                      aTCGGGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGTTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGt                                                                              <  2:61485/251‑1 (MQ=255)                                                                                                                                                                        cGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTg                                                                            <  1:96454/251‑1 (MQ=255)                                                                                                                                                                            aCGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGAGGTTAATGGTACGACCCGGCACGGACATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGc                                                                        >  2:111229/1‑251 (MQ=255)                                                                                                                                                                                   ttGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCtcatc                                                                   <  2:249583/249‑1 (MQ=255)                                                                                                                                                                                   ttGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAAGGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCtcatca                                                                  <  2:104962/250‑1 (MQ=255)                                                                                                                                                                                     ggTGATCAGCTCACCTTCTTCGCCAAAGCGGTTAAGGGTACGACCCGGCACGGTAATGCGGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAgc                                                                <  2:9271/250‑1 (MQ=255)                                                                                                                                                                                          tcatcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGCGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATgg                                                          >  1:454432/1‑251 (MQ=255)                                                                                                                                                                                          tcatcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATgg                                                          >  1:297990/1‑251 (MQ=255)                                                                                                                                                                                                 aaCTTCTTCGCCAAAGCGGTTAATGGGACGACCCGGCGCGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGt                                                   <  2:573614/251‑1 (MQ=255)                                                                                                                                                                                                        tcGCCAAAGCGGTTAATGGTCCGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCAt                                             <  2:442188/250‑1 (MQ=255)                                                                                                                                                                                                          gCCAAAGCGGTTAATGTTACGCCCCGGCCCGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGTTTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCa                                          <  1:540688/251‑1 (MQ=255)                                                                                                                                                                                                             aaaGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGc                                       >  1:262743/1‑251 (MQ=255)                                                                                                                                                                                                                 cGGTTCATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTc                                   >  1:485931/1‑251 (MQ=255)                                                                                                                                                                                                                       aTGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCAGCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGAc                             <  1:576241/251‑1 (MQ=255)                                                                                                                                                                                                                        tGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGAcc                            >  2:161719/1‑251 (MQ=255)                                                                                                                                                                                                                         ggTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCAGCCAGGCGGTCAAAGAcc                            <  1:178238/250‑1 (MQ=255)                                                                                                                                                                                                                           tACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGc                         >  1:411075/1‑251 (MQ=255)                                                                                                                                                                                                                              gACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCTCCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTcc                      >  2:222416/1‑250 (MQ=255)                                                                                                                                                                                                                               aCCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGc                     >  1:424952/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGt                >  1:494145/1‑251 (MQ=255)                                                                                                                                                                                                                                     cACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGCTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTccgcccagc               >  2:397153/1‑243 (MQ=255)                                                                                                                                                                                                                                         gTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGc           >  1:466363/1‑251 (MQ=255)                                                                                                                                                                                                                                             tGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAgg        >  1:32513/1‑250 (MQ=255)                                                                                                                                                                                                                                                 ggAGGTAGGTTTCAGCGCCATAGGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAAc     <  1:466597/249‑1 (MQ=255)                                                                                                                                                                                                                                                  gAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGcc  <  1:130613/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCC  >  NC_000913/2446346‑2446837

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