Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,009,531 C→T 8.9% R441H (CGC→CAC)  PP_5263 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,009,5310CT 8.9% 125.3 / 3.2 45R441H (CGC→CAC) PP_5263GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (2/2);  total (23/22)
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.83e-01

CCAGCTTGCGGTTCTCTTCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAGCTGGCTGTCCTG  >  NC_002947/6009388‑6009673
                                                                                                                                               |                                                                                                                                              
ccAGCTTGCGGTTCTCTTCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGAcc                                                                                                                                            <  8:28634/148‑1 (MQ=255)
 cAGCTTGCGGTTCTCTTCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGAcc                                                                                                                                            >  7:192687/1‑147 (MQ=255)
             ctctTCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATgg                                                                                                                              >  6:7078/1‑149 (MQ=255)
               ctTCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTc                                                                                                                            >  1:126538/1‑149 (MQ=255)
                ttCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACTAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGt                                                                                                                             >  1:271017/1‑147 (MQ=255)
                    cgcATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGc                                                                                                                                                <  8:52264/124‑1 (MQ=255)
                    cgcATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGc                                                                                                                                                >  7:52264/1‑124 (MQ=255)
                        tttccatttccGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGATTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGcac                                                                                                                    >  5:270496/1‑148 (MQ=255)
                             atttccGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTc                                                                                                              <  2:126538/149‑1 (MQ=255)
                             atttccGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTc                                                                                                              >  2:237292/1‑149 (MQ=255)
                                   gCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCTCGTCCATGGTCTGGCGCACCACCTCGTTgcg                                                                                                        >  8:100263/1‑149 (MQ=255)
                                       caaAGCTGCGGTCCACCTTGATCAAGTTGATCGGCAGGCTGTTTAGGTGCAAGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGgcgc                                                                                                    <  6:270496/149‑1 (MQ=255)
                                       caaAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGgcgc                                                                                                    <  5:7078/149‑1 (MQ=255)
                                            cTGAGGACCAACTTGAGCAAGGTGATCTGCAGGCTGTTTAGGTGCAAGAACGAATAAAACCCGGTGCCGAAATCATCGATCGAGAAGCGAACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGTCGCATCAc                                                                                               <  5:176500/149‑1 (MQ=255)
                                               cGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCAccg                                                                                             <  7:100263/148‑1 (MQ=255)
                                                    cACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGtct                                                                                       >  1:24847/1‑149 (MQ=255)
                                                      ccTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCTTCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCg                                                                                     <  7:28589/149‑1 (MQ=255)
                                                      ccTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCg                                                                                     <  5:127455/149‑1 (MQ=255)
                                                          gAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGt                                                                                   <  1:237292/147‑1 (MQ=255)
                                                             cAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCt                                                                              >  1:266503/1‑149 (MQ=255)
                                                                  tGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGTGACCCAGCGCGTCCATGGTCTgg                                                                                                                         <  6:116726/101‑1 (MQ=255)
                                                                  tGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGTGACCCAGCGCGTCCATGGTCTgg                                                                                                                         >  5:116726/1‑101 (MQ=255)
                                                                      cGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCg                                                                                     <  7:5143/133‑1 (MQ=255)
                                                                      cGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCg                                                                                     >  8:5143/1‑133 (MQ=255)
                                                                       ggCAGGCTGTTTAGGTGCACGAACAAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCaa                                                                     <  2:271017/148‑1 (MQ=255)
                                                                                        cACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGTGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCacac                                                   >  3:92178/1‑149 (MQ=255)
                                                                                           gaacgaaGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGt                                                <  3:9688/149‑1 (MQ=255)
                                                                                                 agaaAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTc                                            >  8:40030/1‑147 (MQ=255)
                                                                                                         cGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGt                                                                                   >  6:256916/1‑100 (MQ=255)
                                                                                                         cGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGt                                                                                   <  5:256916/100‑1 (MQ=255)
                                                                                                            tGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGTCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGa                               <  7:40030/149‑1 (MQ=255)
                                                                                                                 aaaTCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCaa                          >  1:226696/1‑149 (MQ=255)
                                                                                                                  aaTCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAAc                         >  2:63891/1‑149 (MQ=255)
                                                                                                                  aaTCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAAc                         >  2:60228/1‑149 (MQ=255)
                                                                                                                    tcatcGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGt                       <  2:266503/149‑1 (MQ=255)
                                                                                                                     catcGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTg                      >  3:233814/1‑149 (MQ=255)
                                                                                                                          agcgagAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCTCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAg                 <  8:192687/149‑1 (MQ=255)
                                                                                                                          agcgagAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGaa                                                                          <  4:224593/92‑1 (MQ=255)
                                                                                                                          agcgagAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGaa                                                                          >  3:224593/1‑92 (MQ=255)
                                                                                                                                 aGCGCACGCCCAGGTGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAgctgg          <  4:92178/149‑1 (MQ=255)
                                                                                                                                   cgcACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAgctggct        <  7:149042/149‑1 (MQ=255)
                                                                                                                                   cgcACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAgctgg          <  2:157705/147‑1 (MQ=255)
                                                                                                                                   cgcACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAgctgg          >  1:157705/1‑147 (MQ=255)
                                                                                                                                          ccAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGgcgc                                                                                                    >  2:217655/1‑50 (MQ=255)
                                                                                                                                          ccAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGgcgc                                                                                                    <  1:217655/50‑1 (MQ=255)
                                                                                                                                          ccAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAGCTGGCTGTCCTg  >  2:232553/1‑148 (MQ=255)
                                                                                                                                          ccAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAGCTGGCTGTCCTg  >  1:4508/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CCAGCTTGCGGTTCTCTTCGCGCATTTCCATTTCCGCAACAAAGCTGCGGTCCACCTTGAGCAAGGTGATCGGCAGGCTGTTTAGGTGCACGAACGAAGAAAACCCGGTGCCGAAATCATCGAGCGAGAAGCGCACGCCCAGGCGACCCAGCGCGTCCATGGTCTGGCGCACCACCTCGTTGCGGCGCATCACCGCCGTCTCGGTCAGCTCGAACTCAAGCCAGCGGGCATCCACACCGTGCTCGAGGATCAGGCGACTCAACGTGGCCAGTAGCTGGCTGTCCTG  >  NC_002947/6009388‑6009673

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