Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,732,992	G→A	I394I (ATC→ATT)  I394I (ATC→ATT)	clpB ← clpB ←	chaperone protein ClpB chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,732,992	0	G	A	100.0%	236.1 / NA	70	I394I (ATC→ATT)  I394I (ATC→ATT)	clpB clpB	chaperone protein ClpB chaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/30);  total (40/30)

TTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCG  >  NC_000913/2732746‑2733233                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  tttAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATc                                                                                                                                                                                                                                                >  2:238918/1‑250 (MQ=255)  ttAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCatcaatca                                                                                                                                                                                                                                               <  2:449762/250‑1 (MQ=255)  ttAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATc                                                                                                                                                                                                                                                >  1:407622/1‑249 (MQ=255)       ggTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCg                                                                                                                                                                                                                                          >  2:34094/1‑250 (MQ=255)            gCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATgg                                                                                                                                                                                                                                      <  2:451159/249‑1 (MQ=255)                 ccAGAAATCGATGCCTTCTCTGCTTTCCACTCTTCTTCTATCTCGGATTACGGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGCCGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTTGATGCTGCTTCATCAATAAGGTCGATGGCTTTATc                                                                                                                                                                                                                               <  2:23641/251‑1 (MQ=255)                      aaGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTCACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGc                                                                                                                                                                                                                           >  2:364598/1‑250 (MQ=255)                               cTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGt                                                                                                                                                                                                                 >  2:259147/1‑251 (MQ=255)                                ttctctgatcttcACTCTTATTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCCGTTCTTCGTTGAGCATATCCAGACGTTTTTTTCTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCACTCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGCTGGCTTTATCCGGCAGCTGACGGTc                                                                                                                                                                                                                >  2:383204/13‑251 (MQ=255)                                                        aaCTCGGAGTCCTGGCGTTCTTTTCCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTTATGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTGCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCGGCAGCAATCAGGTCGATGGCGTTGTCCGGCAGCTTACGGTCAGCAATGTAGCTATGAGACAACg                                                                                                                                                                                         <  2:492831/250‑1 (MQ=255)                                                            cGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTcgcc                                                                                                                                                                                    <  1:34094/251‑1 (MQ=255)                                                               aGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCt                                                                                                                                                                                 <  2:390487/251‑1 (MQ=255)                                                                gTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTg                                                                                                                                                                                <  2:4905/251‑1 (MQ=255)                                                                gTACTGACGTTCTTTGTCGCTCAGTTCTGCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATTAGACAATGTCGGCGCTg                                                                                                                                                                                <  2:268308/251‑1 (MQ=255)                                                                      aCGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCATCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGcaacaa                                                                                                                                                                          <  1:370305/251‑1 (MQ=255)                                                                           cTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGcc                                                                                                                                                                     >  1:103468/1‑251 (MQ=255)                                                                              tGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACACCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGACGCCGCTGCAACAATTGCCggg                                                                                                                                                                  >  2:564331/1‑251 (MQ=255)                                                                                      aGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGAt                                                                                                                                                                                                  <  2:87801/211‑1 (MQ=255)                                                                                      aGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGAt                                                                                                                                                                                                  >  1:87801/1‑211 (MQ=255)                                                                                          cttcGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCCGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGGTGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGAAAACGTCGCCGCGGCACCAATTGCCGGGACAGTAATTTGc                                                                                                                                                      >  2:35255/1‑251 (MQ=255)                                                                                             cGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGcac                                                                                                                                                    <  1:169954/250‑1 (MQ=255)                                                                                                tGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACAtgg                                                                                                                                                >  2:470913/1‑251 (MQ=255)                                                                                                    cATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATggtggt                                                                                                                                            >  1:8660/1‑251 (MQ=255)                                                                                                     atatCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGgtggtg                                                                                                                                           >  1:266299/1‑251 (MQ=255)                                                                                                      catCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGc                                                                                                                                          >  1:70426/2‑251 (MQ=255)                                                                                                          cAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGACTGATGCTCCAGTTTGAGCTGGATAATACGACCATCGAGACGGTAGAACTCTACTGGTTTTGAGCCAATCTGCATACGAATGAAGGATGCTGCTTCATCAATAAGGTCGATGGCTTTATCGGGCAGCTGACGGTCCGCAATGTAGCGAAGAACCAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTACACATGGAGGTGCactt                                                                                                                                      >  1:547613/1‑248 (MQ=255)                                                                                                           aGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTc                                                                                                                                     <  1:411783/251‑1 (MQ=255)                                                                                                                      aCTGGCTTCCTCAGACTCTTTCATTAACGACTGTTGTTCCAGTTTGCGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCCGGTCGATGGCTTTATACGGCAGCTGACGGTCAGCAATGTAGCGATGAGAAAACGTCTCCGCTGCAACAATTGACGGGTAAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCtt                                                                                                                          >  2:339112/1‑251 (MQ=255)                                                                                                                      aCTGGCTTCATCAGACTCTTTCATTACCGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGt                                                                                                                              <  2:190405/247‑1 (MQ=255)                                                                                                                      aCTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGt                                                                                                                              >  1:190405/1‑247 (MQ=255)                                                                                                                         ggCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACACTTGCCGGGCCAGTAATTTGCACATGGTGGTGCCCTTCGTAACGTTCTTTc                                                                                                                        >  1:570113/1‑250 (MQ=255)                                                                                                                          gCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGTGCTGGATGATACGACGATCGAGTCGGTCGAGTTATTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCAAATGGTGGTGCCATTCGTAACGTTCTTTc                                                                                                                        >  2:474141/1‑249 (MQ=255)                                                                                                                          gCTTCATCAGACTCTTTCATTAACGCCTGTTGATCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTc                                                                                                                        >  1:342035/1‑249 (MQ=255)                                                                                                                                 caGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGc                                                                                                               >  1:195790/1‑251 (MQ=255)                                                                                                                                  cgACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGc                                                                                                               >  2:71427/2‑250 (MQ=255)                                                                                                                                         ttCATTAACGCCTGTTGTTCCAGTTTGAGTTGGATGATACGACGATCGCGTCGGTCGAGTTCTTCTGGATTCAAGCAAATCTGCGCACGAATGCTGGCTCATGATTAATCAATCAGAACGCTGGATTTATCCGGCCGCTCACGGTCAGCAATGTAGAGATGAGACAACGTCGCCGCTGCAAAAATTGACGGGTCACCAATTTGCAACTGGTTGTGCAATTCGTAACGATCTTTCAGACCTCGCAGAAACGc                                                                                                       >  1:575761/1‑251 (MQ=255)                                                                                                                                         ttCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGc                                                                                                       >  1:318540/1‑251 (MQ=255)                                                                                                                                         ttCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGc                                                                                                       >  2:475921/1‑251 (MQ=255)                                                                                                                                          tCATTAACGCATGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGCGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCCGCTGACGGTCAGCAATGTAGCGAAGAGACAACGTCGCAGCTGCAACAATTGCCGGATCAGAAATTTGCACATGGTGGTGCAATTCGTAACGTTATTTCAGGCCACGCAGAATCGc                                                                                                       >  2:556700/1‑250 (MQ=255)                                                                                                                                            aTTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTTCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGCCGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATTGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAAt                                                                                                    >  2:60469/1‑251 (MQ=255)                                                                                                                                               aaCGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCCTACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTGTCCGGCAGCTGAAGGTCAGAAATGTAGCGATGAGACAACGTCGCCTCTGCAACAATTGCCGGGTCAGTCATTTGCACATGGTGGTGCACTTGGTAACGTGCTTTCAGGCCCCGCAGAATCGCAATGGt                                                                                                 >  1:534555/1‑251 (MQ=255)                                                                                                                                               aaCGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGt                                                                                                 <  1:364598/251‑1 (MQ=255)                                                                                                                                                  gCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATc                                                                                              >  1:478026/1‑251 (MQ=255)                                                                                                                                                   ccTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTTGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGGCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCt                                                                                             <  1:35255/251‑1 (MQ=255)                                                                                                                                                           tCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACCTGGTGGAGCAATTCGTAACGTTCTTTCAGGCCACGCAGACTCGCAATGGTATCTTCAAc                                                                                        >  1:540835/1‑248 (MQ=255)                                                                                                                                                                   gAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGTTATCTTCAACAGAAGGCGCgg                                                                             <  2:43381/251‑1 (MQ=255)                                                                                                                                                                        ggATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaac                                                                         <  1:524196/250‑1 (MQ=255)                                                                                                                                                                         gatgatACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaacaa                                                                       <  2:478026/251‑1 (MQ=255)                                                                                                                                                                                 gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACACTTGCCGGGTCAGTAATTTGCACCTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCACTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTc                                                               >  2:194205/1‑251 (MQ=255)                                                                                                                                                                                 gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAg                                                                                   >  2:456998/1‑231 (MQ=255)                                                                                                                                                                                 gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGGATCTGCAACAGAAg                                                                                   <  1:456998/231‑1 (MQ=255)                                                                                                                                                                                 gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATAAGAACGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCCGAATCGCAATTGTATCTTCAACAGTAGACTCGGCAACAAACACTTTc                                                               >  1:64708/1‑251 (MQ=255)                                                                                                                                                                                 gacgaTCGAGTCGGTCGAGTTCTTCTGGTGTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTc                                                               <  1:238918/251‑1 (MQ=255)                                                                                                                                                                                   cgaTCGAGTCGTTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGAGGCTGCTTCATCAATCAGTTCGATGGCTTTATCCGGCAGCTGACGGGCAGCAATGTAGCGATGAGACAACGGCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTg                                                             <  1:560714/251‑1 (MQ=255)                                                                                                                                                                                       cGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaa                                                          >  1:216452/1‑250 (MQ=255)                                                                                                                                                                                        gAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                                        <  1:71427/251‑1 (MQ=255)                                                                                                                                                                                            cggtcgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCGCGGCAACAAACACTTTCTGGAAacgacg                                                    <  2:8660/251‑1 (MQ=255)                                                                                                                                                                                             ggtcgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGt                                                   <  1:470913/251‑1 (MQ=255)                                                                                                                                                                                                cgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTcc                                                >  2:324700/1‑251 (MQ=255)                                                                                                                                                                                                cgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGAAGGTCAGCACTGTAGCGATGAGACAACATCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTcc                                                >  1:388488/1‑251 (MQ=255)                                                                                                                                                                                                 gAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTca                                                >  1:530857/1‑249 (MQ=255)                                                                                                                                                                                                 gAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCa                                               <  1:485931/251‑1 (MQ=255)                                                                                                                                                                                                       ttctGGTTTTGAGTCAATCTGCATACGAATGCGGGATGCTGCTGCAGCAAACAGGTCGAGGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGAAAACGTCGCCGATGCAACAATTGACGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTCCAGGCCACGCAGAATCGCAAGGGTATCTTCAACAGAAGGCGCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgca                                          <  2:83183/250‑1 (MQ=255)                                                                                                                                                                                                        tctGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAgcgc                                           <  2:342035/248‑1 (MQ=255)                                                                                                                                                                                                           ggTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcagca                                       <  2:161355/249‑1 (MQ=255)                                                                                                                                                                                                            gTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                    <  1:458553/251‑1 (MQ=255)                                                                                                                                                                                                               ttGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTAGCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCGCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttt                                  <  1:372300/250‑1 (MQ=255)                                                                                                                                                                                                                     cAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                           >  1:104963/1‑251 (MQ=255)                                                                                                                                                                                                                               acggaTGCTGGATGCCGCTGCAGCAAGCAGGGCGATGGCTTGAGCCGGCAGCGGACGGTCAGCAATGTAGCGGTGAGACAATGTCGCCGCTGCAAGGATTGGCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCATAATCGCAATGGTATCTGCAACAGAAGGCTCGGCAACAAACACTGTCTGGGAACGACGTTCCAGCGCAGCTTCTTTTGCAATGTACTGGCGat                 <  2:373609/247‑1 (MQ=255)                                                                                                                                                                                                                                           tgctgcTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCCATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcg     >  2:4818/1‑251 (MQ=255)                                                                                                                                                                                                                                              tgctTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGCTATTCGTCAAGcgtcg  >  2:186608/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  TTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCG  >  NC_000913/2732746‑2733233

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