Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,852,597 G→C 9.1% R245R (CGC→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,5970GC 9.1% 130.7 / 4.8 44R245R (CGC→CGG) PP_2505GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base C (2/2);  total (23/21)
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.60e-01

ACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCGGTTCTGAAATGCT  >  NC_002947/2852452‑2852744
                                                                                                                                                 |                                                                                                                                                   
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                                                                              >  2:194556/1‑144 (MQ=255)
                                                                         ccGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACg                                                                              <  1:194556/144‑1 (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                                  >  4:31045/1‑149 (MQ=255)
                                                                                                                cagcagGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGctcct                                  >  1:299930/1‑149 (MQ=255)
                                                                                                                  gcagGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCa                                <  5:84381/149‑1 (MQ=255)
                                                                                                                                   gagtgcgCCCGCCGCCTGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCg               >  7:189631/7‑149 (MQ=255)
                                                                                                                                                ggcggcGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCGGTTCTGAAAtgct  <  5:176820/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCCTCTCCTCCCTGGCGTACCAGCAGGTCATCTGGCCGCAGCGTAGCCCGCAGGCGGCGGGCGCACTCGCAGAGCACTTTATCGCCGGCATCGTGGCCATGCTCGTCATTGACCGATTTGAAACGGTCGAGGTCTGCGTAGACGGCACCTAGCTGTAAACCTCGCTCCTCAGCCAAACGAAGTTCTCGGTTCTGAAATGCT  >  NC_002947/2852452‑2852744

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