Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,001,864 G→A 9.4% R182Q (CGG→CAG)  PP_2625 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,001,8640GA 9.4% 130.7 / 2.0 51R182Q (CGG→CAG) PP_2625hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base A (2/3);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACA  >  NC_002947/3001717‑3002003
                                                                                                                                                   |                                                                                                                                           
gCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                           <  2:891617/150‑1 (MQ=255)
        gCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGccc                                                                                                                                   >  1:990373/1‑150 (MQ=255)
        gCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGccc                                                                                                                                   >  2:792982/1‑150 (MQ=255)
               cccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCg                                                                                                                            >  2:809746/1‑150 (MQ=255)
                ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGa                                                                                                                           <  1:312881/150‑1 (MQ=255)
                ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                             <  1:889231/148‑1 (MQ=255)
                ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                             >  2:889231/1‑148 (MQ=255)
                 cGGCTACATCTGATTGCCAGGGGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAg                                                                                                                          <  2:983066/150‑1 (MQ=255)
                  ggCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGCGCAGGCCCCGCCCGGCATCGAGc                                                                                                                         <  2:550687/150‑1 (MQ=255)
                   gCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCa                                                                                                                        <  2:714895/150‑1 (MQ=255)
                          cTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAgg                                                                                                                 >  1:116089/1‑150 (MQ=255)
                             attaacAGGCGTCCGATCGAGGCAAGCTGGGTCGAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGCGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGCCGGGGCCTGCGCGCCATCGAGCAGATCAGGAAg                                                                                                              >  2:448686/6‑150 (MQ=255)
                              ttGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGc                                                                                                             <  1:364134/150‑1 (MQ=255)
                                      ggTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGt                                                                                                     <  2:439186/149‑1 (MQ=255)
                                      cGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGTTCAGGAAGCTGATCGGt                                                                                                     <  2:651637/150‑1 (MQ=255)
                                        tCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGc                                                                                                   >  1:1108761/1‑150 (MQ=255)
                                                 ggCAAGCTGGGTCCAGGTACGCCTGAGCGGGGAGGGGTGGTCTACACGTGAGACGGCCGGGCAAGAGATGGTGGAGGGCCAAGTGCTGGGTGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCGGGAAGCTCATCGGTGCGCAGCGAAt                                                                                          >  2:51316/1‑150 (MQ=255)
                                                 ggCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTAt                                                                                          <  1:567568/150‑1 (MQ=255)
                                                  gCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGCCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGACCGGTGCCCAGCGTATc                                                                                         >  2:163845/1‑150 (MQ=255)
                                                  gCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATc                                                                                         <  2:389997/150‑1 (MQ=255)
                                                    aaGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGa                                                                                       >  1:385209/1‑150 (MQ=255)
                                                       cTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGATGt                                                                                    >  2:543513/1‑150 (MQ=255)
                                                            tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTTATCGGTGCCCAGCGTATCGAGGTTgtga                                                                               <  1:712076/150‑1 (MQ=255)
                                                            tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                               <  2:385209/150‑1 (MQ=255)
                                                               aGGTACGCCTGAGCGGGGAGGGCTGGTCTATCCGTGAGACGCCCTGGTAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgagtg                                                                            <  1:909739/150‑1 (MQ=255)
                                                                     gCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGAAGGAGAACGGCCAAGTGCCCGATGCCGGGCGGGGCATGCACGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTAGCAAGGTTGGGAGTGATCACg                                                                      >  2:266603/1‑150 (MQ=255)
                                                                       cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                           >  2:1005035/1‑79 (MQ=255)
                                                                       cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                           <  1:1005035/79‑1 (MQ=255)
                                                                       cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAGCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGc                                                                    >  1:10162/1‑150 (MQ=255)
                                                                         gAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGa                                                                  <  1:239952/150‑1 (MQ=255)
                                                                            cGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGt                                                               <  1:576265/150‑1 (MQ=255)
                                                                                   ggCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCGGAAAAGCACTTTCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACggg                                                        <  1:809746/150‑1 (MQ=255)
                                                                                   ggCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACggg                                                        <  1:792982/150‑1 (MQ=255)
                                                                                        gTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGAc                                                   >  1:427922/1‑150 (MQ=255)
                                                                                           tACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                           >  2:813627/1‑59 (MQ=255)
                                                                                           tACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                           <  1:813627/59‑1 (MQ=255)
                                                                                             cccGTGAGACGCCCTGGCAATTGAAGGAGCCCGGCCAAGTGCTCGATGCCGGCCCGGGCCTGGCCGCCATCGAGCAGATCAGGCAGCTGATCGGTGCCCAGCGCATCGGGCTTGTGAGCGCGCACGGCGAGGTCGACGGGCCGCCCTTCg                                              >  2:517406/1‑150 (MQ=255)
                                                                                             cccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCg                                              >  1:88414/1‑150 (MQ=255)
                                                                                               cGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTAATGGATCACGGCGAGGTCGACGGGCTGACCTTCGAt                                            >  2:261891/1‑150 (MQ=255)
                                                                                                   agaCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATgcgc                                        >  1:842515/1‑150 (MQ=255)
                                                                                                     aCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGAGGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCaa                                      <  2:88414/150‑1 (MQ=255)
                                                                                                           tGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCAGGCCCTTCCAGACAAATATCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                >  2:879564/1‑150 (MQ=255)
                                                                                                            ggCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTgg                               <  2:990373/150‑1 (MQ=255)
                                                                                                                   gATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCCGATCGTTGCCCAGCGTATCGAGGTTGTTAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCtt                        <  1:476687/150‑1 (MQ=255)
                                                                                                                       gagaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATc                    <  1:30967/150‑1 (MQ=255)
                                                                                                                        agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCa                   >  1:947991/1‑150 (MQ=255)
                                                                                                                        agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCa                   <  2:902720/150‑1 (MQ=255)
                                                                                                                               ccAAGTGCTCGATGTCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc            >  1:636213/1‑150 (MQ=255)
                                                                                                                                  agggCTCGAGGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTAGCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGcgcg         <  2:234230/147‑1 (MQ=255)
                                                                                                                                    tGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGaa       >  2:452190/1‑150 (MQ=255)
                                                                                                                                     gCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAg      <  2:842515/150‑1 (MQ=255)
                                                                                                                                      cTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGc     >  1:970621/1‑150 (MQ=255)
                                                                                                                                       tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca    >  1:622631/1‑150 (MQ=255)
                                                                                                                                       tCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGcc                <  2:671617/138‑1 (MQ=255)
                                                                                                                                       tCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGcc                >  1:671617/1‑138 (MQ=255)
                                                                                                                                         gATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca  >  1:761420/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACA  >  NC_002947/3001717‑3002003

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: ‑