Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,755 G→A 9.3% G396G (GGC→GGT)  PP_2413 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,759,7550GA 9.3% 118.8 / 2.0 43G396G (GGC→GGT) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (2/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGGCGCAGTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGC  >  NC_002947/2759607‑2759901
                                                                                                                                                    |                                                                                                                                                   
gTTGGCGCAGTACCTGGTCTCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAgg                                                                                                                                                     <  3:232341/149‑1 (MQ=255)
         gTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGc                                                                                                                                            <  2:91942/149‑1 (MQ=255)
         gTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTCGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGc                                                                                                                                            <  5:347859/149‑1 (MQ=255)
            ccTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGa                                                                                                                                         <  3:65702/149‑1 (MQ=255)
               ggTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCg                                                                                                                                      >  5:200951/1‑149 (MQ=255)
                gTCGCCCACGTCGTGTCCATGGCTATGATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCg                                                                                                                                      <  2:368279/148‑1 (MQ=255)
                    cccACGTAGTGTCTATGGCTATAATTTATCCGCTTGAAGTGGTCGATATCCAGCTAAATGGCTGAAAACTCCCTCTTTTAAGCCTATAGTAACGAGAGGGTGATTTCCATGCTGCGGCGGTTGCCCAGAGCGATCCGCTGACCGGCCTg                                                                                                                                 <  2:111436/149‑1 (MQ=255)
                     ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACTGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                                <  2:268680/149‑1 (MQ=255)
                      cACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGACCGATCCGCTGACCGGCCTggg                                                                                                                               <  4:217634/149‑1 (MQ=255)
                       aCGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGc                                                                                                                              <  2:218995/149‑1 (MQ=255)
                              gtCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGa                                                                                                                       >  1:28934/1‑149 (MQ=255)
                               tCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAt                                                                                                                      <  8:244089/149‑1 (MQ=255)
                                     gCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgc                                                                                                                 <  2:28934/148‑1 (MQ=255)
                                              tATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATa                                                                                                        <  6:200951/148‑1 (MQ=255)
                                                tCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGc                                                                                                     >  5:42979/1‑149 (MQ=255)
                                                 ccGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAAACTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGct                                                                                                    >  1:117341/1‑149 (MQ=255)
                                                        aaGTGGTCGATACCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCAGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTgca                                                                                              <  2:44808/148‑1 (MQ=255)
                                                             gTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAg                                                                                         <  2:145097/148‑1 (MQ=255)
                                                              tggATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGGCCAGGCCGGTCAGCGGATCGGTCTGGGCGTGCCGATTGAGCCGTCCGATACGCTCTTGCAGCAAGtt                                                                                       <  4:199950/147‑1 (MQ=255)
                                                                      cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTg                                                                                >  6:214761/1‑148 (MQ=255)
                                                                                      aaCTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCtt                                                              <  6:42979/149‑1 (MQ=255)
                                                                                         tCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAg                                                           >  2:50898/1‑149 (MQ=255)
                                                                                                    aGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGgccgcc                                                 <  3:191242/148‑1 (MQ=255)
                                                                                                         ccAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAg                                           >  2:41481/1‑149 (MQ=255)
                                                                                                          cAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGaa                                            <  7:327314/147‑1 (MQ=255)
                                                                                                           aGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCAGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTa                                         >  8:260161/1‑149 (MQ=255)
                                                                                                            gTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATCGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTAc                                        >  6:58150/1‑149 (MQ=255)
                                                                                                               aCGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCATACGCGCGCTTGAGCTCGGCCGCCTCGAAGTAcc                                       <  1:50898/148‑1 (MQ=255)
                                                                                                                    aGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCg                                >  5:367311/1‑149 (MQ=255)
                                                                                                                    aGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCg                                >  5:38015/1‑149 (MQ=255)
                                                                                                                       cTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCAc                             >  3:324499/1‑149 (MQ=255)
                                                                                                                       cTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCAc                             <  6:367311/149‑1 (MQ=255)
                                                                                                                         gATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGCACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCt                           >  8:310079/1‑149 (MQ=255)
                                                                                                                             tCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATg                       >  3:322487/1‑149 (MQ=255)
                                                                                                                             tCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATg                       >  6:200454/1‑149 (MQ=255)
                                                                                                                             tCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATg                       >  7:21944/1‑149 (MQ=255)
                                                                                                                                aGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCTGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCa                     >  1:293616/1‑148 (MQ=255)
                                                                                                                                aGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCa                     >  3:297695/1‑148 (MQ=255)
                                                                                                                                            ttGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTcggcgg         >  8:350682/1‑148 (MQ=255)
                                                                                                                                              gCCCAGACCGATCCGCTGACCGGCCTGCGCGTCCCGATTGACCCGGCCCATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGACCTTGAGCTCCGCCGCCTCGAACTACCAGGACCGCACCTGATGCAAATGCTCGGCGGTAc      >  2:92678/1‑149 (MQ=255)
                                                                                                                                               cccAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTAc      <  7:280519/148‑1 (MQ=255)
                                                                                                                                                ccAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCg    >  7:233504/1‑149 (MQ=255)
                                                                                                                                                ccAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCg    >  3:354158/1‑149 (MQ=255)
                                                                                                                                                  aGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGgc  >  7:234374/1‑149 (MQ=255)
                                                                                                                                                    gCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATg                 >  2:286550/1‑132 (MQ=255)
                                                                                                                                                    gCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATg                 <  1:286550/132‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGGCGCAGTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGC  >  NC_002947/2759607‑2759901

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑