Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,463,556 A→G 8.3% W268R (TGG→CGG)  PP_3078 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,463,5560AG 8.3% 133.4 / 2.1 48W268R (TGG→CGG) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base G (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCAGGCT  >  NC_002947/3463419‑3463703
                                                                                                                                         |                                                                                                                                                   
cGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATg                                                                                                                                                >  8:38057/1‑143 (MQ=255)
cGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATg                                                                                                                                                <  7:38057/143‑1 (MQ=255)
     cGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGt                                                                                                                                     <  4:120310/149‑1 (MQ=255)
            cTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAgg                                                                                                                              >  3:130109/1‑149 (MQ=255)
            cTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAgg                                                                                                                              >  3:110460/1‑149 (MQ=255)
            cTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAgg                                                                                                                              >  2:173822/1‑149 (MQ=255)
               gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCGACCATGCCGGGTcgc                                                                                                                                       <  5:64541/137‑1 (MQ=255)
               gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGa                                                                                                                           <  4:132083/149‑1 (MQ=255)
               gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGa                                                                                                                           >  1:251528/1‑149 (MQ=255)
               gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTATGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTcgc                                                                                                                                       >  6:64541/1‑137 (MQ=255)
               gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCCCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTg                                                                                                                            <  3:99786/148‑1 (MQ=255)
                 ggTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAtg                                                                                                                         <  1:75422/149‑1 (MQ=255)
                            tAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                              >  2:153406/1‑149 (MQ=255)
                            tAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                              >  2:148557/1‑149 (MQ=255)
                               ccACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTc                                                                                                           <  4:110460/149‑1 (MQ=255)
                                 aCCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCACGTTCGAAGGTGATGCTGCGGCCGTTGTCGa                                                                                                         >  6:199615/1‑149 (MQ=255)
                                  ccGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCACGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                                        <  5:199615/149‑1 (MQ=255)
                                    gCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATgc                                                                                                      <  8:240587/149‑1 (MQ=255)
                                      cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATgcgc                                                                                                    <  5:13732/149‑1 (MQ=255)
                                       tttGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGcc                                                                                                   >  6:184195/1‑149 (MQ=255)
                                               ccTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                              <  7:92910/130‑1 (MQ=255)
                                               ccTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                              >  8:92910/1‑130 (MQ=255)
                                                 tGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCAGGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTa                                                                                         >  1:194691/1‑149 (MQ=255)
                                                      cgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGcc                                                                                    <  1:153406/149‑1 (MQ=255)
                                                      cgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGcc                                                                                    <  1:148557/149‑1 (MQ=255)
                                                       gcCACTTCGGTGTCGCCGAAGTACTCGATGCTCAGTTTGTCGAAATTGTAGCGGCCACGGTTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGccc                                                                                   <  7:183913/149‑1 (MQ=255)
                                                        ccACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                                    >  8:188734/1‑99 (MQ=255)
                                                        ccACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                                    <  7:188734/99‑1 (MQ=255)
                                                                     gCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGc                                                                      <  2:251528/148‑1 (MQ=255)
                                                                            tACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGaa                                                              <  1:173822/149‑1 (MQ=255)
                                                                             aCTCGATGCGCAGTTTGTCGAAATTGTAGCGGCGACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGaa                                                              <  4:130109/148‑1 (MQ=255)
                                                                                     cgcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCAcc                                                     >  3:123130/1‑149 (MQ=255)
                                                                                      gcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                                    >  5:176890/1‑133 (MQ=255)
                                                                                      gcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                                    <  6:176890/133‑1 (MQ=255)
                                                                                                gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGt                                          >  7:138211/1‑149 (MQ=255)
                                                                                                              cACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                                 >  6:200484/1‑128 (MQ=255)
                                                                                                              cACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                                 <  5:200484/128‑1 (MQ=255)
                                                                                                                cGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTc                          <  8:138211/149‑1 (MQ=255)
                                                                                                                 ggCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCg                         >  5:152781/1‑149 (MQ=255)
                                                                                                                 ggCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCg                         >  6:217018/1‑149 (MQ=255)
                                                                                                                   atggtggGCACGTGACCTGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTg                       >  3:122771/2‑149 (MQ=255)
                                                                                                                    tGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGc                      >  7:180975/1‑149 (MQ=255)
                                                                                                                      gCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTc                    >  7:128024/1‑149 (MQ=255)
                                                                                                                       tgggCACGTGACCTGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCg                   >  6:54135/2‑149 (MQ=255)
                                                                                                                          gCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGc                <  4:123130/149‑1 (MQ=255)
                                                                                                                              gTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTa            >  2:307174/1‑149 (MQ=255)
                                                                                                                               tgacctGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTAc           <  5:54135/144‑1 (MQ=255)
                                                                                                                                      cccGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCAg     <  4:122771/148‑1 (MQ=255)
                                                                                                                                        caccaTGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCAGGCt  >  6:140880/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
CGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCAGGCT  >  NC_002947/3463419‑3463703

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