Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	313,726	A→G	pseudogene (80/90 nt)	ykgP ←	pseudogene, oxidoreductase family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	313,726	0	A	G	100.0%	154.7 / NA	41	pseudogene (80/90 nt)	ykgP	pseudogene, oxidoreductase family
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/24);  total (17/24)

ACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGG  >  NC_000913/313591‑313870                                                                                                                                        |                                                                                                                                                 acacTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAAccagccag                                                                                                                                   <  1:24996/151‑1 (MQ=255)         ggCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTc                                                                                                                           <  1:632627/151‑1 (MQ=255)         ggCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTc                                                                                                                           <  2:324133/151‑1 (MQ=255)                   aTGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGa                                                                                                                  <  2:63839/150‑1 (MQ=255)                    tGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAAt                                                                                                                >  2:423983/1‑151 (MQ=255)                              tGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCAt                                                                                                      >  1:140034/1‑151 (MQ=255)                                   cAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGttt                                                                                                 <  1:272395/151‑1 (MQ=255)                                         cGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCgg                                                                                           <  1:77203/151‑1 (MQ=255)                                         cGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCgg                                                                                           >  2:451146/1‑151 (MQ=255)                                           ggTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGaa                                                                                         <  1:68406/151‑1 (MQ=255)                                               aaTTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATgg                                                                                      <  2:575305/150‑1 (MQ=255)                                                   aTCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAAt                                                                                 >  1:502839/1‑151 (MQ=255)                                                      tGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTc                                                                               >  2:365799/1‑150 (MQ=255)                                                      tGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTc                                                                               >  2:190723/1‑150 (MQ=255)                                                       gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAAtt                                                                                <  1:594034/148‑1 (MQ=255)                                                       gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCAt                                                                             <  1:339492/151‑1 (MQ=255)                                                                   gAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTtgtg                                                                 <  2:489896/151‑1 (MQ=255)                                                                   gAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTtgtg                                                                 <  1:524369/151‑1 (MQ=255)                                                                   gAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTtgt                                                                  <  1:603634/150‑1 (MQ=255)                                                                               gTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATg                                                     >  2:191632/1‑151 (MQ=255)                                                                                  aCTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATgaga                                                  <  1:229845/151‑1 (MQ=255)                                                                                        tgatttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCAt                                             <  1:90662/150‑1 (MQ=255)                                                                                         gatttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATaa                                           <  1:228813/151‑1 (MQ=255)                                                                                               atAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTaa                                      >  1:67641/1‑150 (MQ=255)                                                                                                  gTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAggg                                  <  1:222745/151‑1 (MQ=255)                                                                                                   ttttACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAggg                                  >  1:436953/1‑150 (MQ=255)                                                                                                             gggAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGAt                       <  1:142861/151‑1 (MQ=255)                                                                                                                   ataataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACg                 <  2:593107/151‑1 (MQ=255)                                                                                                                    taataGATATTTAATCCATTGACGGAAAACAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACgg                >  2:418974/1‑151 (MQ=255)                                                                                                                     aataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACgg                <  1:365806/150‑1 (MQ=255)                                                                                                                     aataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGc               >  2:295467/1‑151 (MQ=255)                                                                                                                     aataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGc               >  2:592642/1‑151 (MQ=255)                                                                                                                      ataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGc               >  1:553117/1‑150 (MQ=255)                                                                                                                           tatTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAtt         >  2:151756/1‑151 (MQ=255)                                                                                                                            atTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAtt         <  1:340393/150‑1 (MQ=255)                                                                                                                            atTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAtt         <  2:233883/150‑1 (MQ=255)                                                                                                                              ttAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGt      >  1:236847/1‑151 (MQ=255)                                                                                                                               tAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTc     >  1:612395/1‑151 (MQ=255)                                                                                                                                aaTCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCa    <  1:407112/151‑1 (MQ=255)                                                                                                                                aaTCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCa    <  2:553106/151‑1 (MQ=255)                                                                                                                                  tCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAgg  >  2:237227/1‑151 (MQ=255)                                                                                                                                        |                                                                                                                                                 ACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGG  >  NC_000913/313591‑313870

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