Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	253,835	G→A	intergenic (+674/+424)	yafP → / ← pepD	putative N‑acetyltransferase YafP/peptidase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	253,835	0	G	A	100.0%	170.0 / NA	46	intergenic (+674/+424)	yafP/pepD	putative N‑acetyltransferase YafP/peptidase D
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/24);  total (22/24)

AACATGCAGCGCGTAGTGGGGATGCCAGATTTACACCCCGGACGCGGCTACCCGATTGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCGGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGA  >  NC_000913/253590‑254077                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aaCATGCAGCGCGTAGTGGGGATGCCAGATTTACACCCCGGACGCGGCTACCCGATTGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAGACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCtc                                                                                                                                                                                                                                               >  2:7762/1‑251 (MQ=255)            cgTAGTGGGGATGCCAGATTTACACCCCGGACGCGGCTACCCGATTGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGCCACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGgcgc                                                                                                                                                                                                                                    >  2:246703/1‑251 (MQ=255)                               ttACACCCCGGACGCGGCTACCCGATTGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGc                                                                                                                                                                                                                  >  2:152720/1‑250 (MQ=255)                                    ccccGGACGCGGCTACCCGATTGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAGACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTggg                                                                                                                                                                                                            <  1:7762/251‑1 (MQ=255)                                                         tGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCa                                                                                                                                                                                        <  1:399440/250‑1 (MQ=255)                                                             gcAGCGTTCTTCTCTGTTGGTCGTTTTTGCCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAAc                                                                                                                                                                                   <  2:211976/251‑1 (MQ=255)                                                               aGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACa                                                                                                                                                                                  >  2:192615/1‑250 (MQ=255)                                                                  gTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCg                                                                                                                                                                              <  2:12506/251‑1 (MQ=255)                                                                   ttcttcTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCgg                                                                                                                                                                             >  1:300375/1‑251 (MQ=255)                                                                          ctGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCCGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAGGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGc                                                                                                                                                                         <  2:279538/248‑1 (MQ=255)                                                                                 tCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTAt                                                                                                                                                               >  2:88308/1‑251 (MQ=255)                                                                                    tttttACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATg                                                                                                                                                              >  2:120433/1‑249 (MQ=255)                                                                                         aCCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTGTCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCt                                                                                                                                                       <  2:33009/251‑1 (MQ=255)                                                                                           ccGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGc                                                                                                                                                     >  1:316344/1‑251 (MQ=255)                                                                                              gCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGtt                                                                                                                                                  <  2:220259/251‑1 (MQ=255)                                                                                                aaGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTtcg                                                                                                                                                <  1:206580/251‑1 (MQ=255)                                                                                                         cgGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCg                                                                                                                                       <  2:300375/251‑1 (MQ=255)                                                                                                                     tGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACAt                                                                                                                           >  2:373917/1‑251 (MQ=255)                                                                                                                      ggAAACAGAAACGGGCCTCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATg                                                                                                                          >  2:169607/1‑251 (MQ=255)                                                                                                                              aaaCGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAAGCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATaaa                                                                                                                  <  1:20022/251‑1 (MQ=255)                                                                                                                                 cGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAAcc                                                                                                               <  2:304628/251‑1 (MQ=255)                                                                                                                                      cGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGGCAACATGTCAATACAACCGACTc                                                                                                          >  2:103681/1‑251 (MQ=255)                                                                                                                                         tACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTcggc                                                                                                       >  2:389411/1‑251 (MQ=255)                                                                                                                                                                   ggTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCg                                                                             >  1:250282/1‑251 (MQ=255)                                                                                                                                                                      ttCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGt                                                                           <  1:192615/250‑1 (MQ=255)                                                                                                                                                                         tctcGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTAtt                                                                        >  1:67065/1‑250 (MQ=255)                                                                                                                                                                            cGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGTATGCAATCCGTTATTAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAATATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAg                                                                      <  1:406150/249‑1 (MQ=255)                                                                                                                                                                            cGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGc                                                                     <  2:186420/250‑1 (MQ=255)                                                                                                                                                                              aTGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTg                                                                  >  2:33900/1‑251 (MQ=255)                                                                                                                                                                              aTGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCCCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTg                                                                  >  2:198693/1‑251 (MQ=255)                                                                                                                                                                                        cGCATGGGCATTAAGCGAAAGCAGGAGGGGCACTAATCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTtcagt                                                        >  1:212858/1‑251 (MQ=255)                                                                                                                                                                                          cATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCAc                                                                                           >  1:332612/1‑214 (MQ=255)                                                                                                                                                                                          cATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTtcagtc                                                       >  2:75229/1‑250 (MQ=255)                                                                                                                                                                                          cATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCGGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCAc                                                                                           <  2:332612/214‑1 (MQ=255)                                                                                                                                                                                           aTGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTtcagtc                                                       >  2:72087/1‑249 (MQ=255)                                                                                                                                                                                                      gCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGGGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAg                                          <  1:189570/251‑1 (MQ=255)                                                                                                                                                                                                            gCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCt                                    <  2:37066/251‑1 (MQ=255)                                                                                                                                                                                                               ggTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCAGCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTGGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCAGGCTAAc                                 <  2:183052/251‑1 (MQ=255)                                                                                                                                                                                                                  gtgGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTTAGTCAGAGCGTAGTCAGCATGCTAACaa                               <  2:294041/250‑1 (MQ=255)                                                                                                                                                                                                                  gtgGCACTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAg                              <  1:218477/251‑1 (MQ=255)                                                                                                                                                                                                                       aCTATTCAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGGATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCt                         <  1:33900/251‑1 (MQ=255)                                                                                                                                                                                                                          aTTCATTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGGCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTg                        <  2:55980/249‑1 (MQ=255)                                                                                                                                                                                                                             cAGTGGATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACg                   <  1:120433/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGAtt                <  1:152720/250‑1 (MQ=255)                                                                                                                                                                                                                                  gATTTGTCCGAGTCCGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGc              >  1:123865/1‑251 (MQ=255)                                                                                                                                                                                                                                                cGTATCAGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCCCGCTTGATGATTGCCTGGa  >  2:325858/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   AACATGCAGCGCGTAGTGGGGATGCCAGATTTACACCCCGGACGCGGCTACCCGATTGGCGCAGCGTTCTTCTCTGTTGGTCGTTTTTACCCGGCAAGACGTCGCGGTAACGGTGCTGGAAACAGAAACGGGCCGCTACTCTGACACACTGCGTTCGGCGCTGGTTTCTCTCGATGGCGATAACGCATGGGCATTAAGCGAAAGCTGGTGTGGCACTATTCAGTGGATTTGTCCGAGTCCGTATCGGCCTCATCATGGGCGCAAAAACTGGTTTCTGGGCATTGGGCGTTTTACCGCTGATGAGCAGGAACAATCGGATGCAATCCGTTATGAGACGCTGCGTTCGTCGGGGCCGGGCGGTCAACATGTCAATAAAACCGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGA  >  NC_000913/253590‑254077

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