Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	100.0%	237.2 / NA	70	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/35);  total (35/35)

CCAAGTGTGCTGATGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337151‑1337638                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ccAAGTGCGCTGATGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGCCACCGGATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGACACGGCAATCAAAAAAGGCGAttt                                                                                                                                                                                                                                               >  2:172214/1‑251 (MQ=255)     gggtgCTGATGGCCTGAGGCGTTCCTGG‑GAAAAAAGCCGTAACTCCTGGCCTAAATCCGCAACCATGGGGCAAAGCCTCGCCGGCACCGGGTTCGAAAGTTGTTTCTGAGTATCGTGCAACAGCGAAACCGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAcc                                                                                                                                                                                                                                           <  2:221631/249‑1 (MQ=255)      tgtgCTGATGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GCCCGCTGACCGATCCTATCGAAACGGAAATCAAAAAAGGCGCTTTAACCg                                                                                                                                                                                                                                          >  2:477837/1‑251 (MQ=255)          cTGATGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGTGTTCGAAAGTCGTTTCTGCTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCCCCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg                                                                                                                                                                                                                                      >  1:192193/1‑251 (MQ=255)            gATGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTTATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggt                                                                                                                                                                                                                                     <  2:570685/250‑1 (MQ=255)              tGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCATATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcg                                                                                                                                                                                                                                  >  1:719731/1‑251 (MQ=255)                 ccGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcggtg                                                                                                                                                                                                                               >  1:290698/1‑251 (MQ=255)                   gCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACTAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCt                                                                                                                                                                                                                             <  2:534686/251‑1 (MQ=255)                   gCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCt                                                                                                                                                                                                                             >  1:575911/1‑251 (MQ=255)                    cAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAACAAAGGCGATTTAACCGTCGGTGCGGTGCTg                                                                                                                                                                                                                            >  2:231741/1‑251 (MQ=255)                      ggC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTc                                                                                                                                                                                                                          <  1:733382/251‑1 (MQ=255)                       gC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTcc                                                                                                                                                                                                                         >  1:841528/1‑251 (MQ=255)                         gttGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGTGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGAGCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCg                                                                                                                                                                                                                        <  2:94705/250‑1 (MQ=255)                                                 cctGGGACTCAAGCGGCATCCATGGGTAAAAGCGTCGCTGGCACCGTGTTCGAAAGTCGTTTCTTATTATCTGGCAAAAGCGAAACTGACACCGTTTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCTTAACTTTGAAGGCCGTa                                                                                                                                                                                                <  2:575911/250‑1 (MQ=255)                                                   tGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                              >  1:594650/1‑251 (MQ=255)                                                   ggggCCGCAAGCGGCATCCATTGGTCAAAGGGTCGCAGGCACCGGGTTCGAAAGTCTTTTCTGATTATCTGGCAAAAGCTAAACTGACACCGTATCTCGACGAACTGGTGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                              <  2:269729/250‑1 (MQ=255)                                                      gCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                                           >  1:503487/1‑251 (MQ=255)                                                      gCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGAGGTGCTGTCCGGAAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                                           >  2:639109/1‑251 (MQ=255)                                                              cggaaaCCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTACCCTTGTGGGATACGGTTGTACCACCTGTCTTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGACACGGCCATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGACGTATCCATCAGCTGGt                                                                                                                                                                                   >  2:756284/5‑251 (MQ=255)                                                                 cAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTa                                                                                                                                                                                 >  1:660675/1‑250 (MQ=255)                                                                             aaaGCTTCGCTGGCACCGGGGTCGAAACTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTCTCTCGTCGAACTGGGTTTTAACCTTGTGGGATCCGGTTGTACCTCCTGTAGTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACTGCAATCAAAAAAGGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTTAAGGTCGTATCCATCCGCTGGTTAAAACTAActggct                                                                                                                                                                    <  2:636019/251‑1 (MQ=255)                                                                               aGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTTACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCTGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt                                                                                                                                                                                  <  2:627442/235‑1 (MQ=255)                                                                               aGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCTGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt                                                                                                                                                                                  >  1:627442/1‑235 (MQ=255)                                                                               aGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGGGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctgg                                                                                                                                                                  <  2:628001/251‑1 (MQ=255)                                                                                  gTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCt                                                                                                                                                               >  1:130217/1‑251 (MQ=255)                                                                                  gTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCt                                                                                                                                                               >  1:72720/1‑251 (MQ=255)                                                                                        ggCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGCCGTACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAACAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAATCTAACTGGCTGGACTcgccgc                                                                                                                                                         >  2:172450/1‑251 (MQ=255)                                                                                               ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggt                                                                                                                                                  <  1:319371/251‑1 (MQ=255)                                                                                                gTTCTAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                 <  1:393794/251‑1 (MQ=255)                                                                                                    gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                             <  1:146054/251‑1 (MQ=255)                                                                                                          cGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATg                                                                                                                                       <  1:308416/251‑1 (MQ=255)                                                                                                                tGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGACGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTgg                                                                                                                                 >  1:519321/1‑251 (MQ=255)                                                                                                                 gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTTTGGGATACGGTTGTACCACCTGTATTGGTAACCCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                                                                                                                   <  2:385544/248‑1 (MQ=255)                                                                                                                   ttATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCgg                                                                                                                              >  1:203659/1‑251 (MQ=255)                                                                                                                       cTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGACCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGACGCTGGTGGTTGCCTATGCGcgggcgggga                                                                                                                           >  2:631890/1‑248 (MQ=255)                                                                                                                        tGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGAAACTTTGAAGGCCGTATCCATCCGCTGGTTCAAACTAACTGGCTGGCCTCGCCGCCGCCGGTGGTTGCCTAAGCGCTGGCGGGaaa                                                                                                                          >  2:838904/1‑250 (MQ=255)                                                                                                                               aGCGAAAGTGACACCGTATCTCGACGAACTGGGGTTTAGCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGGGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCTTAACTTGGAAGGCCGTATCCATCCGCTGGTTAAAGCTAGCAGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAata                                                                                                                  <  1:71579/251‑1 (MQ=255)                                                                                                                                gCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAatat                                                                                                                 >  1:343947/1‑251 (MQ=255)                                                                                                                                  gAAACTGACCCCGTAGCTCGGCGAACTGGGGTTGAACCGTGTGGGATCCTGTTGTACCACCGGGATTGGTAACGCTG‑GGCCGTGGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGATGTACCCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCa                                                                                                               <  2:192193/251‑1 (MQ=255)                                                                                                                                  gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATCTGAATATc                                                                                                                >  1:785460/1‑250 (MQ=255)                                                                                                                                    aaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTCTTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGACTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATCTGAATATCAAc                                                                                                             >  2:671258/1‑251 (MQ=255)                                                                                                                                       tGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAATATGAATATCAACCTg                                                                                                          >  1:694587/1‑251 (MQ=255)                                                                                                                                           acCGTACCTCGACGAACTGGGGTTTAACCTTGTGGGTTACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGTTCCGATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCtt                                                                                                      <  2:810646/251‑1 (MQ=255)                                                                                                                                              gcactgtgacgaATTGGGGTTTAACTGTTTGGTATACGGTTGTACCCCCTGTATTGGTACCTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTg                                                                                                   <  2:531897/244‑1 (MQ=255)                                                                                                                                               tATCCCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATCGGTAACTCTG‑GGCCGCTGCCCGATCCGATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCt                                                                                                    <  2:130217/249‑1 (MQ=255)                                                                                                                                                 tctcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGc                                                                                                >  1:479388/1‑251 (MQ=255)                                                                                                                                                      acgaACTGGTGTTTAACCTTGTGGGAGTCGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGTCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACGTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGTCTTCTGAGCCTATc                                                                                           <  2:100375/251‑1 (MQ=255)                                                                                                                                                        gaactgggtgtaacccTTGTGGGATACGGTTGTACCACTTGCATTTTTAACTATG‑GGCCTCTGCCTGATCCTACCGAAACGGCAATCAAAAAAGGCGTTTGAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCgg                                                                                         <  2:742578/237‑1 (MQ=255)                                                                                                                                                        gaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTACGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGACATATGACACTCAACCTGGCTTATAAACCTATCgg                                                                                         >  1:199142/1‑251 (MQ=255)                                                                                                                                                              gggTTTAACCTTGTTGGATACGGTTGTACCCCCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGa                                                                                   <  1:671258/251‑1 (MQ=255)                                                                                                                                                              gggTTTAACATTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGAGCCGTTCGAAACGGCAATCAACAACGGCGCTTTAACCATCGGTGCGGTTCTGTCCGGCAACCGTAACTTTGATGGCAGTATCCATACGTTGGTTAAGACTAACTAGCTGGCGTCGCCACGGCTGGTGGTTGCCTATGCGCTGGCGGGTACAATTAATATCAAccccgcttgtgccgctatcgcccgac                                                                                    >  2:7642/1‑226 (MQ=255)                                                                                                                                                                gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATc                                                                                 >  1:639553/1‑251 (MQ=255)                                                                                                                                                                      ccTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCaa                                                                             >  1:828430/1‑249 (MQ=255)                                                                                                                                                                          gtgGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGTCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCg                                                                        >  1:135124/1‑250 (MQ=255)                                                                                                                                                                             ggATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATcc                                                                    <  1:369982/251‑1 (MQ=255)                                                                                                                                                                              gATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGGGATCCg                                                                   >  2:342159/1‑251 (MQ=255)                                                                                                                                                                               atacggtgggtccaccTGGATTGGGAACTCTG‑GTCCTGTGCCCGATCCCATCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCGGAGCCTATCTGCCATGATCGCAAAGGCTATCCgg                                                                  <  2:412948/240‑1 (MQ=255)                                                                                                                                                                                  cGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGttt                                                               <  1:172214/251‑1 (MQ=255)                                                                                                                                                                                         cccaccTGTAGTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGaa                                                        <  1:756284/250‑1 (MQ=255)                                                                                                                                                                                            cccggTATTGGTA‑CTCTGTGGCCGCTGCCCGATCCAATCTAATCGGAAATAAAAAAAGGCAATTTAACCTTCGGTGCGGTGTGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGa                                                     <  2:765277/247‑1 (MQ=255)                                                                                                                                                                                             cgtGTATTGGTAACTCTG‑GGCAGCTGCACGTCCCTACCGAAACGGCAATAAAAAAAGGCTATGTGACCGTCGGTGCGGTGCTGTCTGGCATCCGTAACTTTGAAGGCCGTATCCATCCGCTGGTAAAATCTACCTGGCTGGCCCCGCGGTCGCTGGGGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGCCCATGATCGCAAAGGCGAGCCGGTTTATCTGAAAGat                                                    <  1:558123/249‑1 (MQ=255)                                                                                                                                                                                               tGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCCTGATCCCAAAGTCGATCCGGGTTATCTGAAAGatat                                                  >  2:823868/1‑251 (MQ=255)                                                                                                                                                                                                              g‑ggCCGCTGCCCGATCGCAACGGAACGGTAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTTCTGTCCGGCAACCGTAACTTTGAAGGCCGTAGCCATCCGCGAGTTAAAACAAACTGGCTGGCCTCTCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCAGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCTGcaa                                    <  2:223244/250‑2 (MQ=255)                                                                                                                                                                                                              g‑ggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGca                                     <  1:477837/249‑1 (MQ=255)                                                                                                                                                                                                                      tGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCCCAAGAAAtt                            >  1:670388/1‑251 (MQ=255)                                                                                                                                                                                                                            aTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCACCCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGCTCGCCAAGGCGCTCCGGTTTTTCTGAAAAATATCCTGCCATAGGCCCCAGAACTTGCCCGt                      >  2:275239/1‑251 (MQ=255)                                                                                                                                                                                                                             tCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACGGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGt                      <  2:244292/250‑1 (MQ=255)                                                                                                                                                                                                                                   cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGGGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCGCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAg               <  2:290698/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAg         <  1:631890/251‑1 (MQ=255)                                                                                                                                                                                                                                           cAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtc       <  1:590432/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGATGGTGGTTGCCTATGCGCTGGCGGGAAATATGAAAATCAACCCGGATTCTGAGCCTATCGGCCATGATAGCCAAGGCGATCCGGTTTTTCTGAAAGATATCTTGCCATCGGCCCAAGAAAATGACCGTgaggacagacaagtctcccc  >  2:571427/1‑232 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CCAAGTGTGCTGATGGCCGCAGGC‑TTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTG‑GGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337151‑1337638

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