Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,933	T→C	pseudogene (346/914 nt)	gatC ←	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,933	0	T	C	100.0%	128.2 / NA	36	pseudogene (346/914 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)

CAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGTCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTCGAAATTTTCCGCCATCG  >  NC_000913/2173791‑2174079                                                                                                                                               |                                                                                                                                                   cAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACggg                                                                                                                                            >  2:251439/1‑151 (MQ=255)   gTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTc                                                                                                                                          <  2:542954/151‑1 (MQ=255)      gCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGt                                                                                                                                       >  2:506979/1‑151 (MQ=255)       ccGAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTa                                                                                                                                      <  2:2768/151‑1 (MQ=255)         gAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAAc                                                                                                                                    >  1:150111/1‑151 (MQ=255)          aGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAAc                                                                                                                                    >  2:591041/1‑150 (MQ=255)              tataAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGc                                                                                                                               >  2:304513/1‑151 (MQ=255)                  aacaaacGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTCACCACCCGCGCCATACGGGTCAGTAACATCGCCAc                                                                                                                            >  1:218228/1‑150 (MQ=255)                       acGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTaa                                                                                                                       >  1:150995/1‑150 (MQ=255)                       acGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAAc                                                                                                                      >  1:67034/1‑151 (MQ=255)                                     aCCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCg                                                                                                        >  1:220382/1‑151 (MQ=255)                                          acCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATc                                                                                                   <  2:641949/151‑1 (MQ=255)                                          acCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATc                                                                                                   >  1:637182/1‑151 (MQ=255)                                           ccTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATc                                                                                                   <  2:451716/150‑1 (MQ=255)                                            cTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGc                                                                                                 >  1:665299/1‑151 (MQ=255)                                                cATCCCTATCATCAATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGccac                                                                                              >  2:229950/1‑150 (MQ=255)                                                         catcCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAAtt                                                                                    <  1:334315/151‑1 (MQ=255)                                                         catcCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAAt                                                                                     <  1:678568/150‑1 (MQ=255)                                                                 aaCCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGc                                                                            >  1:557403/1‑151 (MQ=255)                                                                    cGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCa                                                                         <  1:673375/151‑1 (MQ=255)                                                                         gCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCa                                                                    <  2:86870/151‑1 (MQ=255)                                                                         gCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCa                                                                    <  1:309493/151‑1 (MQ=255)                                                                              aTGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGt                                                               <  2:15066/151‑1 (MQ=255)                                                                                           ccGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGgcc                                                  >  1:428887/1‑151 (MQ=255)                                                                                             ggTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGgcc                                                  >  1:573553/1‑149 (MQ=255)                                                                                                      cATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATTCGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCACGGCCAGCCAACATc                                        <  2:457302/150‑1 (MQ=255)                                                                                                      cATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCg                                       <  2:629794/151‑1 (MQ=255)                                                                                                         atGCCAGATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGAcc                                    >  1:295448/1‑151 (MQ=255)                                                                                                               gATATTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGc                              <  2:563890/151‑1 (MQ=255)                                                                                                                 tatTCCAGATATCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGc                              >  2:451221/1‑149 (MQ=255)                                                                                                                           atCAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTc                   <  2:417684/150‑1 (MQ=255)                                                                                                                             cAACATTTACCGCCCGCGCCATACGGTTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCTCCAGCGCAATTTCCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGt                    <  2:212874/147‑1 (MQ=255)                                                                                                                             cAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTCGa                 >  2:579099/1‑150 (MQ=255)                                                                                                                             cAACATTTACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTCGa                 >  2:317186/1‑150 (MQ=255)                                                                                                                                  tttACCACCCGCGCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTCGAAAtttt           >  2:36861/1‑151 (MQ=255)                                                                                                                                           cgcgCCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTCGAAATTTTCCGCCATCg  >  1:447343/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                   CAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACCTGCCATCCCTATCATCCATGAACCGGTTGCCAGATGCAGCAACGCGCCGGTGAAGGTCATATGCCAGATATTCCAGATATCAACATTTACCACCCGCGTCATACGGGTCAGTAACATCGCCACGTTAACCAGAATCGCAATCGGAATCGCCACCAGCGCAATTTGCGACGCCCAGGTCATTGGTGAAGAGCCCGGCCAGCCAACATCGACCACATGCAGATTCAGGTCGAAATTTTCCGCCATCG  >  NC_000913/2173791‑2174079

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