Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,503,789 C→T 13.0% intergenic (+38/‑60) xanP → / → tusD xanthine permease/sulfur transfer protein complex subunit TusD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,503,7890CT13.0% 103.1 / 4.5 45intergenic (+38/‑60)xanP/tusDxanthine permease/sulfur transfer protein complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base C (15/21);  new base T (2/3);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGG  >  NC_002947/4503641‑4503937
                                                                                                                                                    |                                                                                                                                                    
gcgcACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGc                                                                                                                                                      >  2:284601/1‑149 (MQ=255)
       tGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGcac                                                                                                                                                    <  8:145921/144‑1 (MQ=255)
       tGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGcac                                                                                                                                                    >  7:145921/1‑144 (MQ=255)
          cggcggCGATTCGCAATGTGCTTGAGGCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggc                                                                                                                                            <  8:364931/149‑1 (MQ=255)
            gcggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgc                                                                                                                                          <  5:341099/149‑1 (MQ=255)
            gcggagATTAGCAATGTGCTGGTGTCGGGGGTGGCTACCGGTGGGATTTGTGCGCTGTTGCTGAACTGGTTCTTGCCGTAGGGCAAGGAGCAGGCTTGAACCTTGAGACTTAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgc                                                                                                                                          <  8:311830/149‑1 (MQ=255)
              ggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTTCCGGGGGGTTTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg                                                                                                                                        <  2:203465/149‑1 (MQ=255)
                       cAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCCGCCCCGCTGGGCCCAGCGGGCGGGGGGgtttttt                                                                                                                               >  6:42746/1‑148 (MQ=255)
                         aTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGGGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCt                                                                                                                             <  5:36821/149‑1 (MQ=255)
                               tGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGGCCAAGCGGGCGCGGGCTTTTTTGCTTTAACa                                                                                                                       >  4:236606/1‑149 (MQ=255)
                                ggAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAt                                                                                                                      <  8:365932/149‑1 (MQ=255)
                                      gggTGTGGCTACCTGGGGGATTTGTGCGCGGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTAGAACCTTGAGCCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCAc                                                                                                                 <  4:37224/148‑1 (MQ=255)
                                        gtgtGGCTACCGGGGGGATTTGTGCGCTGGTGTTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATg                                                                                                                     >  1:243003/1‑142 (MQ=255)
                                        gtgtGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATg                                                                                                                     <  2:243003/142‑1 (MQ=255)
                                              cTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCttt                                                                                                        <  4:163264/149‑1 (MQ=255)
                                                 ccGGGGGGATTTGTGCGCTGGTGCTGAACTTGTTCTTGCCGGATGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGACTAGCCCGCGCCCGCTTTCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                     <  5:363527/149‑1 (MQ=255)
                                                 ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                     <  4:147693/149‑1 (MQ=255)
                                                    gggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCttt                                                                                                                           <  3:151426/124‑1 (MQ=255)
                                                    gggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCttt                                                                                                                           >  4:151426/1‑124 (MQ=255)
                                                            tgtgCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGcc                                                                                           <  1:284601/148‑1 (MQ=255)
                                                                     gTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTc                                                                                  <  7:91082/148‑1 (MQ=255)
                                                                          gAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATc                                                                            >  4:120322/1‑149 (MQ=255)
                                                                          gAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATc                                                                            >  7:203312/1‑149 (MQ=255)
                                                                               ggTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCAt                                                                                                    <  5:226673/120‑1 (MQ=255)
                                                                               ggTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt                                                                       >  6:310415/1‑149 (MQ=255)
                                                                                   cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggc                                                                                                                                            >  8:124147/1‑76 (MQ=255)
                                                                                   cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggc                                                                                                                                            <  7:124147/76‑1 (MQ=255)
                                                                                        cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCCCGCCCCCTTGGACAAGCCGGCGCCGGGTTTTTTTGGTTATAATGGCCCCAGTCCATTGCATGAGTTAGCCAATAAATTAGGAAAGGCGGTGTTATgccc                                                              >  2:220315/1‑149 (MQ=255)
                                                                                               cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGc                                                                                                                                                      >  5:44379/1‑54 (MQ=255)
                                                                                               cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGc                                                                                                                                                      <  6:44379/54‑1 (MQ=255)
                                                                                               cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGGACAAGCGGGCGCGGGGGTTTTTGCTTTATCATGGCACCATTGCTTTGGATGAGTTAGGCATGAAATTCGCGATCGCGGTGTTCTCCCCGGCCCAt                                                       >  2:94292/1‑149 (MQ=255)
                                                                                                aaGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCTGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATg                                                      >  2:213417/1‑149 (MQ=255)
                                                                                                    aGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGcc                                                  <  3:120322/149‑1 (MQ=255)
                                                                                                         gCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCtt                                             <  3:236606/149‑1 (MQ=255)
                                                                                                           ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTcg                                           >  5:169617/1‑149 (MQ=255)
                                                                                                             gAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTcg                                           <  5:310415/147‑1 (MQ=255)
                                                                                                               aCCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCGGGCGCGGGGTTTTTTGCTTTATCATGGCACCATTGCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGGGTTGTCCCCGGCCCATGCGCCCTCTTCGCGAc                                       >  2:127110/1‑149 (MQ=255)
                                                                                                                         tGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGcc                                                  <  3:128053/128‑1 (MQ=255)
                                                                                                                         tGAATGAGAGTAGCCCGCGCCCCCTTGCACAAACGGGCGCCGGCTTTTTTGCTTTATCATGGCACCATTCCATTGCATGAGTTAGCCATGAAATTAGCTATCGAGGTGTTCTGCCCGGCCCAAGCGcc                                                  >  4:128053/1‑128 (MQ=255)
                                                                                                                              gagaGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCg                        >  5:392769/1‑149 (MQ=255)
                                                                                                                                gagTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGcc                      <  7:312172/149‑1 (MQ=255)
                                                                                                                                   tAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAg                   >  4:325032/1‑149 (MQ=255)
                                                                                                                                       ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCgg               <  6:313492/149‑1 (MQ=255)
                                                                                                                                        cgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACACGCCCTGCGCTTCGCCGAGGCGGt              >  3:274712/1‑149 (MQ=255)
                                                                                                                                         gcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg             <  6:392769/149‑1 (MQ=255)
                                                                                                                                            cccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCATTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGGTTCGCCGAGGCGGGGCTg          >  8:191801/1‑149 (MQ=255)
                                                                                                                                             ccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTgg         >  3:74432/1‑149 (MQ=255)
                                                                                                                                              cgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCTCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGc        >  8:171620/1‑149 (MQ=255)
                                                                                                                                               gcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCATTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCCAGGCGGTGCTGGcc       >  8:177697/1‑149 (MQ=255)
                                                                                                                                                   gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcgg     <  8:203312/147‑1 (MQ=255)
                                                                                                                                                    cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg  <  6:169617/149‑1 (MQ=255)
                                                                                                                                                    cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcg   <  3:325032/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGG  >  NC_002947/4503641‑4503937

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