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%	174.6 / NA	56	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/34);  total (22/34)

CCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337151‑1337638                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ccAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAttt                                                                                                                                                                                                                                                >  1:35661/1‑251 (MQ=255)   aaGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                                                                                                                                                                                                                                               >  1:168263/1‑250 (MQ=255)   aaGTGTGCGGATGGCCGCAGGCTTGCGGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCCACCATGGTTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCGCGACGAACTGGGGTTTAACCTTGTGGGATACTGGTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGTCGATTTaa                                                                                                                                                                                                                                              <  2:614339/251‑1 (MQ=255)         gatgaggccgtcAGGCTTGCGTGAGAAAAAAGCCGTAACTCTGTGCCGCAAGCGGCAACAATTGGTGAAA‑GTGTCGCGGGCACCGGGTTTGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCTACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGTGCCGCTGCCCGATCCGATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                                                                         <  2:495343/239‑1 (MQ=255)               ggCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg‑tgcg                                                                                                                                                                                                                                  <  2:273818/250‑1 (MQ=255)                              ggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTCACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGCTCCTATCGAAACGGCAATCAAAAAAGGAGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCaccc                                                                                                                                                                                                                  >  2:310699/1‑248 (MQ=255)                               gagAAAAAATCTGTAACCCTGGGCCCCAAGCGGCAACCATGGGTCAAA‑GCGTCCCTGGCACCGGGTTCGAAAGTCGTTTCTTATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGATCTGGGGTGTAACCTTGTGGGATACGGTTGGACCACCTGTATTGGTAACTCTGGGCCGCTGACCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCATCCg                                                                                                                                                                                                                 <  2:404626/249‑1 (MQ=255)                                  aaaaaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTa                                                                                                                                                                                                               >  1:697177/1‑250 (MQ=255)                                     aaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACtt                                                                                                                                                                                                           <  1:151365/251‑1 (MQ=255)                                      aaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACttt                                                                                                                                                                                                          <  1:137429/251‑1 (MQ=255)                                              aCTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑GGCGGTGCTGTCCGGCAACCGAAACTTTGAAGGCCg                                                                                                                                                                                                  >  2:512651/1‑251 (MQ=255)                                                      ccTCAAGCGGCAGCCATGGGTCAAATGGGT‑GCTGGCACCGGGTTCGAAAGTCGTTGCTGATTATCTGGCAAAAGCGAAACTTACTCCGTATCCCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCAGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatc                                                                                                                                                                                          <  2:154779/251‑1 (MQ=255)                                                                                        ggCACCGGGTTCGAAAGTCGGTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                                                                                                                           <  2:224835/249‑1 (MQ=255)                                                                                          cACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgccg                                                                                                                                                       <  1:227329/251‑1 (MQ=255)                                                                                          cACCGGGGTCGAAAGTCGTTTCTGCTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAAATTTGAAGGCCGGATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgccg                                                                                                                                                       >  2:425935/1‑251 (MQ=255)                                                                                             cGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtg                                                                                                                                                    >  1:529731/1‑251 (MQ=255)                                                                                             cGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTACGGCAACCGTAACTTTGACGGCCGTATCCATCTGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCTCtg                                                                                                                                                    >  2:261304/1‑251 (MQ=255)                                                                                                   cGAAAGTCGTTTCTGATTACCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTggtggt                                                                                                                                               <  1:397003/250‑1 (MQ=255)                                                                                                           gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGACACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                            >  2:109874/1‑213 (MQ=255)                                                                                                           gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                            <  1:109874/213‑1 (MQ=255)                                                                                                                 gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTgg                                                                                                                                 <  1:512651/250‑1 (MQ=255)                                                                                                                         ggCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGaaa                                                                                                                          <  1:708999/249‑1 (MQ=255)                                                                                                                                   aaaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgccgc                                                                                                                                                      <  1:658110/211‑1 (MQ=255)                                                                                                                                   aaaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCa                                                                                                               >  1:230964/1‑250 (MQ=255)                                                                                                                                        gACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCt                                                                                                           <  1:105449/249‑1 (MQ=255)                                                                                                                                          cacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggt                                                                                                                                                  <  2:9695/208‑1 (MQ=255)                                                                                                                                          cacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggt                                                                                                                                                  >  1:9695/1‑208 (MQ=255)                                                                                                                                            cctaaTCTCGACGAATTGGGTTTTAACCTTGTGGGATACGGTTGTACCACCTGTATGGGTAACTGTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTACCCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCt                                                                                                       <  2:381131/245‑1 (MQ=255)                                                                                                                                                        gaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAATATGAATTTCAAACTgggatcgggagagcgtcg                                                                                         >  2:60836/1‑235 (MQ=255)                                                                                                                                                        gaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTg                                                                                                          <  1:60836/234‑1 (MQ=255)                                                                                                                                                          aCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGcc                                                                                       >  1:441986/1‑251 (MQ=255)                                                                                                                                                          aCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATCAATATCAACCTGGCTTCTGAGCCTCTCGGcc                                                                                       >  2:328494/1‑251 (MQ=255)                                                                                                                                                                gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGACATGATc                                                                                 >  1:94479/1‑251 (MQ=255)                                                                                                                                                                      cctggtgGAATACGGTTGTACCACCTGTATTGGTAACTTTGGGCCTCTTCCCGATCCTACCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TTCGGTGCGGACCGGAACCCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAg                                                                           <  2:512305/247‑1 (MQ=255)                                                                                                                                                                      ccTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAg                                                                           <  1:265886/251‑1 (MQ=255)                                                                                                                                                                         tgtgGGATACGGTTGTACCCCCTGTAGTGGTAACTCGGTGACGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAgg                                                                          <  2:169044/249‑1 (MQ=255)                                                                                                                                                                         tgtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAgg                                                                          <  1:661774/249‑1 (MQ=255)                                                                                                                                                                                  ccgttttaccctcttgtgtGGTAACTCTGGGCAGCTGCCTGATCCTATGTAAAGTGCAATCAAAAAAGGCGAGTTAGCCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAATGCCTTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCTGTAAATATGAATACCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGtgt                                                               <  2:180837/233‑3 (MQ=255)                                                                                                                                                                                    gTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTAt                                                             <  1:492291/251‑1 (MQ=255)                                                                                                                                                                                       gTACCACCGGTATTGGTAACTCTGGGCCTCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCt                                                           <  2:140895/250‑1 (MQ=255)                                                                                                                                                                                        tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTg                                                          >  1:4913/1‑250 (MQ=255)                                                                                                                                                                                                    tGGTAACTCGTGGACGCTGCCCGATCCTAACGAAACGACAATCAAAAAAGGCGAGTTAACCGGCGG‑TGCGGTGCTGTCCGGCAACCGTACCTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGATGGCCTCGCAGCCGCTGGTGGTTGCCTATGCGCTGGCGTGAAATATGAATATCAACCTGGCTTTTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTgg                                              <  2:689258/250‑1 (MQ=255)                                                                                                                                                                                                           cctGGGCCGCTGCCCGATCCTACCGAAACGGCAATCAAAAAAGGCGATGTAACCGTCGG‑TGCGGTGCTGTCCTGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGTCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGc                                      <  1:297468/250‑1 (MQ=255)                                                                                                                                                                                                              gggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAACTATGAATATCAACCCGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGCTATCTGGCCATCGGcac                                    >  1:621469/1‑250 (MQ=255)                                                                                                                                                                                                                  cGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGaa                               <  1:76842/251‑1 (MQ=255)                                                                                                                                                                                                                   gCTGCCCGTTCCTATCGAAACGGCAATCAAAAAAGGCGAGTGAACCGTCGG‑TGCGGTGCTGTCCGTCAACCGTAACGTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGAACGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGaaa                              <  2:42887/251‑1 (MQ=255)                                                                                                                                                                                                                      gCCCGATCCTATCGAAACGGCAAACAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTg                           <  1:15340/251‑1 (MQ=255)                                                                                                                                                                                                                        ccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCAATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATATGGCCATCGGCACAAGAAATTggc                         >  1:149592/1‑249 (MQ=255)                                                                                                                                                                                                                          gATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCg                       >  1:70794/1‑251 (MQ=255)                                                                                                                                                                                                                               tagcGAAACGGCAATCAAAAAAGGAGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCg                   <  2:230964/247‑1 (MQ=255)                                                                                                                                                                                                                                aTCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCgg                  <  1:433249/250‑1 (MQ=255)                                                                                                                                                                                                                                   gAAACGGCAATCAAAAAAGGCGATTTACCCGTCGG‑TGCGGTGCTGTCCGGCATCCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGa              <  2:441986/251‑1 (MQ=255)                                                                                                                                                                                                                                   gAAACGGAAATCAAAAAAGTCGTTTAAACCGACGGGTG‑GGTGCTGTCCGGCATCCGGAACTTTGAGGGCCGTATCCATCCGCGGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAAGATGAATATCATCCTGGCTTCTGATCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTAGCTGAAAGATATCTGGCCAGCGGCACAAGAAATTGCCCTTGCGGTAGa              <  2:32690/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggCAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAAGCCGTATCCATCCCCTGGTTCCAACTAACTGGCTGGACTCGCCGCCGCTGGTAGTTGCCTATGCGCTGGCGGGGAATATGAATAACAACCTGGCTTCTGCGCCTACCGGCCATGATCGCAACGGCGACCCGGTTTCTCTTAAAGATATCTGGCCATCCGCCCCAGCACTTGCCCCTGCGGTCGAACaa          >  1:658932/1‑250 (MQ=255)                                                                                                                                                                                                                                          cAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtc       >  1:472883/1‑251 (MQ=255)                                                                                                                                                                                                                                          cAATCAAAAAAGGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGGAGAACAAGtc       >  1:201746/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaaagGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCcac  >  1:78334/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAA‑GCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGG‑TGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  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: ‑