Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,772 T→C 9.5% T391A (ACC→GCC)  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,7720TC 9.5% 110.7 / 3.1 42T391A (ACC→GCC) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (22/16);  new base C (2/2);  total (24/18)
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.23e-01

CCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGC  >  NC_002947/2759627‑2759919
                                                                                                                                                 |                                                                                                                                                    
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                                         >  7:21944/1‑149 (MQ=255)
                                                                                                         tCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATg                                         >  6:200454/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)
                                                                                                                                          ggATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGt        <  1:41481/149‑1 (MQ=255)
                                                                                                                                                gTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCGTGCAGCAAGTTCAGGCTGAACCGCA‑ACGCGCGCTGGAGCTCGGCCGCCGCGAAGTACCCGGCCCGCACCGGAGGCAAATGCTCGGCGGTACCGGCGCGGTCCAGCGTGCGGGc  >  4:286452/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACAGGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCA‑ACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGC  >  NC_002947/2759627‑2759919

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