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

CAAGTGTGCTGATGGCCGCAGGCTTGCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCA  >  NC_000913/1337152‑1337637                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    ccagttagCTGATGGCCGCATGCGTG‑TGGACAAAAAAAACCGTAACTCTGGTCCTCA‑GCGGGCAACCATGTGGCAAAGCGTCGCTGGCACCTGGTGCGAAAGTCGATTCTGAGTATCTTGCAAAAGCGAAACTGACACCGTATCTCGATGATCTGGGTTTTACCCTTGTGGGATACGGTTGTACCACCTTTAGTGGTAACTCTGAGCCGCTGCCCGATCCTATCGAAACGGCAATCAAGAAAGGCGATTTa                                                                                                                                                                                                                                             <  2:525427/244‑1 (MQ=255) cAAGTGTGCTGATGGCCGCAGGCTTGCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                                                                                                                                                                                                                                             <  1:525996/251‑1 (MQ=255)          tGATGGCCGCAGGCTTGCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggt                                                                                                                                                                                                                                    >  1:524858/1‑251 (MQ=255)                 cGCAGGCTTGCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCggtgcggtgc                                                                                                                                                                                                                             >  1:215860/1‑251 (MQ=255)                  ccAGGCTTTCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTAGCTGGCAAAAGCGAAACTGCCACCGGATCTCGACGAACTGGGGTTTAACCTTTTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGAGCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGTGGTGCt                                                                                                                                                                                                                            <  2:421543/250‑1 (MQ=255)                          gCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGc                                                                                                                                                                                                                     >  1:385899/1‑250 (MQ=255)                           cTGG‑CGAAAAAAGCAGGAAATCTGGGAAACAAGGGG‑AAAAGATGGGTCAAAGCGACGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCa                                                                                                                                                                                                                    >  1:90423/1‑250 (MQ=255)                                cGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGGCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGGCCGGCAACCGt                                                                                                                                                                                                               <  2:590159/251‑1 (MQ=255)                                  aaaaaaGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTa                                                                                                                                                                                                              >  1:104401/1‑250 (MQ=255)                                   aaaaaGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGGGGGATACGGTTGTACCTCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAAc                                                                                                                                                                                                            <  1:588277/251‑1 (MQ=255)                                     aaaGCCGTAACTCTGGGCCTCAATCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAAGCAAAAAAGGCGATTTaa                                                                                                                                                                                                                                            <  2:4538/217‑1 (MQ=255)                                     aaaGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTaa                                                                                                                                                                                                                                            >  1:4538/1‑217 (MQ=255)                                       aGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTg                                                                                                                                                                                                        <  1:52916/251‑1 (MQ=255)                                           gTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAgg                                                                                                                                                                                                    >  1:128192/1‑251 (MQ=255)                                             aaCTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGcc                                                                                                                                                                                                  <  1:129922/251‑1 (MQ=255)                                              cccctGGGCTTCAAGCAG‑CAACCATGGGTCAAAGCGACGCTGGCACCGGGTTCGAAAGCCGTGTCTGAGTATCTGGCAAAAGCGAAACTTACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATGCGGTTGTACCACCGTTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCg                                                                                                                                                                                                 <  2:510310/248‑1 (MQ=255)                                              aCTCTGAGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCAACGGGTTCGAACGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTGCCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGAGATTTAATCGTCGGTGCGGTGATGTCCGGCACCCGTAACTTTGAAGGCCg                                                                                                                                                                                                 >  2:321728/1‑251 (MQ=255)                                               ctctGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATATCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCTATCAAAAAAGGCGATTTAAACGTCGGTGCGGTGCTGTCCGGCAACAGTAACTTTGAAGGCCGt                                                                                                                                                                                                >  2:172364/1‑251 (MQ=255)                                                   gggCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTat                                                                                                                                                                                              <  1:172364/249‑1 (MQ=255)                                                           aaacggcagccttgtggcAAAGCGCCGTTGGCACCGGGTTCGAAAGGCTGTTCTGATTATTTGGCGAAAGCGTAACTGACTCCGGATCTCGACTAACTGGGGTGTAACGTTGGGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCAGCTGCCCTATCCTATCGAACCTGCAATCAAAAAAGGGTATTTAACCGTCTGTGCGGTGCTGTCCGGCAACCGTAACTTTGAATGCCGTATCCATCCGca                                                                                                                                                                                     <  2:201778/234‑2 (MQ=255)                                                             cGG‑CAACCATGGGTAAAAGCGTCGCTGGCACCGGGTTCGACAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAGCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTCTCGAACCGGCATTCAAAAAAGGCGATTTAACCGTCGGGGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCCGATCCATCCGCTggg                                                                                                                                                                                  >  2:303020/1‑250 (MQ=255)                                                                    ccATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt                                                                                                                                                                                                                                                >  2:404732/1‑183 (MQ=255)                                                                    ccATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt                                                                                                                                                                                                                                                <  1:404732/183‑1 (MQ=255)                                                                     cATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAAATTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                           >  1:58879/1‑251 (MQ=255)                                                                      aTGGGTCAAAGCGTCGCTGGCAACGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                           >  2:163709/1‑250 (MQ=255)                                                                          gTCAAAGCGTCGCTGGCCCCGGGTTCGAAAGTCGTTGTTTATTATCTGGCAAAAGAGAAACTGACTCCGTCTCTCGACGCACTGGGGTTTAACCTTGTGTGATCCGTTTGTCCCACCTGTATTGGTAACTCTGGGCCGCGGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAATTAACCGTCGTTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAActg                                                                                                                                                                      <  2:269360/251‑1 (MQ=255)                                                                                   tCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                             >  1:403226/1‑219 (MQ=255)                                                                                   tCGCTGGCACCGGGTTCGAAAGTCGTTTCTGAGTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGATCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATGTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                             <  2:403226/219‑1 (MQ=255)                                                                                    cgccggCACCGGGTTCGATAGTCGTTTCTGATTATCTGTAAAAAGCGAAACTGACACCGTATCTCGACGCTCTGGGGTTTAACCTTGGGGGTTACGGTTGTACCACCAGTATGGGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTACACGTCGGTGCGGTGCTGTCCGGCACCCGTAACTTTTAAGGCCGTCACCATCCGCTGGTTAAAACGAACTGGCGGGccgc                                                                                                                                                             <  2:90423/246‑3 (MQ=255)                                                                                    cGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcg                                                                                                                                                            <  2:104401/251‑1 (MQ=255)                                                                                         gCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTAGTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgcc                                                                                                                                                       <  1:459441/251‑1 (MQ=255)                                                                                                gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTTTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                <  1:11083/251‑1 (MQ=255)                                                                                                    gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                            >  1:372111/1‑251 (MQ=255)                                                                                                      aaGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGcc                                                                                                                                          <  1:486279/251‑1 (MQ=255)                                                                                                           gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgc                                                                                                                                     <  1:437350/251‑1 (MQ=255)                                                                                                             ttCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAAACTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTACGGTGCTGTCCGGCAACCGTAACTTTGAGGGCCGTGTCCATCCGCTGGTTACAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcgc                                                                                                                                   >  1:540252/1‑251 (MQ=255)                                                                                                                              aaGCGAAACTGACACCGTATCTCGACGAAATGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                  >  1:463142/1‑251 (MQ=255)                                                                                                                                             cGTATCACGACGAACTGGGGTTTAACCTTGTGGGCTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCACTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGgaatatgaatatacaccgtgcttct                                                                                                   >  2:98970/1‑238 (MQ=255)                                                                                                                                                    cgacgaACTGGGGTATAACCTTGTGGGAGACGGTTGTACACCCTGTATTGGTAACTTTGGGCCGCTGCCCTATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACAGTAACTGTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGc                                                                                               <  2:130443/248‑1 (MQ=255)                                                                                                                                                               ggattAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGAGATTTAACCGTCGGTGCGGAGCTGTACGGCACCAGTGAATTTGAAGGCCGTATCCATCCGCTGTTTACAACTAACTGGATGGCCTCGCCGCCGCTGGTGGGTGCCTCTGCGCTGGCGGGAAATTAGAATATCTACCGGGCATCTGAggcttacggccctgt                                                                                  >  2:666843/4‑236 (MQ=255)                                                                                                                                                                gTTTAACCTTTGGGGAAACGGTTGTACCACCCGTATTGGTAACTCTGGGCCGCGGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCAATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATACCAACCTGGCTTCTGAGCCTATCGGCCATGATc                                                                                <  1:226326/251‑1 (MQ=255)                                                                                                                                                                   cacccTTGTGGGATGCGGTTGTACCATCTGTAATGGTAACTCTGGGCCGCTGCCCGAGCCTATCGAAACGGCAATCAAAAAAGTCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAGCCGGAACTTTGAAGGCCGTATCCCTCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGc                                                                              <  2:200115/247‑1 (MQ=255)                                                                                                                                                                          gtgGGATACGGTTGTACCTCCTTTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAGAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTGGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGa                                                                      <  2:590021/251‑1 (MQ=255)                                                                                                                                                                           ggggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAt                                                                     <  1:83084/250‑1 (MQ=255)                                                                                                                                                                             ttccgttcttttACCACCTGTATTGGTAACTCTGGGCCCCTGCCCGATCCTATCGAAACGGCAAACAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATc                                                                    <  2:373108/239‑1 (MQ=255)                                                                                                                                                                              gattcgggtGTACCACCTGTATTGGTAACTTTGGGCCGCTGCCGGACCCTAACGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGTCGATCCg                                                                  <  2:670293/243‑1 (MQ=255)                                                                                                                                                                                      tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAACAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGAAACTTTGCAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATCTGAATATCAACCTGGCTTCTGAGCCTATCGGCCAGGATCGCAAAAGCGATCCGGTTTATCt                                                          >  1:325651/1‑251 (MQ=255)                                                                                                                                                                                      tGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATc                                                                    >  1:373108/1‑241 (MQ=255)                                                                                                                                                                                       gTACCACCTGTATTGGTAACTCTGGGCCGATGCCCGCTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATCCGCTGGCGGAGAATATGAAAATCTACCTGGCGTCTGAGCCTAAAGGCCATGATCGCAAAGGAGATCAGGtgtttctg                                                         >  2:183683/1‑244 (MQ=255)                                                                                                                                                                                       gTACCACCTGTATGGGTAACTCGGGGCCGCTGCCTGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGATCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCt                                                          <  2:361824/250‑1 (MQ=255)                                                                                                                                                                                          ccaccTGTAGTGGTAACTCTGGGCCGCTGCCATAGCCGAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGCCCGGCAACCGTAACTTTGAAGGCCGTATCCATCAGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGaaa                                                      <  2:163068/251‑1 (MQ=255)                                                                                                                                                                                          ccaccGGTATTGTTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGaaa                                                      <  2:33735/251‑1 (MQ=255)                                                                                                                                                                                           caccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAg                                                     <  1:145809/251‑1 (MQ=255)                                                                                                                                                                                                    tGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGc                                            >  1:275611/1‑251 (MQ=255)                                                                                                                                                                                                                                aTCGAAACTGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcg                                                                                                                                    <  2:510927/135‑1 (MQ=255)                                                                                                                                                                                                                                aTCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcg                                                                                                                                    >  1:510927/1‑135 (MQ=255)                                                                                                                                                                                                                                  cGAAACGGCAATCAGAAACGGCGCTTTAACCGTCGGTGCGGAGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcg                                                                                                                                                            >  2:428328/1‑109 (MQ=255)                                                                                                                                                                                                                                  cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGTCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcg                                                                                                                                                            <  1:428328/109‑1 (MQ=255)                                                                                                                                                                                                                                    aaaCGGCAATCAAAAAAGGCGATTTATCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGaa            <  1:423152/251‑1 (MQ=255)                                                                                                                                                                                                                                     aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCCAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAAc           >  1:215228/1‑251 (MQ=255)                                                                                                                                                                                                                                     aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGACTATGCGCTGGCGGGAAATATGAATATAAACCTGGCTTCTGAACCTATCGGCCATGATCGGAAATGCGATCCGATTTTTTTCAAAGAGATTAGGCCAACGGAACAAGAAATTGCCCgggagggagaac           >  2:368096/1‑240 (MQ=255)                                                                                                                                                                                                                                       cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATACGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACACGAAATTGCCCGTGCGGTAGAACa          >  1:464562/1‑250 (MQ=255)                                                                                                                                                                                                                                              aaaaaaagGCGATTTAACCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCACGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCca  <  2:215228/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    CAAGTGTGCTGATGGCCGCAGGCTTGCTGG‑CGAAAAAAGCCGTAACTCTGGGCCTCAAGCGG‑CAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCA  >  NC_000913/1337152‑1337637

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