Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,909,517 C→A 8.7% G446V (GGC→GTC)  PP_1710 ← phthalate permease family MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,909,5170CA 8.7% 130.7 / 2.4 46G446V (GGC→GTC) PP_1710phthalate permease family MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base A (2/2);  total (21/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 = 2.61e-01

TTGGACTATCCAAGGCACCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGCCGATCACGATCGGTGTGACG  >  NC_002947/1909368‑1909666
                                                                                                                                                     |                                                                                                                                                      
ttggACTATCCAAGGCACCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGc                                                                                                                                                        >  1:1194496/1‑150 (MQ=255)
     cTATCCAAGGCACCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCa                                                                                                                                                   <  2:888788/150‑1 (MQ=255)
             ggCACCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAg                                                                                                                                           >  2:1073309/1‑150 (MQ=255)
                aCCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTc                                                                                                                                        >  2:14425/1‑150 (MQ=255)
                aCCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTc                                                                                                                                        >  2:1226438/1‑150 (MQ=255)
                  cTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGc                                                                                                                                                        >  2:414465/1‑132 (MQ=255)
                  cTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGc                                                                                                                                                        <  1:414465/132‑1 (MQ=255)
                  cTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTc                                                                                                                                             <  2:418051/143‑1 (MQ=255)
                  cTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTc                                                                                                                                             >  1:418051/1‑143 (MQ=255)
                   tgcaCCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGAt                                                                                                                                     <  1:810105/147‑1 (MQ=255)
                     tACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCACAGTCTTCAGCTCGATGc                                                                                                                                   <  2:228342/150‑1 (MQ=255)
                       ccGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGa                                                                                                                                 <  1:1177890/150‑1 (MQ=255)
                       ccGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGa                                                                                                                                 <  2:260248/150‑1 (MQ=255)
                              tttCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCACACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑Ct                                                                                                                         <  1:976716/150‑1 (MQ=255)
                                tCGCGGCTGACGCTGCTCCTACCGGGATAGCGCCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTc                                                                                                                       <  2:88012/150‑1 (MQ=255)
                                                 ccTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAg                                                                                                                            <  2:913200/129‑1 (MQ=255)
                                                 ccTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAg                                                                                                                            >  1:913200/1‑129 (MQ=255)
                                                  cTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTaa                                                                                                     <  1:14425/150‑1 (MQ=255)
                                                      aGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACa                                                                                                 >  2:201067/1‑150 (MQ=255)
                                                        ggATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAgt                                                                                               <  2:886189/150‑1 (MQ=255)
                                                        ggATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGACCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAgt                                                                                               >  1:905225/1‑150 (MQ=255)
                                                        ggATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGACCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAgt                                                                                               >  1:905452/1‑150 (MQ=255)
                                                           tAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTAGCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CCTCCCAAGCAGCAGCGAGTAAGACAGtgtg                                                                                            >  2:376920/1‑150 (MQ=255)
                                                            aGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGACCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGc                                                                                           <  2:905225/150‑1 (MQ=255)
                                                            aGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGACCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGc                                                                                           <  2:905452/150‑1 (MQ=255)
                                                                tCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCGGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGACCAAAACCGGGTGCCCaa                                                                                       >  2:568944/1‑149 (MQ=255)
                                                                tCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCCGCGAGTAAGACAGTGTGCCCAt                                                                                       >  2:22811/1‑150 (MQ=255)
                                                                     cccgttggCAATACCGTCCCCCTGCAGGTACGGCTCTAGCCGCGAAGAGGCCGGTACCGTCAATACTCGCCTTGTCTCAGCCCGTCACCGGCAAGTTGGTGCG‑TACAGCCTGCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACg                                                                                  <  2:685427/144‑1 (MQ=255)
                                                                      ccAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCaa                                                                                                                   <  1:530782/116‑1 (MQ=255)
                                                                      ccAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCaa                                                                                                                   >  2:530782/1‑116 (MQ=255)
                                                                       cAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGcc                                                                                <  2:122373/150‑1 (MQ=255)
                                                                        aGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCca                                                                               <  1:730324/150‑1 (MQ=255)
                                                                          gggCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAagcagc                                                                                                              >  2:191563/1‑117 (MQ=255)
                                                                          gggCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAagcagc                                                                                                              <  1:191563/117‑1 (MQ=255)
                                                                            gCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGCCTTCAGCTCGACGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCAcccc                                                                           >  2:535195/1‑150 (MQ=255)
                                                                                        cccTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCa                                                               >  2:953287/1‑150 (MQ=255)
                                                                                                       tttCGCCGCAAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGaa                                                <  2:230800/150‑1 (MQ=255)
                                                                                                          aGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGa                                             <  1:953287/150‑1 (MQ=255)
                                                                                                                     ggCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGt                                  <  2:806996/150‑1 (MQ=255)
                                                                                                                      gCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTg                                 >  2:1221650/1‑150 (MQ=255)
                                                                                                                        cGGTACCGTCAATACACGCCTTGGCTCAGCGCGCAGTCGTCAGCTCGATGCGATACAG‑CCTGCGAATGCGTTGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCa                               <  1:22811/150‑1 (MQ=255)
                                                                                                                        cGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCa                               >  1:1081192/1‑150 (MQ=255)
                                                                                                                               gTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCGCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACTACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGgc                        >  1:734637/1‑150 (MQ=255)
                                                                                                                                  aaTACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGcc                     >  2:442207/1‑150 (MQ=255)
                                                                                                                                    tACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGCCGa                   <  2:1081192/150‑1 (MQ=255)
                                                                                                                                    tACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGCCGa                   <  1:1073309/150‑1 (MQ=255)
                                                                                                                                     acacGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGCCGAt                  >  1:922787/1‑150 (MQ=255)
                                                                                                                                        cGCCTTGGCTCAGCCCGCAGGCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCCACGCCACACCACCGATGAACCAGAACGCCATGGCGAACCACCCCCTGGCGTGCACGATGCCGCCGAGCAc               >  2:491543/1‑150 (MQ=255)
                                                                                                                                                   aGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGCCGATCACGATCGGTGTGa    <  2:922787/150‑1 (MQ=255)
                                                                                                                                                     cccGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACAACAACGCCATGGCGAACGAACCACTGGCGTGCACCCTTGCCCCCACCACGATCGCTGCGACg  >  2:582780/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTGGACTATCCAAGGCACCTGTACCGGCCCTTTCGCGGCTGAAGCCGCTCCTACAGGGATAGCGTCAAACCCAGGGGCAATACCGTCCCCCTGCAGGCACGGCTTTAGCCGCGAAGAGGCCGGTACCGTCAATACACGCCTTGGCTCAGCCCGCAGTCTTCAGCTCGATGCGATACAG‑CTTCCCAAGCAGCAGCGAGTAAGACAGTGTGCCCATCAACGCCACACCACCGATGAACCAGAACGCCATGGCGAACGAGCCACTGGCGTGCACGATGGCGCCGATCACGATCGGTGTGACG  >  NC_002947/1909368‑1909666

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