Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,914,477 G→C 9.7% R118R (CGC→CGG)  PP_3454 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,914,4770GC 9.7% 153.7 / 8.6 53R118R (CGC→CGG) PP_3454transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (17/30);  new base C (3/2);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

GCGCCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGCG  >  NC_002947/3914329‑3914619
                                                                                                                                                    |                                                                                                                                              
tcgcCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTgcggc                                                                                                                                               >  1:1137024/2‑150 (MQ=255)
gcgcCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTgcggc                                                                                                                                               >  1:46399/1‑150 (MQ=255)
      gCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTgcgg                                                                                                                                                <  2:887906/143‑1 (MQ=255)
      gCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTgcgg                                                                                                                                                >  1:887906/1‑143 (MQ=255)
      gCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCgcagca                                                                                                                                         <  2:913462/150‑1 (MQ=255)
           accaGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTg                                                                                                                                    <  1:506029/150‑1 (MQ=255)
               gcagGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGGCGTTGTCCAGGCGCGGTTCGTTGCGACTGGTGCGACGCAGCAGCGCGCGCa                                                                                                                                >  2:312528/1‑150 (MQ=255)
                   gTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCg                                                                                                                            <  1:849056/150‑1 (MQ=255)
                   gTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCg                                                                                                                            >  2:404398/1‑150 (MQ=255)
                         cTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCTTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCcag                                                                                                                      <  1:404398/150‑1 (MQ=255)
                         cTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCcag                                                                                                                      >  1:467584/1‑150 (MQ=255)
                                 tGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCATGTCGTGGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCAcgc                                                                                                              <  1:782678/150‑1 (MQ=255)
                                   gACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGt                                                                                                            <  2:737714/150‑1 (MQ=255)
                                    aCAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTg                                                                                                           >  2:749181/1‑150 (MQ=255)
                                     cAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTgg                                                                                                          <  2:467584/150‑1 (MQ=255)
                                     cAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTgg                                                                                                          <  2:817568/150‑1 (MQ=255)
                                        cTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGCCGcagcagt                                                                                                                                       >  2:1106422/1‑117 (MQ=255)
                                        cTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGCCGcagcagt                                                                                                                                       <  1:1106422/118‑2 (MQ=255)
                                             cctcctCGCCGCGCCAAATCACCGTGCGCTCGGAAAGGTCGATCCGCCAGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAgg                                                                                                  >  2:849740/1‑150 (MQ=255)
                                                  ccgccggGCCAACCCACCCGTCTCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTg                                                                                             <  2:973692/149‑1 (MQ=255)
                                                   cgccgcGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCTGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTgg                                                                                            <  2:594246/150‑1 (MQ=255)
                                                       gcgcCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTACCAGGGTTTGGTCAc                                                                                        <  2:46399/150‑1 (MQ=255)
                                                                  ccGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCGAGCAGCACGCGTGGCTCACAGGGTGTGGTCACGTAGTCGTCgg                                                                             >  2:263356/1‑150 (MQ=255)
                                                                            ggCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTc                                                                   <  1:159667/150‑1 (MQ=255)
                                                                               cAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAg                                                                >  2:1205324/1‑150 (MQ=255)
                                                                               cAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGAGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAg                                                                <  1:343721/150‑1 (MQ=255)
                                                                               cAGGTCGATACGCAGGCCACCAACCAGCATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAg                                                                <  1:187025/150‑1 (MQ=255)
                                                                                aGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGATGAGGCGCCGGCACGGGCGCACCCGGTCTAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAgg                                                               <  1:749181/150‑1 (MQ=255)
                                                                                          gCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCt                                                     <  2:1106707/150‑1 (MQ=255)
                                                                                                  ccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCCCGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATgg                                             <  2:1137024/150‑1 (MQ=255)
                                                                                                  ccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCAc                                                                                        <  1:224720/107‑1 (MQ=255)
                                                                                                  ccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCAc                                                                                        >  2:224720/1‑107 (MQ=255)
                                                                                                  ccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATgg                                             >  1:224043/1‑150 (MQ=255)
                                                                                                  ccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACGGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATgg                                             <  1:312528/150‑1 (MQ=255)
                                                                                                      tcaggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCAcgc                                                                                                              <  1:621462/81‑1 (MQ=255)
                                                                                                      tcaggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCAcgc                                                                                                              >  2:621462/1‑81 (MQ=255)
                                                                                                      tcaggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTAt                                         >  1:108972/1‑150 (MQ=255)
                                                                                                      tcaggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCCGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTAt                                         >  2:1014115/1‑150 (MQ=255)
                                                                                                               ggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGGGCCTGGTCATGGCTATCGTCGCggg                                <  1:766352/150‑1 (MQ=255)
                                                                                                                    ttGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCCCCCGGGCCAGCAGCACGCGTGGCTCACAgg                                                                                                  >  2:838353/1‑79 (MQ=255)
                                                                                                                    ttGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAgg                                                                                                  <  1:838353/79‑1 (MQ=255)
                                                                                                                     tGTCCAGGTGCGGTTCGTTGACATTGCTGCGGCGCAGCAGCGTTCGCACTCGGGCCTGCAGCACGCGTGGCTCACAGGGTTTGGTTACGTAGTTGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCa                          <  2:1135203/150‑1 (MQ=255)
                                                                                                                      gTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTg                                                                                             >  2:1169301/1‑82 (MQ=255)
                                                                                                                      gTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTg                                                                                             <  1:1169301/82‑1 (MQ=255)
                                                                                                                                 gTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTg              >  2:814529/1‑150 (MQ=255)
                                                                                                                                      ttGACACTGCTGCGGCGCAGCAGCGGGCGCACCCGGGCCAGCAGCACGGGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAgg         <  1:264174/150‑1 (MQ=255)
                                                                                                                                        cacgctgctgCGCCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAgg         <  2:878244/144‑1 (MQ=255)
                                                                                                                                        cacgctgctgCGCCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAgg         >  1:878244/5‑148 (MQ=255)
                                                                                                                                         acgctgctgCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGc                                                              >  2:1186600/4‑94 (MQ=255)
                                                                                                                                         acgctgctgCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGc                                                              <  1:1186600/91‑1 (MQ=255)
                                                                                                                                           acTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCATGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTg    <  1:329573/150‑1 (MQ=255)
                                                                                                                                            ctgcttCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGc   <  1:849740/150‑1 (MQ=255)
                                                                                                                                             tgctgcGCCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGCg  >  2:627304/1‑150 (MQ=255)
                                                                                                                                                tgcGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTcgc                                   <  1:736877/114‑1 (MQ=255)
                                                                                                                                                tgcGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTcgc                                   >  2:736991/1‑114 (MQ=255)
                                                                                                                                                tgcGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTcgc                                   >  2:736877/1‑114 (MQ=255)
                                                                                                                                                tgcGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTcgc                                   <  1:736991/114‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCGCCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGCG  >  NC_002947/3914329‑3914619

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