Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,661 G→C 10.8% intergenic (+27/+94) PP_3662 → / ← PP_3663 hypothetical protein/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,6610GC10.8% 103.2 / 4.7 37intergenic (+27/+94)PP_3662/PP_3663hypothetical protein/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base C (2/2);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

CACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCC  >  NC_002947/4160515‑4160809
                                                                                                                                                  |                                                                                                                                                    
cACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGa                                                                                                                                                    <  6:31129/149‑1 (MQ=255)
   ggTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGc                                                                                                                                                 <  2:217462/149‑1 (MQ=255)
       aCAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCg                                                                                                                                              <  2:42976/148‑1 (MQ=255)
        cAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTc                                                                                                                                            <  1:142050/149‑1 (MQ=255)
                  cTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCgggg                                                                                                                                  <  6:207974/149‑1 (MQ=255)
                     gACCCGCTGGGCCGGCTACTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCggggg                                                                                                                                 <  7:154226/147‑1 (MQ=255)
                          gCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCCGCGCGGCGTATCTACCGTGCTATCAGGTGGCATCCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTg                                                                                                                          >  5:62591/1‑149 (MQ=255)
                                   gctgctGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCggg                                                                                                                                   >  3:262933/1‑131 (MQ=255)
                                   gctgctGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCggg                                                                                                                                   <  4:262933/131‑1 (MQ=255)
                                        tgcAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTg                                                                                                             >  8:64418/1‑148 (MQ=255)
                                                ttCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCa                                                                                                    >  8:68962/1‑149 (MQ=255)
                                                ttCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGt                                                                                                                                             >  8:218062/1‑108 (MQ=255)
                                                ttCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGt                                                                                                                                             <  7:218062/108‑1 (MQ=255)
                                                   gtcgACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCAccc                                                                                                 <  8:101914/149‑1 (MQ=255)
                                                                    gAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGa                                                                                  >  6:47369/1‑147 (MQ=255)
                                                                     aaGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGaaa                                                                                <  7:68962/148‑1 (MQ=255)
                                                                            ccGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTcct                                                                                                                                          >  8:252800/1‑81 (MQ=255)
                                                                            ccGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTcct                                                                                                                                          <  7:252800/83‑3 (MQ=255)
                                                                             cggcggCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCca                                                                                                <  3:177409/124‑1 (MQ=255)
                                                                             cggcggCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCca                                                                                                >  4:177409/1‑124 (MQ=255)
                                                                                  gcgcggcgTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAg                                                                  <  7:64418/149‑1 (MQ=255)
                                                                                        cgTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGCGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACt                                                            >  1:84734/1‑149 (MQ=255)
                                                                                            tGTACGGTGCGATGAGGGCGCCACCTGACCCGCTGCTTTACTGCGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACGCACATGGGTTACATGAAAGGTCCCCCACCTAGGTGACTCGCt                                                        <  6:62591/149‑1 (MQ=255)
                                                                                               aCCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCAc                                                     >  4:52617/1‑149 (MQ=255)
                                                                                               aCCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCAc                                                     >  4:35362/1‑149 (MQ=255)
                                                                                                      tATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCg                                                >  6:251979/1‑147 (MQ=255)
                                                                                                          aGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGtt                                           <  3:35362/148‑1 (MQ=255)
                                                                                                             tGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGc                                                  >  4:179020/1‑138 (MQ=255)
                                                                                                             tGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGc                                                  <  3:179020/138‑1 (MQ=255)
                                                                                                                 aaCCTGACCCGCTGCTTTACTGAGGACGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCaa                                        >  6:264576/1‑143 (MQ=255)
                                                                                                                 aaCCTGACCCGCTGCTTTACTGAGGACGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCaa                                        <  5:264576/143‑1 (MQ=255)
                                                                                                                         ccGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGAc                                                             <  4:257222/115‑1 (MQ=255)
                                                                                                                         ccGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGAc                                                             >  3:257222/1‑115 (MQ=255)
                                                                                                                         ccGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTc                           >  5:128291/1‑149 (MQ=255)
                                                                                                                                   aCTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCGGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCCACTAGGTGACTCGCTCACTGCCGGAGGTCAACAACCGCCCGCTCGAAAAACCg                  >  2:91472/1‑148 (MQ=255)
                                                                                                                                              gCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCAccgc                                 <  6:227496/122‑1 (MQ=255)
                                                                                                                                              gCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCAccgc                                 >  5:227496/1‑122 (MQ=255)
                                                                                                                                                  ggATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGcccc  >  2:249715/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCC  >  NC_002947/4160515‑4160809

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: