Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,046,618 C→A 7.5% P102H (CCT→CAT)  pstB‑I → phosphate ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,046,6180CA 7.5% 166.2 / 3.6 53P102H (CCT→CAT) pstB‑Iphosphate ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base A (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

AGATCACCGGTTTCATCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGCAAGGCG  >  NC_002947/3046471‑3046764
                                                                                                                                                   |                                                                                                                                                  
agaTCACCGGTTTCATCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCt                                                                                                                                                   <  5:25633/149‑1 (MQ=255)
    cACCGGTTTCATCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCtgttg                                                                                                                                               >  3:264903/1‑149 (MQ=255)
     aCCGGTTTCATCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCtgttg                                                                                                                                               >  1:359180/1‑148 (MQ=255)
          tttCATCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTTAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTgc                                                                                                                                         <  2:80413/149‑1 (MQ=255)
              aTCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGcgccg                                                                                                                                     >  8:236602/1‑149 (MQ=255)
              aTCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGcgcc                                                                                                                                      <  8:46111/148‑1 (MQ=255)
                   cccTTCGGGCGGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCtata                                                                                                                                <  4:108706/149‑1 (MQ=255)
                    ccTTCGGGCTGCTGCAAGAGTACCGTGCTGCGTAGCTTGAACAGTATGAATGACCTGTTGAATGGCTTCAGCTTCGAAGGCAATGTGAACTTTCTTGGGCAGGATGTCTACTGCAAGGGTGTCTACCCTGTTGTAGTGCGCAGCTatat                                                                                                                               <  4:253052/149‑1 (MQ=255)
                    ccTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTatat                                                                                                                               >  4:152285/1‑149 (MQ=255)
                    ccTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTatat                                                                                                                               >  6:229077/1‑149 (MQ=255)
                     cTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATc                                                                                                                              >  3:248790/1‑149 (MQ=255)
                          ggCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGTGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCAt                                                                                                                         <  7:15949/149‑1 (MQ=255)
                            cTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATgg                                                                                                                       <  2:371521/149‑1 (MQ=255)
                                             tgctgcGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCtg                                                                                                                                                  >  3:121296/1‑105 (MQ=255)
                                             tgctgcGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCtg                                                                                                                                                  <  4:121296/105‑1 (MQ=255)
                                             tgctgcGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAAc                                                                                                      >  8:233182/1‑149 (MQ=255)
                                                      gCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTtctc                                                                                              <  5:229077/148‑1 (MQ=255)
                                                      gCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTtctc                                                                                              <  6:150791/148‑1 (MQ=255)
                                                           aaCCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGa                                                                                         <  2:359180/148‑1 (MQ=255)
                                                                tATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCt                                                                                   <  8:354555/149‑1 (MQ=255)
                                                                          cTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAAtgt                                                                         >  7:207457/1‑149 (MQ=255)
                                                                          cTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAAtgt                                                                         >  8:346088/1‑149 (MQ=255)
                                                                          cTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAAtgt                                                                         >  4:89322/1‑149 (MQ=255)
                                                                           tGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTtc                                                                                                <  4:220807/125‑1 (MQ=255)
                                                                           tGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTtc                                                                                                >  3:220807/1‑125 (MQ=255)
                                                                               gAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAAc                                                                                                      >  4:222662/1‑115 (MQ=255)
                                                                               gAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAAc                                                                                                      <  3:222662/115‑1 (MQ=255)
                                                                               gAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCt                                                                    >  5:69404/1‑149 (MQ=255)
                                                                                aaGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCATGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCtt                                                                   >  1:44678/1‑149 (MQ=255)
                                                                                aaGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCATGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCtt                                                                   >  1:50428/1‑149 (MQ=255)
                                                                                    gCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGc                                                               <  4:248790/149‑1 (MQ=255)
                                                                                     cTTCCGCTTCGAAGGCCAGGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCAAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGcc                                                              >  7:282534/1‑149 (MQ=255)
                                                                                                aaGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTTGGCCTGCGCCTGAAc                                                   >  5:360363/1‑149 (MQ=255)
                                                                                                aaGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAAc                                                   <  7:233182/149‑1 (MQ=255)
                                                                                                         tgCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACa                                            >  1:9313/1‑147 (MQ=255)
                                                                                                         tgCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACa                                            >  7:348454/1‑147 (MQ=255)
                                                                                                           cACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAgg                                        >  1:70349/1‑149 (MQ=255)
                                                                                                            aCTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGc                                       >  7:72794/1‑149 (MQ=255)
                                                                                                            aCTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCATGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGc                                       <  2:50428/149‑1 (MQ=255)
                                                                                                            aCTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCATGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGc                                       <  2:44678/149‑1 (MQ=255)
                                                                                                              tttCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGa                                     >  8:4560/1‑149 (MQ=255)
                                                                                                                tCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGa                                     <  3:152285/147‑1 (MQ=255)
                                                                                                                    gggCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTgg                                <  8:72794/148‑1 (MQ=255)
                                                                                                                    gggCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTgg                                <  7:236602/148‑1 (MQ=255)
                                                                                                                     ggCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTTGGCCTGCGCCTGAACCGCTACAAAGGCGACTTggg                               <  6:360363/148‑1 (MQ=255)
                                                                                                                     ggCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTggg                               <  3:89322/148‑1 (MQ=255)
                                                                                                                              tCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGt                      <  8:207457/148‑1 (MQ=255)
                                                                                                                                    gCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCAc               >  6:257797/1‑149 (MQ=255)
                                                                                                                                           tgtCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGc        <  7:346088/149‑1 (MQ=255)
                                                                                                                                            gtCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGc        >  3:79588/1‑148 (MQ=255)
                                                                                                                                             tCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGc        >  4:48285/1‑147 (MQ=255)
                                                                                                                                             tCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCATGAAGCACGCCTTGCaa      <  2:70349/149‑1 (MQ=255)
                                                                                                                                               gACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGCAAgg    >  5:105165/1‑149 (MQ=255)
                                                                                                                                                 cccTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGCAAGgcg  >  4:298542/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGATCACCGGTTTCATCGGCCCTTCGGGCTGCGGCAAGAGTACCGTGCTGCGTAGCTTGAACCGTATGAATGACCTGGTGAAGGGCTTCCGCTTCGAAGGCCATGTGCACTTTCTTGGGCAGGATGTCTACGGCAAGGGTGTCGACCCTGTTGTAGTGCGCCGCTATATCGGCATGGTGTTCCAGCAGCCCAACCCGTTCTCGATGAGCATCTTCGACAATGTGGCCTTCGGCCTGCGCCTGAACCGCTACAAAGGCGACTTGGGCGATCGCGTGAAGCACGCCTTGCAAGGCG  >  NC_002947/3046471‑3046764

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