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%	231.8 / NA	62	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 (30/32);  total (30/32)

TTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCG  >  NC_000913/2732747‑2733216                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  ttAATGGTCTGCGTACCAGTAAGCGATGCTTTCTCTGCTGTCCACTCTTCTTCGAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACGCTTTCAGTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGATTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATTCTGGATGCGGCTTCatcaatca                                                                                                                                                                                                                               <  2:109160/250‑1 (MQ=255)                 cAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATcc                                                                                                                                                                                                              <  1:461554/251‑1 (MQ=255)                     aaGCGATGCCTTCTCTTCTGTCCTCGCTTCTTCTATCTCGGAGTACAGACGTTCTTTGTCGCTCAGTGCTTCGTTGAGCATATCCAGACGTTTTTTACTTGCTTCATCAGACTCTGTCATTAACGCCTGGGGTTCCAGGTTTAGCTGTAGGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAAGGCTGGATGCTGCTTCATCAATCAGTTCGATGGCTTTATCCGGCa                                                                                                                                                                                                          <  2:30105/251‑1 (MQ=255)                          aTGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAg                                                                                                                                                                                                         <  2:459874/247‑1 (MQ=255)                                 ctctGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTc                                                                                                                                                                                                >  2:161616/1‑249 (MQ=255)                                        ttCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCAAGTTTGAGCTGGATTATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGATGGATGCTGCTTCATCAAGCAGCCTGTGGGCTTTATCCGGCCGCTGACGATCAGCAATGt                                                                                                                                                                                        >  1:16183/1‑250 (MQ=255)                                                      tAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACg                                                                                                                                                                         >  2:161101/1‑251 (MQ=255)                                                       aaCTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTGTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTGCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGGCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACTGTCAGCAAGGTAGCGATGATACAAc                                                                                                                                                                          <  2:32000/249‑1 (MQ=255)                                                                       gTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTcgcc                                                                                                                                                                    >  1:94411/1‑239 (MQ=255)                                                                       gTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGGCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTcgcc                                                                                                                                                                    <  2:94411/239‑1 (MQ=255)                                                                         tCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTg                                                                                                                                                       <  1:22901/250‑1 (MQ=255)                                                                           tttttCGCTTAGTTTTTTGTTGATCATATTCAGTCGTCTTGTACTGGCTTTATCAGCCTCTTTTATTAACCCCTTTTGTTCCATTTTGAGCTGGATGATACGACGATCGATTCGTACGAGTTCTTCGGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCATTTCGATGGCTTTATCCTGCAGCTGACGGCCATCAATGTACCTATGTGCCAACGTCGACGCTGCAACAAATG‑ct                                                                                                                                                    <  2:505485/246‑2 (MQ=255)                                                                            ttGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCg                                                                                                                                                   <  1:129660/250‑1 (MQ=255)                                                                                 gCTCAGTTCTTCGTTGAGCATCTCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAg                                                                                                                                             <  2:409381/251‑1 (MQ=255)                                                                                   tCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCCGAATCTTTCATTAACGCCTGTTGTTCCAGTTTGAGATGGATGATCCGACGATCGCGCCGGTCGAGTTCTTCTGGTTTTCAGTCAATCTGCATACGCATGCCGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTAGCCGCTGCAACCATTG‑CCGGGTCAGt                                                                                                                                            >  1:245261/1‑250 (MQ=255)                                                                                      gTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAAtt                                                                                                                                        >  1:349708/1‑251 (MQ=255)                                                                                            cggttaGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTGTAGCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGcaca                                                                                                                                  <  2:419834/246‑1 (MQ=255)                                                                                                 aGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACAtggtg                                                                                                                             <  1:306345/251‑1 (MQ=255)                                                                                                  gCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACAtggtgg                                                                                                                            <  1:491283/251‑1 (MQ=255)                                                                                                  gCATATACAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGGTCCAGTTTGAGCTGGATGATACGACGATCGTGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCACCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGAGATGAGACAAAGACGCCGCTGCAACAATTGACCGGG‑CAGTACTTTGCACCtggagg                                                                                                                            >  2:236786/1‑248 (MQ=255)                                                                                                    atatCCAGCCGTTTTTTACTGTCTTCATCAGACTCGTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACTAGCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAAGCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGGCAGCAAGTTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGTTCAGTAATTTGCACATGgtggtg                                                                                                                          <  2:9548/251‑1 (MQ=255)                                                                                                        ccAGACGTTTTTTACTGGCTTCATCAGCCTCCTTCATTATCGCCTGTTGTTCACGTTTGCGCTGGATGATACGACGATCGAGTCGGTCGCGTTCTTCTGGTTTTGCGTCAATCTGCATACGAATGCTGGGTGCTGCTTCATCAATCAGGTCGTTGGCTTTATCCGGCCGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACACTTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAAt                                                                                                                      >  2:575851/1‑251 (MQ=255)                                                                                                         cAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAAtt                                                                                                                     <  1:3395/251‑1 (MQ=255)                                                                                                                tttttACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACg                                                                                                              >  1:147416/1‑251 (MQ=255)                                                                                                                   ttACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGtt                                                                                                            <  1:161101/250‑1 (MQ=255)                                                                                                                     aCTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATCCGACGCTAGAGTCGGTCGAGTTATTCTGGTTTTGAGTCAAACTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTAGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGAGATGAGACAACGTCGCAGCTGAAACAATTG‑CCGGGTCAGTAATTTGCACAGGGTGGTGCACTTCGTAACGTTCtt                                                                                                         >  2:357579/1‑251 (MQ=255)                                                                                                                     aCTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCtt                                                                                                         <  2:332125/251‑1 (MQ=255)                                                                                                                       tGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAAtt                                                                                                                                        <  1:266493/218‑1 (MQ=255)                                                                                                                            tcatcaGCCTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATCCGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGCTGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCACTTCGTAACGTTCTTTCAGGcc                                                                                                  >  1:431105/1‑251 (MQ=255)                                                                                                                            tcatcaGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGcc                                                                                                  >  1:120678/1‑251 (MQ=255)                                                                                                                                caGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACg                                                                                               <  2:230356/250‑1 (MQ=255)                                                                                                                                 aGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGGCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCa                                                                                             <  1:244883/251‑1 (MQ=255)                                                                                                                                 aGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATAAGAAGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGc                                                                                              >  1:563289/1‑250 (MQ=255)                                                                                                                                  gACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGAGCGAGTCGGTCGAGTTCTGCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGGCAGAAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAg                                                                                            <  2:350461/251‑1 (MQ=255)                                                                                                                                         tCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGc                                                                                      >  1:215519/1‑250 (MQ=255)                                                                                                                                              aaCGCCTGTTGTTACAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGt                                                                                <  1:361945/251‑1 (MQ=255)                                                                                                                                                cGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTCTCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCACCAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTACCGTTCTTTCAGGCCACGCAGAATCGCAATGGTAt                                                                              >  2:232011/1‑251 (MQ=255)                                                                                                                                                cGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGtct                                                                              >  1:458390/1‑249 (MQ=255)                                                                                                                                                        gttCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAAc                                                                       >  1:367470/1‑250 (MQ=255)                                                                                                                                                                    gCTGGAGGCTACGACGATCGCATCTGTCGAGTGCTTCTGGGTTTTAGTCAATCTGTAGACGTATGCTGGATGTTGCTTTATCAATCAGTTCTATGTCTTTATCCGGTAGCGGACGGTCAGCAATGTATCGATGAGACAACGTCGCCGCTGCAACAATTG‑CTGGGGCAGTAATTTGCGCATGTTGGTTCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCTATGGTATCTTCAACACAAGGCTCGGca                                                          <  2:506617/251‑1 (MQ=255)                                                                                                                                                                     cTGGATGATCCGACGATCGAGTCGGTCGAGTTCTTCTGGATTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGCATCGCAATGGTATCTTCAAAAGAAGGCTCGGca                                                          >  2:29944/1‑250 (MQ=255)                                                                                                                                                                            atACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGGCGGTAAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTA‑CCGGGTCAGTAATTTGCACCTGGTGGTGCAATTCGTAAAGTTCTTTAAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGATCGGCAACAAacac                                                  >  1:566608/1‑251 (MQ=255)                                                                                                                                                                              acgacgATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCGGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACt                                                 <  1:240147/250‑1 (MQ=255)                                                                                                                                                                               cgacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACttt                                               <  1:444400/251‑1 (MQ=255)                                                                                                                                                                                gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCCACAATTG‑CCGGGTCAGTAATTTGCACCTGGTGGTGCACTTCGTAACGTTCTTTCAGGCCa                                                                                                 >  2:133977/1‑200 (MQ=255)                                                                                                                                                                                gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTGTCAGGCCa                                                                                                 <  1:133977/200‑1 (MQ=255)                                                                                                                                                                                gacgaTCGAGTCGGTCGAGTTCTGCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTc                                              <  1:157065/251‑1 (MQ=255)                                                                                                                                                                                    aTCGAGTCGATCGAGATCTTAAGGTTTGGAGACAATCTGCCTACGAATGATGGATGATGCATCATCAATCAGGTCGATGGCTTTAGCAGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGCCGCCGCTGCACCAATTG‑CCGGGTCAGTAATATGCACATGGTGGTGAAATTCGTAAAGTTCTTTCAGGCCACGCCGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACGACCACTTGCTgg                                           >  2:297088/1‑250 (MQ=255)                                                                                                                                                                                      cGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaa                                         <  2:224837/250‑1 (MQ=255)                                                                                                                                                                                       gAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGacac                                       >  2:112183/1‑248 (MQ=255)                                                                                                                                                                                       gAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                       >  1:274095/1‑251 (MQ=255)                                                                                                                                                                                       gAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                       <  1:117415/251‑1 (MQ=255)                                                                                                                                                                                       gAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                       <  1:114120/251‑1 (MQ=255)                                                                                                                                                                                                 aGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTcc                               >  2:304634/1‑249 (MQ=255)                                                                                                                                                                                                     cttctGGTTTTGGGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTACAACAGAAGGCTCGGCAACAAACACATTCTGGACACGACGTTCCAgcgc                          >  1:61682/1‑250 (MQ=255)                                                                                                                                                                                                           gTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                   >  1:411111/1‑251 (MQ=255)                                                                                                                                                                                                              ttGAGTCAATCTGCATACGAATGCTGGATGCTGCTTAATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGCCGTTCCAGCGCAGCATCtttt                >  2:425863/1‑251 (MQ=255)                                                                                                                                                                                                                gAGTCAATCGGCATACGAATGCTGGATGCTGCGTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGATACAACTTCGCTGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGGGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttt               <  1:315449/250‑1 (MQ=255)                                                                                                                                                                                                                    cAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg          <  1:185024/251‑1 (MQ=255)                                                                                                                                                                                                                      aTCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGt         >  1:204814/1‑250 (MQ=255)                                                                                                                                                                                                                            atactaagtcTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTGCTGGCg  <  2:248836/242‑1 (MQ=255)                                                                                                                                                                                                                            aTACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCg  >  1:64676/1‑251 (MQ=255)                                                                                                                                                                                                                             tACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCg  >  1:405872/1‑250 (MQ=255)                                                                                                                                                                                                                                                 tcatcaATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTACCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCgg                                                                                                                                                  >  1:56120/1‑86 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  TTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTG‑CCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCG  >  NC_000913/2732747‑2733216

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