Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	pseudogene (915‑916/1358 nt)	gatC ←	galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,361	0	C	.	100.0%	70.5 / NA	15	pseudogene (918/1358 nt)	gatC	galactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/5);  total (10/5)
*	NC_000913	2,173,362	0	C	.	100.0%	70.5 / NA	15	pseudogene (917/1358 nt)	gatC	galactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/5);  total (10/5)

ACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTG  >  NC_000913/2173299‑2173414                                                               ||                                                     aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:1241305/1‑64 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:372130/1‑64 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:1357856/1‑64 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:2805321/1‑64 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:2839219/1‑64 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:2249199/1‑64 (MQ=255)   ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   >  1:1005695/1‑62 (MQ=255)   ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   <  1:354993/65‑4 (MQ=255)   ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   >  1:2167506/1‑62 (MQ=255)   ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   <  1:2309256/65‑4 (MQ=255)          aaaaaCAGCAATTAAAATGGTGAGTGGGATAAAAACCGGGCTTGCCGATACCA‑‑CCGCCGTATGTc                                          <  1:3230631/63‑1 (MQ=255)             aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                       >  1:817748/1‑65 (MQ=255)                  cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcaga                                  >  1:3162501/1‑64 (MQ=255)                     tgaaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                               <  1:2863772/63‑1 (MQ=255)                        aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                            >  1:134021/1‑65 (MQ=255)                           tGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                         <  1:2990773/65‑1 (MQ=255)                            ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                        <  1:2979036/65‑1 (MQ=255)                                   gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                 >  1:100249/1‑65 (MQ=255)                                   gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                 >  1:675981/1‑65 (MQ=255)                                   gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                 >  1:598848/1‑65 (MQ=255)                                   gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                 >  1:2537100/1‑65 (MQ=255)                                         aaaaTCGGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg           <  1:716403/65‑1 (MQ=255)                                           aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa         >  1:2184132/1‑65 (MQ=255)                                           aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGa         >  1:2310915/1‑65 (MQ=255)                                              cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt      <  1:2014804/65‑1 (MQ=255)                                              cAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt      <  1:445731/65‑1 (MQ=255)                                                ggCTTGGCCATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc    <  1:18214/65‑1 (MQ=255)                                                 gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt   >  1:1688288/1‑65 (MQ=255)                                                  cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg  >  1:1547716/1‑65 (MQ=255)                                                               ||                                                     ACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTG  >  NC_000913/2173299‑2173414

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