Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,817,197	G→A	R103C (CGC→TGC)	yqaB ←	fructose‑1‑phosphate phosphatase YqaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,817,197	0	G	A	100.0%	151.7 / NA	42	R103C (CGC→TGC)	yqaB	fructose‑1‑phosphate phosphatase YqaB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/22);  total (20/22)

GGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGTA  >  NC_000913/2816958‑2817441                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     ggCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGcc                                                                                                                                                                                                                                           <  2:26157/251‑1 (MQ=255)   cGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCaa                                                                                                                                                                                                                                         <  1:171486/251‑1 (MQ=255)    gTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAAc                                                                                                                                                                                                                                        >  1:339903/1‑251 (MQ=255)                         cgccTGAATACCGCACTCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATcaacacgcg                                                                                                                                                                                                                   >  1:91464/1‑247 (MQ=255)                         cgccTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAagag                                                                                                                                                                                                                   <  2:258590/251‑1 (MQ=255)                                aTACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                                                                                                                                                            >  1:350615/1‑251 (MQ=255)                                      aaaTCTGCATCTTCAAAGACCACTCACTGCGTCGGTTTCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGTCGCGGGTTTATGGAGTTTGACGTGATCGGCAGCTACGACGGCTGCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCATCGATGGCGCTTTCTCTCCCCGTTCCTACAGCCATTGGGCAACTACCATGCCATCTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                                                                                                                                                      <  1:191330/251‑1 (MQ=255)                                            gCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGCCGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCACCGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATc                                                                                                                                                                                                >  2:301261/1‑251 (MQ=255)                                                   cAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACAGCAACAAGAGGAAGCGGTTCGACGCTATCCagcagc                                                                                                                                                                                         <  2:91464/251‑1 (MQ=255)                                                     aaGACCACACACTGCGGCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGGTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCAt                                                                                                                                                                                       <  1:399203/251‑1 (MQ=255)                                                       gACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATAc                                                                                                                                                                                     <  1:9753/251‑1 (MQ=255)                                                         ccACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActt                                                                                                                                                                                   <  1:140605/251‑1 (MQ=255)                                                                   cggctgTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCGTGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATTCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCAGACTTCCTACTGCtt                                                                                                                                                                         <  2:395158/246‑1 (MQ=255)                                                                                                                    gcgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAg                                                                                                                        <  1:17803/251‑1 (MQ=255)                                                                                                                      gggTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCt                                                                                                                      >  2:382095/1‑251 (MQ=255)                                                                                                                        gTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTc                                                                                                                     >  1:387481/1‑250 (MQ=255)                                                                                                                            aTGGTGTTTGACGTGATCGGCAGCGACGTCGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCACTGCCTCAGCGATGGCGCTTTCACTCCCAGTTCCTACAGCCATTGGGCAACGCCCATGCCAACTTTTCACCACATCAACAAGCGGAAGCGGTTCGACGATCTCCAGCAGCCTACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACACATGCGGGTCGAGATCGGCCTGATTCAGCTCaata                                                                                                                 >  2:197634/1‑250 (MQ=255)                                                                                                                              ggTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAAtt                                                                                                              <  1:411393/251‑1 (MQ=255)                                                                                                                                     gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGAATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTg                                                                                                         >  1:250677/1‑249 (MQ=255)                                                                                                                                           aTCGGCAGCGACGCCGGCGTCAAAATAATGGCGTAATCCCTGTGGCGCAAGCAATGCCTCAGCGATGGCGCTTTCACGCCCCGTTCCTACAGCCATTGGGCAACGTCCATGCCAACTTTTCCCCACATCAACAAGAGGAAGCGGTGCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTGTTGACGCGCTTACGCATGCGGGTCGAGATCGGCCTTATTCAGCTCAACAATTGccggcgtaatac                                                                                                 <  2:184073/251‑11 (MQ=255)                                                                                                                                           aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATAc                                                                                                 <  1:214603/251‑1 (MQ=255)                                                                                                                                             cGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGc                                                                                               >  2:189118/1‑251 (MQ=255)                                                                                                                                               gCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCa                                                                                             <  2:279686/251‑1 (MQ=255)                                                                                                                                                cAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCATCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACAGCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAg                                                                                            <  1:50642/251‑1 (MQ=255)                                                                                                                                                cAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGcc                                                                                              >  1:301619/1‑249 (MQ=255)                                                                                                                                                  gCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTAGCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGt                                                                                          <  1:412703/251‑1 (MQ=255)                                                                                                                                                    gacgacGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTgg                                                                                        >  2:78681/1‑251 (MQ=255)                                                                                                                                                             gTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGAGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCAt                                                                               >  1:368959/1‑251 (MQ=255)                                                                                                                                                                 aaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTAAGCTCAATAATTGCCTGAGCAGTACGCCAGGAGGGCGATCAATTa                                                                             >  1:220237/1‑249 (MQ=255)                                                                                                                                                                        ggCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAatc                                                                    <  2:159260/251‑1 (MQ=255)                                                                                                                                                                                cccAGGTGCGCCAGCAAGGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACGGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCATCTCAATAATTGCCTGATCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCt                                                             <  2:364825/250‑1 (MQ=255)                                                                                                                                                                                       gcgcCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATc                                                      <  1:74623/250‑1 (MQ=255)                                                                                                                                                                                              cAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGa                                               >  1:21005/1‑250 (MQ=255)                                                                                                                                                                                                 tGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGAc                                           >  1:187139/1‑251 (MQ=255)                                                                                                                                                                                                             ggCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttttt                                                                                                                                                                >  1:37678/1‑122 (MQ=255)                                                                                                                                                                                                             ggCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttttt                                                                                                                                                                <  2:37678/122‑1 (MQ=255)                                                                                                                                                                                                                                 cTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCttt            >  2:105348/1‑250 (MQ=255)                                                                                                                                                                                                                                   aCAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACg         <  1:154333/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTg      >  1:39141/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGc     <  2:196380/251‑1 (MQ=255)                                                                                                                                                                                                                                          ttGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCATTAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTCGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGtg  >  2:98426/1‑250 (MQ=255)                                                                                                                                                                                                                                           tGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGt   >  2:410637/1‑249 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     GGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGTA  >  NC_000913/2816958‑2817441

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