Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,128,188 C→G 9.3% A171A (GCC→GCG htrG → signal transduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,128,1880CG 9.3% 159.8 / 7.6 54A171A (GCC→GCGhtrGsignal transduction protein
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base G (2/3);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

CATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCGGCAAG  >  NC_002947/4128041‑4128336
                                                                                                                                                   |                                                                                                                                                    
cATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCg                                                                                                                                                    >  2:1342260/1‑150 (MQ=255)
          aGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGAt                                                                                                                                          <  1:1074275/150‑1 (MQ=255)
            cTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCGGGCCTgggtatcgc                                                                                                                                        <  2:1250689/150‑7 (MQ=255)
              ggAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATgaga                                                                                                                                      >  2:1070931/1‑150 (MQ=255)
                     ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCa                                                                                                                               >  1:483497/1‑150 (MQ=255)
                          tACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTGGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGAGCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCGGGGGGGTGAGAACAAGCAGGtgg                                                                                                                          >  2:1235708/1‑149 (MQ=255)
                              gggCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAg                                                                                                                              >  1:103964/1‑142 (MQ=255)
                              gggCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAg                                                                                                                              <  2:103964/142‑1 (MQ=255)
                                              cTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTAt                                                                                                      <  1:151758/150‑1 (MQ=255)
                                                  aTTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACCAGTAGCTTATACAGAGCGATTCGGACCTGCGCGATACCCCGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcg                                                                                                  <  1:324185/150‑1 (MQ=255)
                                                   ttCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcggc                                                                                                >  2:53247/1‑150 (MQ=255)
                                                      gcgcAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgcag                                                                                              >  1:928707/1‑150 (MQ=255)
                                                               gcTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCgg                                                                                     >  2:978001/1‑150 (MQ=255)
                                                                   attgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCa                                                                                 <  2:1090628/150‑1 (MQ=255)
                                                                      gatgaGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCg                                                                                      <  2:1282691/142‑1 (MQ=255)
                                                                      gatgaGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCg                                                                                      >  1:1282691/1‑142 (MQ=255)
                                                                           gTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTgatgat                                                                                                                      >  2:1254147/1‑105 (MQ=255)
                                                                           gTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTgatgat                                                                                                                      <  1:1254147/105‑1 (MQ=255)
                                                                               gAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCg                                                                     <  1:289210/150‑1 (MQ=255)
                                                                                    aCGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCt                                                                >  1:275964/1‑150 (MQ=255)
                                                                                    aCGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCt                                                                >  1:1036039/1‑150 (MQ=255)
                                                                                      gCAACAAGGCACTGACAGAGCAGCTCGACCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACAGGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGa                                                              <  1:1235708/150‑1 (MQ=255)
                                                                                        aataaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATc                                                            <  1:20944/147‑1 (MQ=255)
                                                                                           aaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTg                                                         >  1:953826/1‑150 (MQ=255)
                                                                                                tcTGATCGAGCTGCTCGTTCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGc                                                    >  2:621820/2‑150 (MQ=255)
                                                                                                         gcagcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                                                                        >  1:716329/1‑57 (MQ=255)
                                                                                                         gcagcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                                                                        <  2:716329/57‑1 (MQ=255)
                                                                                                            gcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTgg                                        <  2:1011380/150‑1 (MQ=255)
                                                                                                              tCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCGGACGGGGCGGATCCTGTCGGCAATGAACCGCgggg                                      >  2:944854/1‑149 (MQ=255)
                                                                                                               cGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGgcgc                                    <  1:53247/150‑1 (MQ=255)
                                                                                                                     gagCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCaagaa                               >  2:743175/1‑150 (MQ=255)
                                                                                                                     gagCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCaagaa                               >  2:751194/1‑150 (MQ=255)
                                                                                                                      agCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCaagaag                              >  1:837335/1‑150 (MQ=255)
                                                                                                                      agCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGGGGGCGCaagaag                              >  1:487669/1‑150 (MQ=255)
                                                                                                                        ccAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAgaagaa                            <  1:1147514/150‑1 (MQ=255)
                                                                                                                        ccAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAgaagaa                            >  1:825604/1‑150 (MQ=255)
                                                                                                                        ccAGTCGGACCTGCGCGATACCAACGGGTTGCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAgaagaa                            <  1:1070931/150‑1 (MQ=255)
                                                                                                                                aCCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGcc                                                                                       <  2:1184984/83‑1 (MQ=255)
                                                                                                                                aCCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGcc                                                                                       >  1:1184984/1‑83 (MQ=255)
                                                                                                                                aCCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTACGGCGGCAGCCTTGGCCGTGCAGGCTTGCTGGCGGGGCTGGCCCTGCCGGCACTGAGCCGGGTGCGCAATAAGAACGATCGCt                    <  2:873964/150‑1 (MQ=255)
                                                                                                                                   tGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgt                                                                                                          >  1:288543/1‑61 (MQ=255)
                                                                                                                                   tGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgt                                                                                                          <  2:288543/61‑1 (MQ=255)
                                                                                                                                   tGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGt                 >  2:246255/1‑150 (MQ=255)
                                                                                                                                          tACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGc          >  2:1297275/1‑150 (MQ=255)
                                                                                                                                          cccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGc                                                                 >  2:469602/3‑95 (MQ=255)
                                                                                                                                          cccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGc                                                                 <  1:469602/93‑1 (MQ=255)
                                                                                                                                           ccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCt         >  2:217571/2‑150 (MQ=255)
                                                                                                                                             ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCg       >  1:802465/1‑150 (MQ=255)
                                                                                                                                             ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCg       <  1:1342260/150‑1 (MQ=255)
                                                                                                                                             ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTAGTTc               >  1:118585/1‑142 (MQ=255)
                                                                                                                                             ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTAGTTc               <  2:118585/142‑1 (MQ=255)
                                                                                                                                                ggCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTg                                         >  2:238197/1‑113 (MQ=255)
                                                                                                                                                ggCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTg                                         <  1:238197/113‑1 (MQ=255)
                                                                                                                                                 gCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCGGCaa   <  2:1036039/150‑1 (MQ=255)
                                                                                                                                                  cccGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCGGCAAg  <  1:621820/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCGGCAAG  >  NC_002947/4128041‑4128336

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