Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,217,310	A→G	S68G (AGC→GGC)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,217,310	0	A	G	100.0%	256.9 / NA	69	S68G (AGC→GGC)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/33);  total (36/33)

ACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATT  >  NC_000913/4217063‑4217556                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        aCCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGcc                                                                                                                                                                                                                                                      >  2:377163/1‑251 (MQ=255)  ccGCCTGTTAGAGTCAACCCCCCCATAACTATGGAGCATCTGCACATGACAACCCGTACCCAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCt                                                                                                                                                                                                                                                     >  2:534487/1‑251 (MQ=255)   cGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTc                                                                                                                                                                                                                                                    >  2:468236/1‑251 (MQ=255)      cTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGgc                                                                                                                                                                                                                                                 >  1:275606/1‑251 (MQ=255)      cTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGgc                                                                                                                                                                                                                                                 >  2:384729/1‑251 (MQ=255)                 aaCCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTgg                                                                                                                                                                                                                                      >  2:575851/1‑251 (MQ=255)                   ccaccaCATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCg                                                                                                                                                                                                                                    <  2:176317/251‑1 (MQ=255)                   ccaccaCATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCg                                                                                                                                                                                                                                    >  1:481073/1‑251 (MQ=255)                    caccacATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCgg                                                                                                                                                                                                                                   >  2:55681/1‑251 (MQ=255)                        acaTAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAg                                                                                                                                                                                                                               <  2:53648/251‑1 (MQ=255)                        acaTAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAg                                                                                                                                                                                                                               <  1:408624/251‑1 (MQ=255)                         caTAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAgg                                                                                                                                                                                                                              <  1:521050/251‑1 (MQ=255)                         caTAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAgg                                                                                                                                                                                                                              >  1:230095/1‑251 (MQ=255)                            aaCTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGgcgc                                                                                                                                                                                                                           >  1:320537/1‑251 (MQ=255)                            aaCTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGgcgc                                                                                                                                                                                                                           >  2:458927/1‑251 (MQ=255)                                aTGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑a                                                                                                                                                                                                                      <  2:425577/251‑1 (MQ=255)                                  ggAGCATCTGCACATGAAAACCCGTACACAACACATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑aa                                                                                                                                                                                                                     >  2:466142/1‑250 (MQ=255)                                   gAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AAc                                                                                                                                                                                                                    >  1:30241/1‑250 (MQ=255)                                           gCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGaa                                                                                                                                                                                                             <  1:464965/249‑1 (MQ=255)                                            cacaTGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGATACATCAACGGCCTCGGCGCACTGACTGGCGGTCCGGCGCTGC‑AACAGGCGaaaa                                                                                                                                                                                                           >  1:563337/1‑249 (MQ=255)                                                                    aTTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTc                                                                                                                                                                                  <  1:169453/251‑1 (MQ=255)                                                                       gaagaaTTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCggg                                                                                                                                                                               >  1:93754/1‑251 (MQ=255)                                                                         agaaTTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGAt                                                                                                                                                                             <  1:54076/251‑1 (MQ=255)                                                                            aTTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGc                                                                                                                                                                          <  1:326749/251‑1 (MQ=255)                                                                            aTTACAGAAAGAGTGGACGCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCGATCTGTCGGGATGGc                                                                                                                                                                          <  2:211033/251‑1 (MQ=255)                                                                                 aGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGGATTGAAGCAGTCTATCTGTCGGGATGGCAgggg                                                                                                                                                                     >  2:553352/1‑249 (MQ=255)                                                                                  gAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAg                                                                                                                                                                    <  2:530853/251‑1 (MQ=255)                                                                                    aaGAGTGGACTCAACCGCGTTGGGAACGCATTACTCGCACCTACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAACTGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCCGGCGCTGCCAAC‑GGCGAAAGCGGGTCTTGAAGCAGTTTATCTCTCGGGATGGCAGGtcacg                                                                                                                                                                  >  2:186700/1‑247 (MQ=255)                                                                                     agagGGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACGGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAgcgg                                                                                                                                                                 <  1:237776/251‑1 (MQ=255)                                                                                           gACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGCAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGATGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑CACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTAGGGAAGGCAGGAAGCGGCGGACg                                                                                                                                                           >  1:50878/1‑251 (MQ=255)                                                                                            aCTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGCATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAACGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑ACCAGGCGAAAACGGGGCTTGAAGCAGTCTATCTGTCGGGCTGGCCGATAGcggcggcacc                                                                                                                                                          >  2:585287/1‑247 (MQ=255)                                                                                                       ttGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAACGGGTATTGAAGCAGTCTATCTGTCCGGCTGGCAGGGAGcggcgggccccgacccggcgg                                                                                                                                               >  2:569319/1‑236 (MQ=255)                                                                                                         gggAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCAGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGGAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAACGGCGGACGCTAACATTGCGGcc                                                                                                                                             >  1:191237/1‑251 (MQ=255)                                                                                                           gAAGGCATTACTCGCACATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGGAGCGGCGGACGCGAACCCGGCGGcc                                                                                                                                             >  1:165338/1‑249 (MQ=255)                                                                                                            aaGGCATTCCTCGCCCATACAGTGCGGAAGATGTTTTGAAATCACGCGGGTCAGTCAATCCTGAATGCGCGCTGTCGCAACTGTGGGCAGCGAAAATGTGGCTTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGTTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTATCGGCGGACGCTAACCTGGCGGCCAGc                                                                                                                                          <  1:275646/251‑1 (MQ=255)                                                                                                             aGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGAAGCGGCGGACGCTAAACTGGCGGcccgcc                                                                                                                                         >  2:548580/1‑247 (MQ=255)                                                                                                                      tCGCCCATACAGAGCGGTAGATGTGGCGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGAGAAACTGAGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGGAGCGGGGGACGCTAAACCGGGGGCCCGCATGTATCCgg                                                                                                                                >  2:266967/1‑251 (MQ=255)                                                                                                                            aTACAGTGCGGAAGATGTGGTGAAATTGCGCGGTTCAGTCAATCCTGAATGCCCGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAATCGGGTATTGAAGCAGTCTATCTGGCGGGATGGTAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGt                                                                                                                          <  2:275606/251‑1 (MQ=255)                                                                                                                              aCAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCCGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCg                                                                                                                        >  1:47576/1‑251 (MQ=255)                                                                                                                                aGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACACCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCAAGCATGTATCCGGATCAGTCGct                                                                                                                      <  2:195809/251‑1 (MQ=255)                                                                                                                                   gCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGctct                                                                                                                    <  2:251057/250‑1 (MQ=255)                                                                                                                                         gATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCGGCATGTATCCGGATCAGTCGCTCTATCCGGc                                                                                                             <  2:67956/251‑1 (MQ=255)                                                                                                                                           tgtgGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCa                                                                                                            >  1:289638/1‑250 (MQ=255)                                                                                                                                            gtggtgAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCaa                                                                                                           <  1:575851/250‑1 (MQ=255)                                                                                                                                                    aTTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGc                                                                                                  <  1:369903/251‑1 (MQ=255)                                                                                                                                                       aCGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAg                                                                                               <  2:430053/251‑1 (MQ=255)                                                                                                                                                            gTTCAGTCAATCCTGCATGCACGCTGGCGCAACTGGTCGCCGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACACGCGAAAGCGGGTATTGAAGCAGTCCATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCCGGCGGCCAGCCTGTATCCGGGTTAGTCGGTCTATCCGGCAAACTCGGGGCCCGCtggg                                                                                          >  2:502037/1‑249 (MQ=255)                                                                                                                                                                aGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGGAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCCGCATGTTTCCGGATCAGGCGCCCTATCCGGCAAAATCGGTGCCAGCTGtggtgg                                                                                      >  2:535538/1‑251 (MQ=255)                                                                                                                                                                aGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGtggtgg                                                                                      <  1:560863/251‑1 (MQ=255)                                                                                                                                                                    aaTCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCt                                                                                                                                                    >  1:522824/1‑185 (MQ=255)                                                                                                                                                                    aaTCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCt                                                                                                                                                    <  2:522824/185‑1 (MQ=255)                                                                                                                                                                       ccTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGGCTATCTGTCGGGATGGCAGGGAGCGGCGGACGCTCACCCGGCGGCCAGCATGGATCCGGGTCAGTCGCTCTATCCGGCAAAATCGGGGCCAGCTGGGGGGGAgggggg                                                                               >  2:223671/1‑246 (MQ=255)                                                                                                                                                                           aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATcaac                                                                           <  1:59662/251‑1 (MQ=255)                                                                                                                                                                               cACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAacaaca                                                                       <  2:68217/251‑1 (MQ=255)                                                                                                                                                                                aCGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGGCGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCCCTATTCGGCCAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAacac                                                                      >  2:41299/1‑251 (MQ=255)                                                                                                                                                                                 cGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACAcc                                                                     <  2:12156/251‑1 (MQ=255)                                                                                                                                                                                   cTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCtt                                                                   >  1:266191/1‑251 (MQ=255)                                                                                                                                                                                       cgcAACGGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCGATCCGGCAAACTCGGTGCCAGCTGTGGTGGATCGGATCAACATCACCTGCcgt                                                               <  2:556902/251‑1 (MQ=255)                                                                                                                                                                                         cAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCcgtcg                                                             <  1:81462/251‑1 (MQ=255)                                                                                                                                                                                                           aTGTGGCGTCTGCTGCTCGGTGAGTCGAAAAAAGGCTCCATCAACGGCCTCGGCTCACTGACTGGCGGTCAGGCGCTGC‑AACAGGAGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATg                                           <  1:436642/251‑1 (MQ=255)                                                                                                                                                                                                                  gTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcgcg                                    <  1:534487/251‑1 (MQ=255)                                                                                                                                                                                                                      gctgcACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGAGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGGCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAGGGTCCGCGGGc                                 <  1:427984/250‑1 (MQ=255)                                                                                                                                                                                                                        tgcACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAAACTGGCGGCCAGCATGTATCCGGGTCAGTCGCTCTATTCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCAtt                              >  1:8329/1‑251 (MQ=255)                                                                                                                                                                                                                          cACGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGa                            <  1:217372/251‑1 (MQ=255)                                                                                                                                                                                                                           aCGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCAttt                             >  1:86264/1‑248 (MQ=255)                                                                                                                                                                                                                            cGGTGAGTCGAAAAAAGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGACTATCTGTCGGGATGGCAGGAAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGACGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                          >  1:187219/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaaaaaGGCTACATCAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATcc                <  1:69497/251‑1 (MQ=255)                                                                                                                                                                                                                                            ggCTACATCAACGGCCTCGGCGCACTGACGGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTa          <  2:30241/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cAACGGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGGATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTCACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATACGGCAAACACGGTGCCCCCTGTGGGGGAGCGGGTCCACAACACCCTCCGTCGTGCCGAGCCGCTCCCATGGTCCCCGGGCATTGCGCCGGGCGATACGCGCCATGTccttt  >  2:501403/1‑247 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        ACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGC‑AACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATT  >  NC_000913/4217063‑4217556

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