Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,606,996 C→A 8.0% V707V (GTG→GTT)  PP_3182 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,606,9960CA 8.0% 154.1 / 4.5 50V707V (GTG→GTT) PP_3182GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (2/2);  total (25/25)
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.22e-01

TCAGTACACGCATGCCTTGTCCTCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATGCAGGTGGG  >  NC_002947/3606852‑3607144
                                                                                                                                                   |                                                                                                                                                     
tCAGTACACGCATGCCTTGTACTCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCAcc                                                                                                                                                    >  1:43430/1‑148 (MQ=255)
              ccTTGTCCTCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCcaccac                                                                                                                                      <  2:43430/148‑1 (MQ=255)
                   tCCTCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTccag                                                                                                                                >  2:282196/1‑149 (MQ=255)
                      tCGATACCGTTGAATGGCTGGTAAAGCCACAGCTGCGCGAGTGACTGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTcc                                                                                                                                                       <  8:145383/123‑1 (MQ=255)
                      tCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTccgggcc                                                                                                                             >  4:246146/1‑149 (MQ=255)
                      tCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGAGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTCGGCAATATCGCAGCCGGCAAGATA‑CAGGTACACCAC‑‑CTGTACACGGGTTTcc                                                                                                                                                       >  7:145383/1‑123 (MQ=255)
                        gATACCGTTTACGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCGGGCTTATGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCAccccc                                                                                                                                                 <  4:256262/127‑1 (MQ=255)
                        gATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTccaggcc                                                                                                                             >  5:283761/1‑147 (MQ=255)
                             cGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGgcgc                                                                                                                      >  3:310330/1‑149 (MQ=255)
                               ttGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCtt                                                                                                                    >  6:65400/1‑149 (MQ=255)
                               ttGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCtt                                                                                                                    >  6:70541/1‑149 (MQ=255)
                                    gggCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCa                                                                                                               <  5:70541/149‑1 (MQ=255)
                                    gggCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCa                                                                                                               <  5:65400/149‑1 (MQ=255)
                                          gTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGgatga                                                                                                         >  5:162264/1‑149 (MQ=255)
                                               cccACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGAtggtg                                                                                                    <  3:204827/149‑1 (MQ=255)
                                                ccACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTg                                                                                                                   >  5:172690/1‑133 (MQ=255)
                                                ccACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTg                                                                                                                   <  6:172690/133‑1 (MQ=255)
                                                  acaGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCGGGCCCAGGCGCTTGGCCAGGATGATGGTggag                                                                                                 <  3:246146/149‑1 (MQ=255)
                                                  acaGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTggag                                                                                                 >  5:64821/1‑149 (MQ=255)
                                                   caGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTggagg                                                                                                <  5:320117/149‑1 (MQ=255)
                                                      cTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGaca                                                                                             <  7:186718/149‑1 (MQ=255)
                                                       tcggcCACCTACTGTGGCGGGCTTAAGTGTTAGCCCTGGGAAATATCGCCGCC‑GCAAACTTCCCGGTAACCCACTACCCTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGac                                                                                              <  1:282196/146‑1 (MQ=255)
                                                            cAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTg                                                                                       <  6:283761/149‑1 (MQ=255)
                                                             aGCGACAGTGGCCGGATGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTgg                                                                                      >  4:335507/1‑149 (MQ=255)
                                                              gCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTggg                                                                                     <  4:310330/149‑1 (MQ=255)
                                                               cGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACTCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGc                                                                                    >  5:205531/1‑149 (MQ=255)
                                                                  cAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCCCGGGTTTCCACCCCCTCGGCCCCCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTg‑t                                                                                >  5:283296/1‑149 (MQ=255)
                                                                    gTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACTCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTg‑tcg                                                                              <  6:205531/149‑1 (MQ=255)
                                                                        ccGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCAt                                                                          >  2:274535/1‑149 (MQ=255)
                                                                         cGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATg                                                                         >  8:89977/1‑149 (MQ=255)
                                                                                    gTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGAc                                                               <  6:162264/148‑1 (MQ=255)
                                                                                                    tCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCg                                               >  2:9948/1‑148 (MQ=255)
                                                                                                          ccGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAAcgatc                                        >  2:61456/1‑149 (MQ=255)
                                                                                                             gCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAAcgatcgat                                     <  1:61456/149‑1 (MQ=255)
                                                                                                              cAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACgatcgatc                                    <  6:64821/149‑1 (MQ=255)
                                                                                                                  cTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCa                                <  1:9948/149‑1 (MQ=255)
                                                                                                                      cAGGTACACCAC‑‑CTGTTCACCGGTTTCCACCCCCTCCGCCACGACTTCCAGGCCCAGGCTCTTGGCCAGGATTATGGTGTAGGACTCGATTGGGCTG‑TCGTGATGGCTGTGGGTCAGCGGCGCGATCAGCGATCGCTCGATCTTCAGc                              <  6:283296/148‑1 (MQ=255)
                                                                                                                             accac‑‑cTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCTGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACa                                                              <  6:22229/109‑1 (MQ=255)
                                                                                                                             accac‑‑cTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACa                                                              >  5:22229/1‑109 (MQ=255)
                                                                                                                               cac‑‑cTGTTCACGGGTTTCAACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAg                      >  4:246680/1‑147 (MQ=255)
                                                                                                                               cac‑‑cTGTTCACGGGTTTCAACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAg                      >  4:242009/1‑147 (MQ=255)
                                                                                                                               cac‑‑cTGTTCACGGGTTTCAACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAg                      <  3:246680/147‑1 (MQ=255)
                                                                                                                               cac‑‑cTGTTCACGGGTTTCAACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAg                      <  3:242009/147‑1 (MQ=255)
                                                                                                                                ac‑‑cTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTGCCCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAggg                    <  1:274535/149‑1 (MQ=255)
                                                                                                                                          aCGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAAt           >  3:103959/1‑149 (MQ=255)
                                                                                                                                            gggTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGaa            >  7:226992/1‑146 (MQ=255)
                                                                                                                                            gggTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGaa            <  8:226992/146‑1 (MQ=255)
                                                                                                                                            gggTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATg          >  8:147596/1‑148 (MQ=255)
                                                                                                                                              gTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAgg                     >  7:145762/1‑135 (MQ=255)
                                                                                                                                              gTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAgg                     <  8:145762/135‑1 (MQ=255)
                                                                                                                                              gTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATGCAg       <  3:335507/149‑1 (MQ=255)
                                                                                                                                               tttCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATGCAgg      >  5:323210/1‑149 (MQ=255)
                                                                                                                                                   cACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATGCAGGTggg  >  4:17604/1‑149 (MQ=255)
                                                                                                                                                   cACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATGCAGGTgg   <  7:147596/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAGTACACGCATGCCTTGTCCTCGATACCGTTGAAGGGCGGGTAAACCCACAGCTGCGCCAGCGACAGTGGCCGGCTGAAGTGGTAGCCCTGGGCAATATCGCAGCCGGCAAGCTT‑CAGGTACACCAC‑‑CTGTTCACGGGTTTCCACCCCCTCGGCCACCACTTCCAGGCCCAGGCGCTTGGCCAGGATGATGGTGGAGGACACGATTGGGCTG‑TCGTCATGGCTGTTGGACAGCGCCGCGATCAGCGAACGATCGATCTTCAGCTTGCTCAGGGGCAGTGAATGCAGGTGGG  >  NC_002947/3606852‑3607144

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