Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,258	A→G	pseudogene (105/442 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,258	0	A	G	95.7%	135.9 / NA	46	pseudogene (105/442 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 A (1/1);  new base G (29/15);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCA  >  NC_000913/2173121‑2173404                                                                                                                                           |                                                                                                                                                   cAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGAtggtg                                                                                                                                         >  1:439559/1‑149 (MQ=255)         gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATc                                                                                                                                 >  1:498993/1‑149 (MQ=255)               aGTTGGGTTGTCAAACCAATAGCTTGCGTCGGTATCCACAGGGTGATGCTCATAATGATGACATCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCaaa                                                                                                                           <  2:716318/149‑1 (MQ=255)                    ggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGgca                                                                                                                      >  1:204279/1‑149 (MQ=255)                    ggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCCGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTAGCAAGATCGCCCAACGcca                                                                                                                      >  2:131103/1‑149 (MQ=255)                     gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGgcca                                                                                                                     >  2:746298/1‑147 (MQ=255)                      tgtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGacagc                                                                                                                    >  1:640057/1‑149 (MQ=255)                          cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGGTTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCt                                                                                                                >  1:389897/1‑149 (MQ=255)                           aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTg                                                                                                               >  2:853999/1‑149 (MQ=255)                            aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGa                                                                                                              <  2:748215/149‑1 (MQ=255)                                   aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCg                                                                                                       <  1:853999/149‑1 (MQ=255)                                    gTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCgg                                                                                                      <  2:72719/149‑1 (MQ=255)                                                      aGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACGGCAATTaa                                                                                    >  1:98631/1‑149 (MQ=255)                                                         gTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAt                                                                                 <  2:439559/149‑1 (MQ=255)                                                         gTGATGATCAGAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAt                                                                                 <  2:270563/149‑1 (MQ=255)                                                            aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGt                                                                              <  2:204279/149‑1 (MQ=255)                                                             tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGAAATTAAAATGGTg                                                                             >  1:522670/1‑149 (MQ=255)                                                                 cATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAAATAAAATGGTGAGTg                                                                         >  1:5213/1‑149 (MQ=255)                                                                  aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTgg                                                                        >  1:846301/1‑149 (MQ=255)                                                                   tAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATAAGCGGACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTggg                                                                       >  1:672925/1‑149 (MQ=255)                                                                             acCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCa                                                             <  2:266570/149‑1 (MQ=255)                                                                             acCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTCAGTGGGATAAAAATCa                                                             >  2:345983/1‑149 (MQ=255)                                                                                cGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGc                                                          >  1:527129/1‑149 (MQ=255)                                                                                    aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGc                                                      >  1:335393/1‑149 (MQ=255)                                                                                    aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGc                                                      <  2:646419/149‑1 (MQ=255)                                                                                                    gATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                      >  2:816381/1‑143 (MQ=255)                                                                                                    gATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                      <  2:335393/149‑7 (MQ=255)                                                                                                     aTTTCCACGATGCACGGCGACCGCCATCGCCCAGCCAGAA‑CGGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                     >  1:906673/1‑142 (MQ=255)                                                                                                          cACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTgg                                                                        >  2:270953/1‑109 (MQ=255)                                                                                                          cACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTgg                                                                        <  1:270953/109‑1 (MQ=255)                                                                                                          cACGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                              >  1:780985/1‑149 (MQ=255)                                                                                                           aCGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACGCAAAGAGAAATTACAACGGGTATTGTGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                             >  1:860849/1‑149 (MQ=255)                                                                                                            ggATGCAGGGCGAGCGCCATGG‑CCACCAAGAAGCCTATGGTGGCAAAATCGCCAAAGGGCAGTACCTAATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                            <  2:815868/148‑1 (MQ=255)                                                                                                            cGATGCACGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                            <  1:19489/149‑1 (MQ=255)                                                                                                                  aCGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa                      >  2:664936/1‑149 (MQ=255)                                                                                                                  aCGGCGACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa                      >  2:665169/1‑149 (MQ=255)                                                                                                                       gACCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                 <  2:846301/149‑1 (MQ=255)                                                                                                                        aCCGCCATCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                <  2:780985/149‑1 (MQ=255)                                                                                                                              aTCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg          >  1:816445/1‑149 (MQ=255)                                                                                                                               tCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATg                                                                                >  1:299881/1‑80 (MQ=255)                                                                                                                               tCG‑CCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATg                                                                                <  2:299881/80‑1 (MQ=255)                                                                                                                                   ccACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACACCATTGTGCTGCACCAACATCGCAGTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa      >  1:450538/1‑149 (MQ=255)                                                                                                                                    cACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa     >  2:828625/1‑149 (MQ=255)                                                                                                                                     aCGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCTGATCAAGGCCAAt    >  2:39705/1‑149 (MQ=255)                                                                                                                                     aCGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt    >  1:159637/1‑149 (MQ=255)                                                                                                                                      cGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc   >  2:326827/1‑149 (MQ=255)                                                                                                                                       gaaggaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa  <  2:350955/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                   CAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCG‑CCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCA  >  NC_000913/2173121‑2173404

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