Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,911,964 C→A 8.9% L104M (CTG→ATG) ‡ 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,9640CA 8.9% 129.6 / 3.9 45L104M (CTG→ATG) ‡PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (17/24);  new base A (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 = 4.94e-01

CCTCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCT  >  NC_002947/3911816‑3912103
                                                                                                                                                    |                                                                                                                                            
ccTCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGcc                                                                                                                                              >  2:55493/1‑149 (MQ=255)
  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                                                                                                                            <  5:76480/122‑1 (MQ=255)
                                            cccTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATc                                                                                                                            >  6:76480/1‑122 (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)
                                                                                                                                                    |                                                                                                                                            
CCTCGGCGACCAAGGCTTCACCTGGCTGTCCTTCCTGTTCATCGCCCTGCTCGGCAGCCTGATCGTAGTGCTGCATTTCACCCGCCGCAGCCATCACTTCGACGTGCTCAACCTGCTGTTCATTACCGTCATCACCCTGGGCATGCGCCTGGTGAT‑CGAAGGCATCCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCT  >  NC_002947/3911816‑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: ‑