Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 165,269 C→T 10.7% G63D (GGT→GAT)  pntB ← pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947165,2690CT10.7% 150.3 / 7.7 56G63D (GGT→GAT) pntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base T (3/3);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

CGGCGATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGCACAC  >  NC_002947/165123‑165414
                                                                                                                                                  |                                                                                                                                                 
cggcgATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCg                                                                                                                                                 >  2:111610/1‑149 (MQ=255)
    gATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGGTATCGGCta                                                                                                                                             <  5:187385/149‑2 (MQ=255)
        gCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGc                                                                                                                                         <  8:80471/149‑1 (MQ=255)
          aaTGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCaa                                                                                                                                       <  1:315834/149‑1 (MQ=255)
                cacCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAg                                                                                                                                 >  6:4871/1‑149 (MQ=255)
                          aGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCtt                                                                                                                       <  3:138169/149‑1 (MQ=255)
                              cGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGACCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTtata                                                                                                                    >  2:356842/1‑148 (MQ=255)
                                aTCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGACCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAAt                                                                                                                 <  1:356842/149‑1 (MQ=255)
                                   aTGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACTAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGa                                                                                                               <  1:111610/148‑1 (MQ=255)
                                     gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGGTATCGGCTATGTCa                                                                                                                                        <  6:50401/121‑1 (MQ=255)
                                     gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGGTATCGGCTATGTCa                                                                                                                                        >  5:50401/1‑121 (MQ=255)
                                     gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGa                                                                                                               <  7:347642/146‑1 (MQ=255)
                                     gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGa                                                                                                               >  8:347642/1‑146 (MQ=255)
                                       tgtgCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGcc                                                                                                           >  2:22002/1‑148 (MQ=255)
                                           cATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCg                                                                                                      >  2:369589/1‑149 (MQ=255)
                                            aTGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCg                                                                                                      >  7:250583/1‑148 (MQ=255)
                                             tGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCaa                                                                                                                           >  1:131423/1‑126 (MQ=255)
                                             tGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCaa                                                                                                                           <  2:131423/126‑1 (MQ=255)
                                              gAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAg                                                                                                              >  5:316380/1‑138 (MQ=255)
                                              gAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAg                                                                                                              <  6:316380/138‑1 (MQ=255)
                                               aaaGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAg                                                                                                   <  6:221486/148‑1 (MQ=255)
                                                 aGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGa                                                                                                <  4:272534/149‑1 (MQ=255)
                                                             gggCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTg                                                                                                 <  4:214114/136‑1 (MQ=255)
                                                             gggCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTg                                                                                                 >  3:214114/1‑136 (MQ=255)
                                                             gggCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGcc                                                                                    >  2:197027/1‑149 (MQ=255)
                                                                 aTTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCAt                                                                                >  7:215393/1‑149 (MQ=255)
                                                                    ttGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGc                                                                              <  1:369589/148‑1 (MQ=255)
                                                                        tCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGAt                                                                          >  7:85071/1‑148 (MQ=255)
                                                                           tttCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATg                                                                      <  1:197027/149‑1 (MQ=255)
                                                                              cGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGcc                                                                    <  5:4871/148‑1 (MQ=255)
                                                                                        ggCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTg                                                          <  8:250583/148‑1 (MQ=255)
                                                                                               aTCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTgccgcgccg                                                  >  4:111058/1‑149 (MQ=255)
                                                                                                   aaCCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAg                                              >  4:217167/1‑149 (MQ=255)
                                                                                                   aaCCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAg                                              >  5:288230/1‑149 (MQ=255)
                                                                                                   aaCCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAg                                              <  8:215393/149‑1 (MQ=255)
                                                                                                    gtcggcggCACCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGAt                                                            <  5:9182/132‑1 (MQ=255)
                                                                                                    gtcggcggCACCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGAt                                                            >  6:9182/3‑134 (MQ=255)
                                                                                                        gcggTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAggtggt                                         <  3:111058/149‑1 (MQ=255)
                                                                                                        gcggTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAggtggt                                         >  8:210624/1‑149 (MQ=255)
                                                                                                           gTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAAc                                                                <  2:32402/123‑1 (MQ=255)
                                                                                                           gTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAAc                                                                >  1:32402/1‑123 (MQ=255)
                                                                                                                gccgACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGc                                  >  8:153576/1‑148 (MQ=255)
                                                                                                                        gcagGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGACAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGccc                         >  5:218621/1‑149 (MQ=255)
                                                                                                                            gCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGa                     <  8:33905/149‑1 (MQ=255)
                                                                                                                                   ttgccacaGCCGGTATCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTg              >  6:264987/8‑149 (MQ=255)
                                                                                                                                    tgccacaGCCGGTATCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTg              >  1:44561/7‑148 (MQ=255)
                                                                                                                                      aCATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGa            <  3:217167/148‑1 (MQ=255)
                                                                                                                                       cacaGCCGGTATCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATg          <  5:264987/146‑1 (MQ=255)
                                                                                                                                        aTAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCt               <  5:139032/143‑1 (MQ=255)
                                                                                                                                        aTAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCt               >  6:139032/1‑143 (MQ=255)
                                                                                                                                          aGCCGATACCGGCTGTGGCAAGTTCAGCGCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAg       >  1:211778/1‑149 (MQ=255)
                                                                                                                                          aGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAg       >  7:298971/1‑149 (MQ=255)
                                                                                                                                           gCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAg       >  3:17144/1‑148 (MQ=255)
                                                                                                                                            ccGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGca     <  4:17144/149‑1 (MQ=255)
                                                                                                                                            ccGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGca     <  1:22002/149‑1 (MQ=255)
                                                                                                                                               aTACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGCacac  <  2:211778/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGGCGATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGCACAC  >  NC_002947/165123‑165414

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