Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,068 G→T 15.6% R171R (CGC→CGA PP_4277 ← GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,866,0680GT15.6% 79.5 / 5.9 32R171R (CGC→CGAPP_4277GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base T (2/3);  total (14/18)
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.31e-01

GCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGT  >  NC_002947/4865924‑4866207
                                                                                                                                                |                                                                                                                                           
gCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTg                                                                                                                                        <  1:698686/150‑1 (MQ=255)
gCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTAGATCATGTTGGCCTGGCGCTg                                                                                                                                        <  2:621739/150‑1 (MQ=255)
        ccAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCg                                                                                                                                >  1:753891/1‑150 (MQ=255)
            tGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCGCCAg                                                                                                                            <  2:753891/150‑1 (MQ=255)
            tGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCCATCATGTTGGCCTGTCGCTGGACCAGCGCCAg                                                                                                                            <  1:104521/150‑1 (MQ=255)
                 ccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGTACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAAt                                                                                                                       <  1:634749/150‑1 (MQ=255)
                   aGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAg                                                                                                                     <  2:604657/150‑1 (MQ=255)
                                            cATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTc                                                                                                                                     >  2:722562/1‑109 (MQ=255)
                                            cATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTc                                                                                                                                     <  1:722562/109‑1 (MQ=255)
                                               aGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGTTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGtttt                                                                                         >  1:111872/1‑150 (MQ=255)
                                                   ggTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCAATCATGTTGGCCTGgcgc                                                                                                                                          >  2:741727/1‑97 (MQ=255)
                                                   ggTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCAATCATGTTGGCCTGgcgc                                                                                                                                          <  1:741727/97‑1 (MQ=255)
                                                              cgcgAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCa                                                                          <  1:526437/150‑1 (MQ=255)
                                                                               tCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                               <  1:490175/128‑1 (MQ=255)
                                                                               tCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                               >  2:490175/1‑128 (MQ=255)
                                                                                    cgAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACgg                                                    <  1:358692/150‑1 (MQ=255)
                                                                                                   gAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCgg                                                                                                    >  1:644440/1‑87 (MQ=255)
                                                                                                   gAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCgg                                                                                                    <  2:644440/87‑1 (MQ=255)
                                                                                                    aGCAAGCGGGAGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATg                                    >  1:530567/1‑150 (MQ=255)
                                                                                                                     gTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCCGGTTTTCCACGGTTTCCTGCAGGAACGGGATGTAGGCGAACGGATAAGTCGGACGCAGGACCGCAAAACTCTCgcc                   >  2:528802/1‑148 (MQ=255)
                                                                                                                               tCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATa                                                 >  2:143424/1‑110 (MQ=255)
                                                                                                                               tCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATa                                                 <  1:143424/110‑1 (MQ=255)
                                                                                                                                 gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCATGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCgg                    >  2:777827/1‑137 (MQ=255)
                                                                                                                                 gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCATGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCgg                    <  1:777827/137‑1 (MQ=255)
                                                                                                                                 gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCtt       >  2:436095/1‑150 (MQ=255)
                                                                                                                                   tCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCCTGACCGCAAAACTGTCGGCGACGAAGGTCTTgt     <  2:428525/150‑2 (MQ=255)
                                                                                                                                   tCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGa     <  1:528802/150‑1 (MQ=255)
                                                                                                                                      tGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGt  >  2:354961/1‑150 (MQ=255)
                                                                                                                                          ggCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTc                      <  2:663230/126‑1 (MQ=255)
                                                                                                                                          ggCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTc                      >  1:663230/1‑126 (MQ=255)
                                                                                                                                             cTGTCGCTGGACCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGa     >  1:77716/1‑140 (MQ=255)
                                                                                                                                             cTGTCGCTGGACCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGa     <  2:77716/140‑1 (MQ=255)
                                                                                                                                             cTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCCCGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGa                                                   >  2:154356/1‑94 (MQ=255)
                                                                                                                                             cTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGa                                                   <  1:154356/94‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
GCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGT  >  NC_002947/4865924‑4866207

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