Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,911,965 T→G 8.9% L104R (CTG→CGG) ‡ PP_3452 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,911,9650TG 8.9% 128.7 / 3.8 45L104R (CTG→CGG) ‡PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (17/24);  new base G (2/2);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

TCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCT  >  NC_002947/3911818‑3912103
                                                                                                                                                   |                                                                                                                                           
tCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCATCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTg                                                                                                                                            >  2:117964/1‑149 (MQ=255)
tCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTg                                                                                                                                            >  2:47699/1‑149 (MQ=255)
  ggCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTgg                                                                                                                                           <  5:3570/148‑1 (MQ=255)
  ggCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGt                                                                                                                                          <  8:128341/149‑1 (MQ=255)
         aaGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGaa                                                                                                                                  <  4:195729/149‑1 (MQ=255)
              ttCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCAt                                                                                                                             >  8:289857/1‑149 (MQ=255)
                          tccttcctGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCt                                                                                                                 <  1:117964/149‑1 (MQ=255)
                           ccttcctGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCt                                                                                                                 <  1:47699/148‑1 (MQ=255)
                            cttcctGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTgc                                                                                                               >  7:223476/1‑149 (MQ=255)
                             ttcctGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTgcg                                                                                                              >  2:268827/1‑149 (MQ=255)
                                 tGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCg                                                                                                          >  6:192700/1‑149 (MQ=255)
                                     cATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGgtgt                                                                                                      >  2:234902/1‑149 (MQ=255)
                                     cATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGgtgt                                                                                                      >  2:4098/1‑149 (MQ=255)
                                         gCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGtggt                                                                                                   >  4:62257/1‑148 (MQ=255)
                                         gCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGtggt                                                                                                   <  3:62257/148‑1 (MQ=255)
                                          cccTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATc                                                                                                                            >  6:76480/1‑122 (MQ=255)
                                          cccTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATc                                                                                                                            <  5:76480/122‑1 (MQ=255)
                                             tGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTg                                                                                              >  3:172444/1‑149 (MQ=255)
                                              gCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGc                                                                                             >  4:132715/1‑149 (MQ=255)
                                                  ggCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCgg                                                                                         <  2:147730/149‑1 (MQ=255)
                                                   gctcccTGAGCGTAGGGCTGCAGTTCACCCGCCGCAGCCATCACTTCGACGGGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCggg                                                                                        <  4:259074/145‑1 (MQ=255)
                                                   gCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCggg                                                                                        <  7:241277/149‑1 (MQ=255)
                                                              gTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTAGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGATCCGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCAcc                                                                               <  1:234902/148‑1 (MQ=255)
                                                                     tgcATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTAt                                                                      <  8:223476/149‑1 (MQ=255)
                                                                       cATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAAGCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTAt                                                                      <  6:228915/147‑1 (MQ=255)
                                                                              cccGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCg                                                             <  1:55493/149‑1 (MQ=255)
                                                                                     cAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTTGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTAc                                                      <  1:268827/149‑1 (MQ=255)
                                                                                           tCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCTCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTgc                                                >  2:82620/1‑149 (MQ=255)
                                                                                                   aCGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTgg                                        <  7:289857/149‑1 (MQ=255)
                                                                                                         tCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCAttc                                  >  4:157765/1‑149 (MQ=255)
                                                                                                           aaCCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCAttctt                                <  5:192700/149‑1 (MQ=255)
                                                                                                                    ttCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCAt                       >  3:301300/1‑149 (MQ=255)
                                                                                                                    ttCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCAt                       >  1:207853/1‑149 (MQ=255)
                                                                                                                       aTTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGt                    <  3:31583/149‑1 (MQ=255)
                                                                                                                         gccttcgatcaCCAGGCGCATGCCCAGGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGtc                   >  7:142772/8‑148 (MQ=255)
                                                                                                                         gccttcgatcaCCAGGCGCATGCCCAGGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGtc                   <  8:142772/141‑1 (MQ=255)
                                                                                                                           ccGTCATCACCCTGGGAATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGtcatc                <  3:157765/149‑1 (MQ=255)
                                                                                                                            ttcgatcaCCAGGCGCATGCCCAGGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGtcatca               >  5:94383/5‑149 (MQ=255)
                                                                                                                                 tcaCCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTgg          <  4:301300/149‑1 (MQ=255)
                                                                                                                                   aCCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGt        <  3:132715/149‑1 (MQ=255)
                                                                                                                                         ggCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCTCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGc   <  1:82620/148‑1 (MQ=255)
                                                                                                                                         ggCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATg    <  2:207853/147‑1 (MQ=255)
                                                                                                                                         cgCATGCCCAGGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCt  <  6:94383/148‑1 (MQ=255)
                                                                                                                                            aTGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGtca                  <  4:180116/130‑1 (MQ=255)
                                                                                                                                            aTGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGtca                  >  3:180116/1‑130 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCT  >  NC_002947/3911818‑3912103

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