Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,678,090 G→C 10.1% V365L (GTG→CTG)  PP_3242 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,678,0900GC10.1% 136.1 / 4.3 49V365L (GTG→CTG) PP_3242GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (25/19);  new base C (2/3);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTGCCGAGGTGCTGC  >  NC_002947/3677942‑3678237
                                                                                                                                                    |                                                                                                                                                    
gCTGGCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCg                                                                                                                                                      >  4:371487/1‑149 (MQ=255)
 cTGGCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTTAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTTACCGTCGGCGGCTGGATGGCgt                                                                                                                                                     <  7:236390/149‑1 (MQ=255)
 cTGGCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCgt                                                                                                                                                     <  7:154500/149‑1 (MQ=255)
  tGGCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACAAGACGCTGAAAGCGTTGGACTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCgtg                                                                                                                                                    >  5:336235/1‑149 (MQ=255)
    gCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCCCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGtt                                                                                                                                                  >  1:204166/1‑149 (MQ=255)
                ggTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACt                                                                                                                                      >  3:72866/1‑149 (MQ=255)
                      gtgCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGcc                                                                                                                                 <  3:371487/148‑1 (MQ=255)
                       tgctgcGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGt                                                                                                                               >  8:225752/1‑149 (MQ=255)
                           gcgcCAACTGCGCGAGCGGAGCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGGTGGATGGCGTGTGCGATCCTGCACTGCGCCGAGc                                                                                                                             >  6:5396/1‑147 (MQ=255)
                            cgcCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACg                                                                                                                          >  1:317269/1‑149 (MQ=255)
                                  cTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGgcgcag                                                                                                                    <  1:200212/149‑3 (MQ=255)
                                   tGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGcc                                                                                                                     >  6:281232/1‑147 (MQ=255)
                                     cgcgAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCCTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAgg                                                                                                                 >  4:23034/1‑149 (MQ=255)
                                     cgcgAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCCTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAgg                                                                                                                 >  4:27992/1‑149 (MQ=255)
                                        gAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTa                                                                                                              >  6:71856/1‑149 (MQ=255)
                                               tCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTg                                                                                                       >  1:13936/1‑149 (MQ=255)
                                                   gATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCAc                                                                                                   >  7:123694/1‑149 (MQ=255)
                                                       gAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGt                                                                                               <  8:123694/149‑1 (MQ=255)
                                                         gCGGGCGCTGCTCCAGGCGCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGa                                                                                             <  6:278971/149‑1 (MQ=255)
                                                            ggCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAAGCTATCCCCTCGCACTGCTGATGc                                                                                          >  4:32148/1‑149 (MQ=255)
                                                                   ctCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGctggac                                                                                     <  2:317269/147‑1 (MQ=255)
                                                                      cAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGAGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGATGATGctggacctg                                                                                  <  4:72866/147‑1 (MQ=255)
                                                                             aCCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCCTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTc                                                                        <  3:27992/149‑1 (MQ=255)
                                                                             aCCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCCTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTc                                                                        <  3:23034/149‑1 (MQ=255)
                                                                             aCCAGACGCTGAAAGCGTTGGCCTTGACAGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACtt                                                                         >  1:21572/1‑148 (MQ=255)
                                                                              ccAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTc                                                                        >  6:113001/1‑148 (MQ=255)
                                                                               cAGACGCTGAAAGGGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAgg                                                                                                                 >  6:163905/1‑107 (MQ=255)
                                                                               cAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAgg                                                                                                                 <  5:163905/107‑1 (MQ=255)
                                                                                aGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAg                                                                     >  2:349192/1‑149 (MQ=255)
                                                                                 gACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGCTGATGCTGGACCTGGACTACTTCAAgg                                                                    >  6:291828/1‑149 (MQ=255)
                                                                                    gCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTCGACCTGGACTACTTCAAGGcc                                                                  >  2:21860/1‑148 (MQ=255)
                                                                                      tGAAAGGGTTGGCCTTGACGGACAGCCTTTCAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCTTGCACTGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGTACTACTTCAAGGCGTAc                                                               <  6:336235/149‑1 (MQ=255)
                                                                                      tGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTa                                                                >  1:335046/1‑148 (MQ=255)
                                                                                       gAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTAc                                                               >  2:102157/1‑148 (MQ=255)
                                                                                         aaGCGTTGGCCTTGACCCACAGCCTTACCGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAAc                                                            <  2:373142/149‑1 (MQ=255)
                                                                                                   cTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACg                                                  <  1:171977/149‑1 (MQ=255)
                                                                                                               ccTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTg                                      >  5:146335/1‑149 (MQ=255)
                                                                                                                 ttACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGc                                    >  3:310202/1‑149 (MQ=255)
                                                                                                                   tcAGGGCTGGGTAACCGTCGGCGGGTGGATGGCGTGGTCGATCCTGAAGTGCGCCGTGCACGGCGCGATGGCTATGCCCTGGCAGTGGTGTTGCTGGACCTGGACTACTTCAAGGCCTTCTAGGACCTCTACGGTCACCCGGCTGGCg                                   <  8:347925/147‑1 (MQ=255)
                                                                                                                     aGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCg                                   >  7:347925/1‑146 (MQ=255)
                                                                                                                     aGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGAc                                 >  3:327121/1‑148 (MQ=255)
                                                                                                                      gggCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCa                               <  5:281232/149‑1 (MQ=255)
                                                                                                                        gCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGt                             <  2:335046/149‑1 (MQ=255)
                                                                                                                         cTGGGTAACCGTCGGCGGCTGGATGGCCTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTg                            <  5:5396/149‑1 (MQ=255)
                                                                                                                             gTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTg                        <  2:204166/149‑1 (MQ=255)
                                                                                                                              tAACCGTCGGCGGCTGGATGGCGTGTTCGACCCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCAGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCAACAACGACCACTACGGTCACCCCGCTGGCGCCCAGTGCTAGc                       >  8:146970/1‑149 (MQ=255)
                                                                                                                              tAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGTACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTg                        <  1:102157/149‑1 (MQ=255)
                                                                                                                                   gTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGGGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGc                  <  6:297911/149‑1 (MQ=255)
                                                                                                                                       gcggcTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGAGCTGGACTACTTCAACGCCTAAAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTg              >  7:128557/1‑149 (MQ=255)
                                                                                                                                              gATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTGCCGAGGt       >  1:80142/1‑149 (MQ=255)
                                                                                                                                              gATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTGCCGAGGt       <  4:310202/149‑1 (MQ=255)
                                                                                                                                                  gCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTGCCGAGGtgctg   <  5:113001/149‑1 (MQ=255)
                                                                                                                                                   cGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTGCCGAGGtgctgc  <  1:349192/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGGCAGCGTTTGGCGGTGTGGTGCTGCGCCAACTGCGCGAGCGGATCGTGATCGAGCGGGCGCTGCTCCAGGCTCACCAGACGCTGAAAGCGTTGGCCTTGACCGACAGCCTTACAGGGCTGGGTAACCGTCGGCGGCTGGATGGCGTGTTCGAACCTGAACTGCGCCGTGCACGGCGCCAAGGCTATCCCCTGGCACTGGTGATGCTGGACCTGGACTACTTCAAGGCCTACAACGACCACTACGGTCACCCGGCTGGCGACCAGTGCTTGCGACGCTTTGCCGAGGTGCTGC  >  NC_002947/3677942‑3678237

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: