Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,301,427 G→C 8.3% intergenic (‑71/‑260) sbcD ← / → PP_5497 exonuclease subunit SbcD/rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,301,4270GC 8.3% 141.6 / 4.7 48intergenic (‑71/‑260)sbcD/PP_5497exonuclease subunit SbcD/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base C (2/2);  total (22/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 = 5.80e-01

CGTGTTCGAAGTCGCGTTCCTGGCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGC  >  NC_002947/2301284‑2301565
                                                                                                                                               |                                                                                                                                          
cGTGTTCGAAGTCGCGTTCCTGGCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCgg                                                                                                                                       >  7:122651/1‑149 (MQ=255)
cGTGTTCGAAGTCGCGTTCCTGGCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCgg                                                                                                                                       >  8:238219/1‑149 (MQ=255)
    ttCGAAGTCGCGTTCCTGGCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCACAGCGGcccc                                                                                                                                   <  8:161534/149‑2 (MQ=255)
             gcgTTCCTGGCCGTGCAGGCTTGAGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCggg                                                                                                                           >  7:40236/1‑148 (MQ=255)
                  ccTGGCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGc                                                                                                                      <  8:122651/148‑1 (MQ=255)
                     ctatctGCAGGCTTTGGCCCAGGGGCCAGGCGGAGGGATGAATCAGACGCATGGATGACCTGTTGGTTGTGAGGCGGGCAAGCGTGGAAAGGGCGGTATGGTCACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGcctg                                                                                                                   <  7:172875/143‑3 (MQ=255)
                      gCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCg                                                                                                                                        <  8:216980/126‑1 (MQ=255)
                      gCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCg                                                                                                                                        >  7:216980/1‑126 (MQ=255)
                       ccGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTt                                                                                                                <  4:180267‑M1/149‑2 (MQ=255)
                       ccGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTt                                                                                                                <  7:82523‑M1/149‑2 (MQ=255)
                         gTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtc                                                                                                               >  7:26932‑M1/1‑146 (MQ=255)
                            cAGGCTTTGGCCCAGGTGCCAGTCGGAGGTTTGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatc                                                                                                           <  8:96960‑M1/149‑7 (MQ=255)
                            cAGGCTTTGGCCCAGGTGCCAGTCGGAGGTAAGAAACAGACGCATGGATGACCTGTTGGAGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGGAACCCAAGCGTGGGGCCGCTGGGCGGCCCGAGCGCGCGccagccag                                                                                                                   >  8:172875/1‑141 (MQ=255)
                            aagGATTTGGCCCATGTGTCAGTCGGAGGTATGAAACAGACGCATGGATGAGCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctat                                                                                                            <  4:142006‑M1/147‑6 (MQ=255)
                                       ccAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaaga                                                                                                 >  3:11003‑M1/1‑132 (MQ=255)
                                              ccAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataaccca                                                                                          <  5:232894‑M1/148‑24 (MQ=255)
                                                aGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactg                                                                                       >  5:136585‑M1/1‑123 (MQ=255)
                                                  tCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATTGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgtt                                                                                     <  1:188521‑M1/149‑29 (MQ=255)
                                                   cGGAGGTATGAAACAGACGCATGGATGGCCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGAGTGGGGCCGCTGTGCTGCCCTATCGCTGGCAAGCCAGATtctatcaaacaaaagctaaaccactctt                                                                                     >  4:215381‑M1/1‑120 (MQ=255)
                                                    ggAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttct                                                                                   >  2:154973‑M1/1‑119 (MQ=255)
                                                         tATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttcta                                                                              <  7:238219‑M1/149‑36 (MQ=255)
                                                          aTGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctag                                                                             >  8:171437‑M1/1‑113 (MQ=255)
                                                           tGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctag                                                                             <  8:26932‑M1/148‑37 (MQ=255)
                                                                     gCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaa                                                                    >  3:7888‑M1/1‑102 (MQ=255)
                                                                      cttggtcgacctgttGGTGTCGAGTGGGGCAAGCGGGGAAATGGCGGTATGGTATCCCAAGGGTGGGGCCGCTGTGCGGCCCGATCGCGTGCAAGCCAGCTtctatcaaacaaaagataacgcactgttctttctagcgaaaaaaacag                                                                 <  7:42664‑M1/142‑49 (MQ=255)
                                                                         ggATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcg                                                              >  6:84255‑M1/1‑98 (MQ=255)
                                                                                   ttGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggt                                                      <  2:241008‑M1/147‑60 (MQ=255)
                                                                                       tgtCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatgg                                                >  5:224249‑M1/1‑84 (MQ=255)
                                                                                       tgtCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatgg                                                <  6:136585‑M1/149‑66 (MQ=255)
                                                                                              tctggcaaGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaa                                         <  4:7888‑M1/146‑73 (MQ=255)
                                                                                                   catgcggGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcgggggcaatggcggcgaagtgtt                                    <  3:215381‑M1/146‑78 (MQ=255)
                                                                                                    aaGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttc                                   <  8:40236‑M1/149‑79 (MQ=255)
                                                                                                       cGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgctcact                                <  7:171437‑M1/149‑82 (MQ=255)
                                                                                                        ggggAAAGGGCGGTATGGTAACCCAAGCGCGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggtttcaatggctgcgaagtgttcactc                               <  4:106210‑M1/149‑83 (MQ=255)
                                                                                                                 gCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttgtagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtggtcactcttgtatacg                      >  5:38234‑M1/1‑58 (MQ=255)
                                                                                                                 gCGGTATGGTAACCCAAGCGTGGGGCCGCTGAGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacg                      >  1:27762‑M1/1‑58 (MQ=255)
                                                                                                                   ggTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtg                    <  6:39821‑M1/149‑94 (MQ=255)
                                                                                                                    gTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaaca                                                                  >  3:95182‑M1/1‑55 (MQ=255)
                                                                                                                    gTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaaca                                                                  <  4:95182‑M1/102‑48 (MQ=255)
                                                                                                                         gTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtcac               <  4:109862‑M1/148‑99 (MQ=255)
                                                                                                                         gTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgtactttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtca                >  1:84783‑M1/1‑50 (MQ=255)
                                                                                                                            aCCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtcaatcaaacaaaagataacccacatttctttctagcgaaaaaaaaaggcgcgcggggtcaatggcggcgaagtgttcactcttctattcttggtcacaga            >  8:193995‑M1/1‑47 (MQ=255)
                                                                                                                             cccATGCGTGGGGCCGCTGTGCGTCCCGATCGCGGGCAAGCAAGCTtctatcaaacaaaagataacccactgttctttctagcgaataaaacaggcgcgcggggtcaatggaggcgaagtgttcactcttctatacgtggtgacagcaa          <  5:19796‑M1/149‑104 (MQ=255)
                                                                                                                              ccAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtcacagcaag         >  3:185623‑M1/1‑45 (MQ=255)
                                                                                                                                   gatcgggccgcACAGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtcacagcaagctgtg    >  8:111066‑M1/5‑40 (MQ=255)
                                                                                                                                    gTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtcacagcaagctgtgc   <  1:133788‑M1/149‑111 (MQ=255)
                                                                                                                                    atcgggccgCACAGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtcacagcaagctgtgc   <  7:111066‑M1/146‑111 (MQ=255)
                                                                                                                                     tcgggccGCACAGCGGCCCGATCGCGGGCAAGCCAGCTtctatcaaacaaaagataacccactgttctttctagcgaaaaaaacaggcgcgcggggtcaatggcggcgaagtgttcactcttctatacgtggtcacagcaagctgtgcc  >  6:240686‑M1/3‑38 (MQ=255)
                                                                                                                                               |                                                                                                                                          
CGTGTTCGAAGTCGCGTTCCTGGCCGTGCAGGCTTTGGCCCAGGTGCCAGTCGGAGGTATGAAACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGC  >  NC_002947/2301284‑2301565

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: