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

GCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337159‑1337638                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      gCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                                                                         <  2:429408/250‑1 (MQ=255)   tGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGGTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggg                                                                                                                                                                                                                                      >  2:265050/1‑250 (MQ=255)    gATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggt                                                                                                                                                                                                                                      <  1:42972/250‑1 (MQ=255)              ggCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAATCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTc                                                                                                                                                                                                                           >  1:129259/1‑251 (MQ=255)                               aGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGCTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACAGTCGGTGCGGAGCTGTCCGGCAACCGTAACTTTg                                                                                                                                                                                                          >  2:418122/1‑251 (MQ=255)                                               cTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatcc                                                                                                                                                                                          >  1:269876/1‑251 (MQ=255)                                                 cAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACGCCGTATCTCGACGAACCGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatcc                                                                                                                                                                                          <  2:320090/249‑1 (MQ=255)                                                      ggCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt                                                                                                                                                                                   >  1:288794/1‑251 (MQ=255)                                                       gCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAGGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTa                                                                                                                                                                                  <  1:500904/251‑1 (MQ=255)                                                          aCCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaaaa                                                                                                                                                                               <  1:52869/251‑1 (MQ=255)                                                                      aGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctgg                                                                                                                                                                   >  1:15805/1‑251 (MQ=255)                                                                      aGCGTCGCTGGCACCGGGTTCCAACGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGACCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctgg                                                                                                                                                                   >  2:399307/1‑251 (MQ=255)                                                                           cGCTGGCACCGGGTTCGAAAGTCGGTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTAGCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcg                                                                                                                                                              <  2:286729/251‑1 (MQ=255)                                                                                       gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTCTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                  >  1:379338/1‑251 (MQ=255)                                                                                           gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                              >  2:112242/1‑251 (MQ=255)                                                                                             aaGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGcc                                                                                                                                            <  1:279593/251‑1 (MQ=255)                                                                                                  gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgc                                                                                                                                       <  2:433056/251‑1 (MQ=255)                                                                                                        gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                                                                                                                    <  1:152811/248‑1 (MQ=255)                                                                                                         gttATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTATCCTTGTGGGATACGGTTTTACCACC‑TGTATTGGTAACTCGGGGCTGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGTGGATTTAACCGTCGGTGCGGTGCTGGCCGTCAACCGTAACTGTGAAGTCCGTAGCCAGCCGCTGGTTAAAGCTAACTGGCTGGACTCGCCACCGCTGGTGGTTGCCTATGTTCTGGCg                                                                                                                                <  2:193129/250‑1 (MQ=255)                                                                                                         aTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGACCTGGGGTTTAACCTTCTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCg                                                                                                                                >  2:114229/1‑251 (MQ=255)                                                                                                                         gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATc                                                                                                                 >  1:329839/1‑250 (MQ=255)                                                                                                                          aaaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCa                                                                                                                >  1:361817/1‑250 (MQ=255)                                                                                                                            aCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAcc                                                                                                             >  1:317156/1‑251 (MQ=255)                                                                                                                            aCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAcc                                                                                                             >  1:80821/1‑251 (MQ=255)                                                                                                                                   ccGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAATATGAATAACAACCTGGCTTc                                                                                                      >  2:361813/1‑251 (MQ=255)                                                                                                                                        tctcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCAGGCTTCTTAGc                                                                                                 >  2:484227/1‑251 (MQ=255)                                                                                                                                           cgtctaactttagTTTAACCGTGTGGGATACTTTTGTACCACC‑TGTATTGGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTGAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTATCTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGccg                                                                                               <  2:77418/238‑2 (MQ=255)                                                                                                                                             atgaACTGGGGTCTAACCTTGTGGGATACGCTT‑TAACACCCTGTATTGGTAACTCTGGGCCGCTGCCCGTTCCTATCGACACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTTTACGGCATCCGTAACTTTGAAGGCCGTATCCCTCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGTTTGCCTATGCGCTGGCGGGAAATATGAATATCACCCTGGCTTCTGAGCCTATc                                                                                            <  1:410115/249‑1 (MQ=255)                                                                                                                                                 aCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGcc                                                                                        >  1:207063/1‑251 (MQ=255)                                                                                                                                                       gTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATc                                                                                  <  1:415911/251‑1 (MQ=255)                                                                                                                                                               ttGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGc                                                                          >  1:155048/1‑251 (MQ=255)                                                                                                                                                               ttGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGc                                                                          >  1:274307/1‑251 (MQ=255)                                                                                                                                                                 gtgGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCg                                                                         >  1:356491/1‑250 (MQ=255)                                                                                                                                                                        aCGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGtt                                                                 <  1:343059/251‑1 (MQ=255)                                                                                                                                                                        aCGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGtt                                                                 >  1:143562/1‑251 (MQ=255)                                                                                                                                                                             tGTACCCCC‑TTTATTGGTAACTATGGGCCGCTGCCCGTGCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTTGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCTGCCATGATCTCAAAGGCGATCCGGTTTATc                                                             <  2:164580/250‑1 (MQ=255)                                                                                                                                                                             tGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCTATTTACCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGGGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATc                                                             <  2:183687/250‑1 (MQ=255)                                                                                                                                                                                    cc‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGat                                                     <  2:517975/251‑1 (MQ=255)                                                                                                                                                                                      ctctcttGGCAACTCTGGGCCGCTGCCCGTTCTTATCGAAACGTCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCGGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCAGGCTTCTGAGCCTATCGGCCATGATGGCAAAGGCGATCCGGTTTATCTGAAAGat                                                     >  2:21992/6‑250 (MQ=255)                                                                                                                                                                                                      gggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGca                                     <  2:263472/249‑1 (MQ=255)                                                                                                                                                                                                       ggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGCAGGCCGTCTCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGcaca                                   >  1:193883/1‑250 (MQ=255)                                                                                                                                                                                                        gCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAACGGCGATTTAACCGTCGGTGCGGAGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGCCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGCAATAGGAATATCAAACTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAGCCGGTTTATCTCAAAGATATCTGG‑CAATCGGcaca                                   >  1:389984/1‑249 (MQ=255)                                                                                                                                                                                                        gCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAg                                 <  2:269876/251‑1 (MQ=255)                                                                                                                                                                                                           gCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGaa                               >  1:109901/1‑250 (MQ=255)                                                                                                                                                                                                                 cGATCCTATCGAAACGGCAATCCAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCCAAGGCGATCCGGTTTATCTGCAAGATATCTGG‑CCATCGGCACAAGACATTGccc                        >  2:37927/1‑251 (MQ=255)                                                                                                                                                                                                                   aTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGt                      >  1:505374/1‑251 (MQ=255)                                                                                                                                                                                                                     ccTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGc                    >  1:41387/1‑251 (MQ=255)                                                                                                                                                                                                                        aTCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATATGG‑CCATCGGCACAAGAAATTGCCCGTGCGGt                 >  1:177178/1‑251 (MQ=255)                                                                                                                                                                                                                          cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGGTCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGCggg                 >  2:252711/1‑248 (MQ=255)                                                                                                                                                                                                                           gAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGTTGCTGTCCGGCAACCGTAACTTTGAATGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGCAATATTAATATCAACCTTGCTTCTGAACCTATCGGCCATGATCACAAAGGCGATCCGGTTTATCTGCAAGATATCTGG‑CCATCTGCACAATCAATTGCCCGTTCGGttg               >  1:506896/1‑248 (MQ=255)                                                                                                                                                                                                                           gAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCCACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGCTCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGTGGt                 >  1:195275/1‑248 (MQ=255)                                                                                                                                                                                                                            aaaCGGCAATCAAAAAAGGCGATTTAACCATCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTACAACTCACTGGCTGGCCTTGCCGCCGCTGGTGGTTGGATATGCACTGGCGGGAAAAACGAATAACAACCTGGCATAGGAGCCCATTGGCCAAGAACCCACAGGCGATCCTGTTTATCTGCAAG‑TTTCTGGGCCATCTGCACCAGAAATGGTCCGTGCgtgaga              >  2:424622/1‑245 (MQ=255)                                                                                                                                                                                                                              aCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGCGGTAGAACa           <  1:96726/251‑1 (MQ=255)                                                                                                                                                                                                                                   aaTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtct      >  1:390278/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGGTTATCTGAAAGATATCTGG‑CCATCGGCCCAAGAAATTGCCCGTGCGGTAGAACAAGTCTccca  >  2:146995/1‑248 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      GCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGG‑CCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCAC  >  NC_000913/1337159‑1337638

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