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%	140.7 / NA	42	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 (20/22);  total (20/22)

TAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGC  >  NC_000913/2732748‑2733237                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       tAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑gg                                                                                                                                                                                                                                                 >  1:69821/1‑251 (MQ=255)           cGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCtt                                                                                                                                                                                                                                       >  1:300073/1‑251 (MQ=255)                          ggacTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGTCGTTTTTTCCTGGCTTCTTCAGACTCTTTCATTAACGCCTTTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGTTTTTTAGTCAATCTTCATATTAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATTGCTTTATCCGGCTGCTGAc                                                                                                                                                                                                                        <  2:354984/248‑1 (MQ=255)                                  ctGCTTTCCACCCGTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTTATTAACTCCGGTTGTTCCAGTTTGAGCTTGATGATACGACGCTCGCGTCGGTTGAGTTCTTCTGGTTTTTAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGATGGTCa                                                                                                                                                                                                                   <  1:53216/248‑1 (MQ=255)                                      tttCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCCGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATAGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGAAGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGt                                                                                                                                                                                                            >  1:31601/1‑251 (MQ=255)                                        tCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGGGATCTTCTTGTTTTGCGTCCATCTGCCGACGGATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCACCTAACGGTCCGCAAt                                                                                                                                                                                                              >  2:227438/1‑247 (MQ=255)                                        tCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAAt                                                                                                                                                                                                              <  1:227438/247‑1 (MQ=255)                                         ccgctctTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGTTCAGTTCTGC‑GTTGATCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCAGTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGCCAATCTGCATGCGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAg                                                                                                                                                                                                          <  1:201490/247‑1 (MQ=255)                                          cACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCg                                                                                                                                                                                                        >  1:164853/1‑251 (MQ=255)                                                 cttctAACACGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGAc                                                                                                                                                                                                 >  1:256534/1‑251 (MQ=255)                                                      aaCTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑TGTCGATTGCTTTATCCGGCATCTGACGGTCAGCAATGTTTCGATGAGACAAc                                                                                                                                                                                              <  2:28293/249‑1 (MQ=255)                                                       attCGGAGTACTGACTTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTTCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTc                                                                                                                                                                                           <  2:115836/249‑1 (MQ=255)                                                          cGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCCGTT‑ATCATATCCAGACGTGTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTAGCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTcgcc                                                                                                                                                                                        <  2:129234/251‑1 (MQ=255)                                                                                    aGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAAt                                                                                                                                                              <  1:161442/251‑1 (MQ=255)                                                                                       tcttc‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTg                                                                                                                                                           >  1:222390/1‑251 (MQ=255)                                                                                           c‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGcaca                                                                                                                                                       <  1:8155/251‑1 (MQ=255)                                                                                               tGAGCATATCCAGTCGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTg                                                                                                                                                                           <  2:247368/228‑1 (MQ=255)                                                                                               tGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTg                                                                                                                                                                           >  1:247368/1‑228 (MQ=255)                                                                                                   cATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACAtggtggg                                                                                                                                                >  1:145826/1‑250 (MQ=255)                                                                                                                tttttACTGGCTTCATCAGACTCTTTCAGTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAAc                                                                                                                                    <  2:11189/250‑1 (MQ=255)                                                                                                                        ggCTGCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGGCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTc                                                                                                                            <  2:324145/250‑1 (MQ=255)                                                                                                                             catcaGACTCTTTCATTAACGCCTGTTGTGCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCATTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCa                                                                                                                      <  1:207546/251‑1 (MQ=255)                                                                                                                                  gACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGCATGCTGGATGCTGCTTCCTCAATCA‑GGACGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGCGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGCAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGc                                                                                                                   >  1:171763/1‑249 (MQ=255)                                                                                                                                        ttCATTAACGGCTGTTGTTCCAATTTAAGCCTGATGATACGTCGGCCGAGTCGGTCGAGTTCTTCTGGGTTTGAATCAATCTGCCTTCGAATGCTGGATTCTGCTTCATCAATCACGGT‑GATGGTTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGGGCAATTCGTACCGTTCTTTCAGGCCACGCAGAATca                                                                                                            <  2:220109/250‑2 (MQ=255)                                                                                                                                        ttCATTAACGCCTGTTGTTCCCGTTTGAGCTGGATGATACTACGATCGAGTCGGGCGAGTTCTTCTGGTTGTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCTGGTCAGTAATTTGCACATGGTGTTGCAATTCGTAACGTTCTTTCAGGCCCCGCAGAAt                                                                                                              <  2:297901/248‑1 (MQ=255)                                                                                                                                          cATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGCTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTTCCGGGTCAGTAATTTGCACATGGTGGTGCCATTCGTAACGTTCTTTCAGGCCACGCAGCATCGCa                                                                                                          >  2:92645/1‑250 (MQ=255)                                                                                                                                                cGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGTGTCGGTCTAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCAGCAATCA‑GGTCGATGGCTTTATCCGGCGGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAATAATTGCCGGGTCAGTAAGTTGCACATGGTGGTGCAATGCGTAACGTCCTGTCAGGCCACGCAGAATCGCACTGGTAt                                                                                                   <  2:171636/251‑1 (MQ=255)                                                                                                                                                  ccTGTTGTTCCAGTTTGAGCTGGATGATACGACGCTCGAGTCGGTCGAGTTCTTCTGGTATTGAGTCAATATGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAAACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCt                                                                                                 >  1:41798/1‑251 (MQ=255)                                                                                                                                                   cTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACCATTGCCGGGTCAGTACTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCtt                                                                                                >  2:391727/1‑251 (MQ=255)                                                                                                                                                     gttgttCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGCTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTCGCGATGAGACAACGTCGCCGCTGCAACCATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTc                                                                                               >  1:59293/1‑250 (MQ=255)                                                                                                                                                        gttCCAGTTTGATCTGGGTGATACTTCGCTCGAGTCGGTCGATTTCTTCTTGTTTTGAGTCAATCTGCATACGACTGCTGGATGCTGCCTCATCAATCA‑GGTCGATGGCTTTAGCCGGCAGCTGACTGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCTGGTCAGTTATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGTCCACGTAGAATAGCAATGGTATATTCatc                                                                                            <  2:337550/250‑3 (MQ=255)                                                                                                                                                                  gAGCTGGCTGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGCATGCTGGGTGCTGCTTCCTCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGGGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCgg                                                                                 >  1:387889/1‑251 (MQ=255)                                                                                                                                                                  gAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCgg                                                                                 >  1:349227/1‑251 (MQ=255)                                                                                                                                                                                gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCCTACGAATCCTGGATGCTGCTTCCTCAATCA‑GGTCGATGGCTTTATCCCGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCCACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTc                                                                   >  1:205784/1‑251 (MQ=255)                                                                                                                                                                                                                      aTCTGCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGt                              >  1:15044/1‑250 (MQ=255)                                                                                                                                                                                                                          gCATACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTgg                         <  1:92645/251‑1 (MQ=255)                                                                                                                                                                                                                            aTACGAATGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCg                       <  1:202910/251‑1 (MQ=255)                                                                                                                                                                                                                                   tGCTGGATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCg                >  1:50741/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc              <  1:2025/249‑1 (MQ=255)                                                                                                                                                                                                                                       ggATGCTGCTTCATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCaa            <  2:374014/251‑1 (MQ=255)                                                                                                                                                                                                                                               ctTCATCAATAA‑GGTCGATGGCTTTATCCGGCCGGTGACGCTCAGCAATGTAGCGATGAGACAACGTCGCAGCTGCAACAATGGCCAGGACCGTAATTTGCACATGGTGGTTCAATTCGTAACGTTCTGTCAGGCCCCGCAGAGTCGCAATGGTAGCTTCAACAGAAGGCTCGAAGACACACCCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGCACTGACGACATTCTTCACGcgtcctg    >  2:371521/1‑248 (MQ=255)                                                                                                                                                                                                                                                  cATCAATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc  <  2:222390/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTC‑GTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCA‑GGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGC  >  NC_000913/2732748‑2733237

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