Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,852,595 A→C 8.5% L246R (CTG→CGG)  PP_2505 ← GAF domain/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,852,5950AC 8.5% 141.3 / 4.3 47L246R (CTG→CGG) PP_2505GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base C (2/2);  total (24/23)
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.94e-01

GGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCG  >  NC_002947/2852447‑2852731
                                                                                                                                                    |                                                                                                                                        
gggTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCa                                                                                                                                          <  4:259670/149‑1 (MQ=255)
  gTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAga                                                                                                                                        <  6:210492/149‑2 (MQ=255)
     aCAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAggcgg                                                                                                                                     >  5:287583/1‑149 (MQ=255)
        gAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGgcggcgn                                                                                                                                  <  2:177925/149‑2 (MQ=255)
        gAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGcggcgg                                                                                                                                  <  8:192448/149‑1 (MQ=255)
             ggCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCCGCCtgcgggctac                                                                                                                             <  2:147972/149‑4 (MQ=255)
                gcgcTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTc                                                                                                                          <  5:303375/149‑1 (MQ=255)
                  gcTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCg                                                                                                                         >  8:141577/1‑148 (MQ=255)
                   cTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCa                                                                                                                       <  4:220221/149‑1 (MQ=255)
                      ccGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCagag                                                                                                                    >  3:295695/1‑149 (MQ=255)
                               tgacGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCCGCCtgcgggctacg                                                                                                                            <  8:296928/132‑5 (MQ=255)
                               tgacGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCCGCCtgcgggctacg                                                                                                                            >  7:296928/1‑128 (MQ=255)
                               tgacGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATc                                                                                                           >  6:60018/1‑149 (MQ=255)
                                 acGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATc                                                                                                           <  6:98710/147‑1 (MQ=255)
                                   gCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCg                                                                                                       >  6:216279/1‑149 (MQ=255)
                                   gCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCg                                                                                                       <  5:60018/149‑1 (MQ=255)
                                   gCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCg                                                                                                       >  7:40521/1‑149 (MQ=255)
                                             cGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTg                                                                                               <  6:287583/147‑1 (MQ=255)
                                                  tCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGc                                                                                        >  7:171564/1‑149 (MQ=255)
                                                      cgcCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGc                                                                                        <  4:295695/145‑1 (MQ=255)
                                                      cgcCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTc                                                                                      <  8:40521/147‑1 (MQ=255)
                                                       gcCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACtt                                                                                                               <  1:118592/121‑1 (MQ=255)
                                                       gcCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACtt                                                                                                               >  2:118592/1‑121 (MQ=255)
                                                        ccAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCa                                                                                  >  6:176820/1‑149 (MQ=255)
                                                              cTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGAc                                                                             <  5:216279/148‑1 (MQ=255)
                                                                      tCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTg                                                                      <  8:171564/147‑1 (MQ=255)
                                                                            aCCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTc                                                              >  5:15911/1‑149 (MQ=255)
                                                                             cccGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGCGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTc                                                              >  3:204690/1‑148 (MQ=255)
                                                                              ccGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACg                                                                 <  1:194556/144‑1 (MQ=255)
                                                                              ccGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACg                                                                 >  2:194556/1‑144 (MQ=255)
                                                                               cGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAg                                                           <  7:24565/149‑1 (MQ=255)
                                                                                   gCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCt                                                       >  6:275412/1‑149 (MQ=255)
                                                                                     aaGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGc                                                     <  5:275412/149‑1 (MQ=255)
                                                                                          tACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCTGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGa                                                >  8:346073/1‑149 (MQ=255)
                                                                                            cGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGa                                                <  6:15911/147‑1 (MQ=255)
                                                                                                   ctctcctcCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCAcc                                         >  7:125345/1‑147 (MQ=255)
                                                                                                      tcctccCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCt                                    >  2:219920/1‑149 (MQ=255)
                                                                                                       cctccCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCt                                    >  2:219508/1‑148 (MQ=255)
                                                                                                       cctccCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCt                                    <  1:219508/148‑1 (MQ=255)
                                                                                                           ccTGGCTTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTaaa                               <  7:141577/149‑1 (MQ=255)
                                                                                                           ccTGGCGTACCAGCAGGTCATCTGCCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTaaa                               >  5:160385/1‑149 (MQ=255)
                                                                                                               gCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGCGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAAcc                             <  4:204690/147‑1 (MQ=255)
                                                                                                               gCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTc                           >  2:94663/1‑149 (MQ=255)
                                                                                                                     cagcagGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCCTTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGctcct                     >  6:334365/1‑149 (MQ=255)
                                                                                                                     cagcagGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGctcct                     >  1:299930/1‑149 (MQ=255)
                                                                                                                     cagcagGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGctcct                     >  4:31045/1‑149 (MQ=255)
                                                                                                                       gcagGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCa                   <  5:84381/149‑1 (MQ=255)
                                                                                                                                        gagtgcgCCCGCCGCCTGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCg  >  7:189631/7‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
GGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCG  >  NC_002947/2852447‑2852731

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