Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,076 C→A pseudogene (287/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_0009132,173,0760CA96.5% 187.6 / ‑3.2 57pseudogene (287/442 nt)gatCpseudogene, 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 C (0/2);  new base A (26/29);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCC  >  NC_000913/2172980‑2173171
                                                                                                |                                                                                               
tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                              >  1:976304/1‑100 (MQ=255)
tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                              <  1:953860/100‑1 (MQ=255)
 gAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGc                                                                                               >  2:1425722/1‑98 (MQ=255)
 gAAAATACCGGTCAGATAAATTGCACCAATAATAATGAAACCGGAAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGTGAGAACCGccc                                                                                             <  2:271102/100‑1 (MQ=255)
  aaagTACCGGTCAGATAAATTTCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCt                                                                                            <  1:1196438/100‑1 (MQ=255)
  aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAGACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                            >  2:845478/1‑100 (MQ=255)
  aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                            <  2:901039/100‑1 (MQ=255)
  aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                            <  1:1313846/100‑1 (MQ=255)
  aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                            >  2:1046722/1‑100 (MQ=255)
   aaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTg                                                                                           <  2:937286/100‑1 (MQ=255)
        ncGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCa                                                                                      >  1:934927/2‑100 (MQ=255)
          ggTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCAtt                                                                                    >  2:824649/1‑100 (MQ=255)
           gTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTg                                                                                   <  2:579919/100‑1 (MQ=255)
           gTCAGAGAAATTGCACCGATAATAATGAAACCGTGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACATCCCTGATCCATTg                                                                                   <  2:243408/100‑1 (MQ=255)
            tCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGa                                                                                  <  1:703066/100‑1 (MQ=255)
             cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGaa                                                                                 >  1:896770/1‑100 (MQ=255)
              aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                <  1:447831/100‑1 (MQ=255)
              aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                <  2:1328462/100‑1 (MQ=255)
              aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                <  2:1453546/100‑1 (MQ=255)
               gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGc                                                                               >  2:780723/1‑100 (MQ=255)
                aTAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGcc                                                                              <  1:1017504/100‑1 (MQ=255)
                 tAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGcca                                                                             <  2:117926/100‑1 (MQ=255)
                   aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacc                                                                           <  2:770525/100‑1 (MQ=255)
                    aTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacca                                                                          >  2:1121898/1‑100 (MQ=255)
                      tGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATa                                                                        <  2:366617/100‑1 (MQ=255)
                       gCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAc                                                                       >  2:1502180/1‑100 (MQ=255)
                        cACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAcc                                                                      >  1:1505359/1‑100 (MQ=255)
                           cGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTAATCCATTGAAGCCACCATAccccc                                                                   >  1:996633/1‑100 (MQ=255)
                            gATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCg                                                                  <  2:1319876/100‑1 (MQ=255)
                             ataataATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCgg                                                                 >  2:241026/1‑100 (MQ=255)
                             ataataATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCg                                                                  >  2:242004/1‑99 (MQ=255)
                                ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtt                                                              >  2:663918/1‑100 (MQ=255)
                                  aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                            <  1:900990/100‑1 (MQ=255)
                                   atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTa                                                           <  2:1505350/100‑1 (MQ=255)
                                   atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTa                                                           <  2:482178/100‑1 (MQ=255)
                                      aaaCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACg                                                        >  2:1248454/1‑100 (MQ=255)
                                              aaTATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCa                                                <  1:1567744/100‑1 (MQ=255)
                                                    ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgcc                                          <  2:390869/100‑1 (MQ=255)
                                                     tGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccg                                         <  1:394631/100‑1 (MQ=255)
                                                        ggggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCa                                      >  1:1040631/1‑100 (MQ=255)
                                                         gggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAg                                     <  2:133679/100‑1 (MQ=255)
                                                              aaaaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTggg                                <  1:82657/100‑1 (MQ=255)
                                                                 aCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAACCGCCAGTTGGgtgt                             >  1:1554672/1‑100 (MQ=255)
                                                                  ccTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCCGTTGGGtgtg                            >  2:249295/1‑100 (MQ=255)
                                                                  ccTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGACACCATAGCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgtg                              >  2:321809/1‑98 (MQ=255)
                                                                    tGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCa                          >  1:704579/1‑100 (MQ=255)
                                                                     gAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCaa                         >  1:381336/1‑100 (MQ=255)
                                                                         cAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACACCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGCGTGCAAACCa                     >  2:265359/1‑100 (MQ=255)
                                                                          aGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCACACCa                     >  2:660524/1‑99 (MQ=255)
                                                                                cAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGttt              <  1:824654/100‑1 (MQ=255)
                                                                                  ggTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTgn            <  1:251908/100‑2 (MQ=255)
                                                                                  ggTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc            <  2:304495/100‑1 (MQ=255)
                                                                                    tAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACTCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGt          <  2:1051943/100‑1 (MQ=255)
                                                                                     aaTGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTc         >  2:287500/1‑100 (MQ=255)
                                                                                       tGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTATAACGCCCCAGCATTAGCCGCCAGTCGGGTGTGCAAACCAATAGTTTGCGTcgc       >  2:388626/1‑100 (MQ=255)
                                                                                          gagaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGAt    <  1:639661/100‑1 (MQ=255)
                                                                                            gaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGAttc  <  1:1068033/100‑3 (MQ=255)
                                                                                                |                                                                                               
TGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCC  >  NC_000913/2172980‑2173171

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