Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,037 T→G 9.8% P476P (CCA→CCC)  pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,0370TG 9.8% 144.1 / 6.8 51P476P (CCA→CCC) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base G (3/2);  total (27/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 = 9.34e-01

CTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTA  >  NC_002947/460891‑461179
                                                                                                                                                   |                                                                                                                                               
cTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTgg                                                                                                                                               <  4:245107/149‑1 (MQ=255)
cTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTgg                                                                                                                                               <  3:19822/149‑1 (MQ=255)
  aCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGTGCCCTGATGGGGATCATGGACGGGGGCAGCGCCGTGCATGGCAGCCAGCAACTCGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTgg                                                                                                                                               >  1:205176/1‑147 (MQ=255)
   ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGc                                                                                                                                              <  6:177861/147‑1 (MQ=255)
    cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAgg                                                                                                                                           >  5:243176/1‑149 (MQ=255)
            gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGgc‑g                                                                                                                                    >  5:217967/1‑147 (MQ=255)
             cgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAg                                                                                                                                 <  6:217967/149‑1 (MQ=255)
               tCTGGCAGTTGCTTGAGCAGGGGCCTGATGCGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCTTCAAGGCCAAGGGTTGTTCCACTGCCCGCAGC‑TCTGCTGCTGGGAGGCCGGC‑GGAGTg                                                                                                                               <  6:167983/149‑1 (MQ=255)
                   gCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCaa                                                                                                                            <  3:142286/148‑1 (MQ=255)
                   gCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAAc                                                                                                                           >  3:218235/1‑149 (MQ=255)
                      gTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGc                                                                                                                        <  7:139956/149‑1 (MQ=255)
                             gAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTAt                                                                                                                  <  4:218235/148‑1 (MQ=255)
                              aGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTc                                                                                                                >  2:70016/1‑149 (MQ=255)
                                cccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTccc                                                                                                              >  5:233990/4‑149 (MQ=255)
                                cAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTccc                                                                                                              >  3:2175/1‑149 (MQ=255)
                                  cggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCa                                                                                                                                  <  5:8079/126‑1 (MQ=255)
                                  cggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCa                                                                                                                                  >  6:8079/2‑127 (MQ=255)
                                        tGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                      >  5:42155/1‑149 (MQ=255)
                                            ggggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCActgctg                                                                                                   >  4:178344/1‑148 (MQ=255)
                                             gggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCAACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCActgctg                                                                                                   <  1:70016/148‑1 (MQ=255)
                                              ggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTAc                                                                                                >  5:275808/1‑149 (MQ=255)
                                                aTCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTg                                                                                              <  5:197404/149‑1 (MQ=255)
                                                     ggCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGtt                                                                                         >  2:270621/1‑149 (MQ=255)
                                                      gCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGtt                                                                                         <  7:169649/148‑1 (MQ=255)
                                                          gtggCATCGCGGTGCAGGGCAGCCAGCAACTGTCAGGTGGATTTAGCATCAAGGCCTAGTCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTTCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACg                                                                                    <  2:205176/147‑1 (MQ=255)
                                                            ggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGa                                                                                  <  7:164854/149‑1 (MQ=255)
                                                             gCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAg                                                                                 <  3:178344/149‑1 (MQ=255)
                                                             gCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAg                                                                                 <  6:42155/149‑1 (MQ=255)
                                                                   ccGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCgg                                                                           <  4:2175/149‑1 (MQ=255)
                                                                          gggCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATACCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAAc                                                                    >  8:146226/1‑149 (MQ=255)
                                                                           ggCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACg                                                                   >  4:191340/1‑149 (MQ=255)
                                                                              agccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTc                                                                <  6:233990/149‑1 (MQ=255)
                                                                                ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                              <  3:218953/149‑1 (MQ=255)
                                                                                              gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAgg                                                >  2:196939/1‑149 (MQ=255)
                                                                                                   ttGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTg                                           >  7:218951/1‑149 (MQ=255)
                                                                                                      aGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGtt                                         <  1:270621/148‑1 (MQ=255)
                                                                                                            aaggccaaggcTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGca                                   >  2:18256/1‑148 (MQ=255)
                                                                                                               gccaaggcTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCgg                                                                           >  2:261953/1‑105 (MQ=255)
                                                                                                               gccaaggcTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCgg                                                                           <  1:261953/105‑1 (MQ=255)
                                                                                                                          ttCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCt                    >  2:70511/1‑149 (MQ=255)
                                                                                                                             cACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGCGGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTgg                  >  1:36440/1‑149 (MQ=255)
                                                                                                                             cACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGgc                 >  3:87058/1‑149 (MQ=255)
                                                                                                                               cTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGgcg                >  4:38760/1‑148 (MQ=255)
                                                                                                                                tGCCAGCAGC‑AGTGCTGCGGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTaa                         >  2:228758/1‑138 (MQ=255)
                                                                                                                                tGCCAGCAGC‑AGTGCTGCGGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTaa                         <  1:228758/138‑1 (MQ=255)
                                                                                                                                 gCCAGCAGC‑AGTGCTGCGGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCt                             <  1:188643/133‑1 (MQ=255)
                                                                                                                                 gCCAGCAGC‑AGTGCTGCGGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCt                             >  2:188643/1‑133 (MQ=255)
                                                                                                                                  cccGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCAt             <  8:243142/148‑1 (MQ=255)
                                                                                                                                  ccagcagc‑aGTGCTGCGGGCAGTCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCa              >  7:7712/4‑147 (MQ=255)
                                                                                                                                   ccGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAg           >  7:217308/1‑149 (MQ=255)
                                                                                                                                      cagc‑aCTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGt         >  1:173970/1‑148 (MQ=255)
                                                                                                                                               ctgctgGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTa  <  7:58175/147‑1 (MQ=255)
                                                                                                                                               ctgctgGCAGGCCGGC‑GCAGTGCAACTGCCAGTATACCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTa  <  7:146226/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGC‑ACTGCTGCTGGCAGGCCGGC‑GCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTA  >  NC_002947/460891‑461179

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: ‑