Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,602 A→C 10.4% intergenic (‑257/‑87) pydP ← / → pydB NCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,547,6020AC10.4% 157.3 / 5.4 57intergenic (‑257/‑87)pydP/pydBNCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base C (3/3);  total (29/28)
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.36e-01

AGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCCC  >  NC_002947/4547461‑4547742
                                                                                                                                             |                                                                                                                                            
aGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACa                                                                                                                                        <  3:42893/148‑1 (MQ=255)
  tCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGa                                                                                                                                     <  6:245706/149‑1 (MQ=255)
  tCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGa                                                                                                                                     >  3:170899/1‑149 (MQ=255)
   cAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGAt                                                                                                                                    >  1:372990/1‑149 (MQ=255)
   cAAATACCAGGAAAACCCCGGCAACCATACGCCCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAAAGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTCTCAGGACTGACAGGaa                                                                                                                                    >  6:336660/1‑148 (MQ=255)
       tACCAGGAAAACCCCCGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTAAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGTCCTGACAAGAttcga                                                                                                                                <  8:139507/149‑4 (MQ=255)
         ccAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccg                                                                                                                              >  3:77952/1‑149 (MQ=255)
         ccAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcgag                                                                                                                               <  1:51275/148‑5 (MQ=255)
           aGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                                <  8:222522/145‑4 (MQ=255)
           aGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                                >  7:222522/1‑142 (MQ=255)
               aaaCCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgc                                                                                                                        >  3:202980/1‑149 (MQ=255)
                  ccccCGCGACCATGGGCGCGTTAAATTTTGCATTTTTATATCAACAATTTAAAAACCACAAATTTAATGGTAAAATTTTTTTGCTTAATCCAGAATCAGCATCTAGCCTCGATTATTGTCAGGACTGACAGTATTGAccgccattcccg                                                                                                                     <  8:311825/149‑1 (MQ=255)
                   aacGGCGCCCATTCTCGCCTTAATTATTAACTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTACTTGCTGAATCCAGAATCAGCATCTAGCATCGAATCTATTAAGTAATGACACGATTGAacgccctcaccgc                                                                                                                    >  4:75555/3‑149 (MQ=255)
                    ccGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgcccgcc                                                                                                                   <  7:51289/149‑1 (MQ=255)
                       gcgcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACa                                                                                                                                        >  2:174872/1‑125 (MQ=255)
                       gcgcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACa                                                                                                                                        <  1:174872/125‑1 (MQ=255)
                        cgcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgcccgccctgc                                                                                                               <  7:335763/149‑1 (MQ=255)
                        agcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgcccgccctg                                                                                                                >  4:52559/2‑148 (MQ=255)
                         gcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCa                                                                                                              >  6:103371/1‑149 (MQ=255)
                                 gcgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGa                                                                                                       >  2:68892/1‑148 (MQ=255)
                                 gcgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGa                                                                                                       >  8:13622/1‑148 (MQ=255)
                                     gatAATTTTTTACTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATg                                                                                                  <  8:375945/147‑1 (MQ=255)
                                          ttttttCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCa                                                                                                >  6:415425/1‑146 (MQ=255)
                                             tttCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAAc                                                                                          >  5:33412/1‑149 (MQ=255)
                                              ttCCTATTTATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAAc                                                                                          <  7:333822/148‑1 (MQ=255)
                                                ccTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCg                                                                                        <  1:68892/148‑1 (MQ=255)
                                                 cTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGc                                                                                      <  4:77952/149‑1 (MQ=255)
                                                    ttGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCa                                                                                    <  4:170899/148‑1 (MQ=255)
                                                             acaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                            >  1:403574/1‑147 (MQ=255)
                                                                  ttAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                                     >  7:283758/1‑149 (MQ=255)
                                                                    aaaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGc                                                                    <  3:52559/148‑1 (MQ=255)
                                                                     aaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCCAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCg                                                                   <  5:415425/148‑1 (MQ=255)
                                                                                tttAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCa                                                       <  8:25947/149‑1 (MQ=255)
                                                                                tttAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCa                                                       <  8:20064/149‑1 (MQ=255)
                                                                                    aTGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTc                                                   >  1:11646/1‑149 (MQ=255)
                                                                                             ttttCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGACCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGAt                                          >  4:149260/1‑149 (MQ=255)
                                                                                                     cTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgcg                                  >  4:208082/1‑149 (MQ=255)
                                                                                                     cTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGCACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgc                                   <  8:106387/148‑1 (MQ=255)
                                                                                                        aaTCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCa                               >  6:189902/1‑149 (MQ=255)
                                                                                                          tCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATTACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTGGATCCGTGGCGCCAcc                             >  2:301583/1‑149 (MQ=255)
                                                                                                          tCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTTTTGATCCGTGGCGCCAcc                             <  4:202980/149‑1 (MQ=255)
                                                                                                          tCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGACCGCCCTGCATCCATGACAATGCAAAGAACCGGCCTAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCAcc                             >  6:409144/1‑149 (MQ=255)
                                                                                                           ccAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTg                                     <  5:341032/140‑1 (MQ=255)
                                                                                                           ccAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTg                                     >  6:341032/1‑140 (MQ=255)
                                                                                                             aGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTg                          <  6:33412/149‑1 (MQ=255)
                                                                                                                 tCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTa                      >  6:174296/1‑149 (MQ=255)
                                                                                                                      aTCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCAc                 <  7:13622/149‑1 (MQ=255)
                                                                                                                           gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGa               <  3:208082/146‑1 (MQ=255)
                                                                                                                            ccTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagag           <  3:75555/149‑1 (MQ=255)
                                                                                                                             cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagaga          <  2:372990/149‑1 (MQ=255)
                                                                                                                             cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagaga          <  8:283758/149‑1 (MQ=255)
                                                                                                                              gtcaaTCCTGTCAGTCCTGACAAGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagagag         >  6:240097/4‑149 (MQ=255)
                                                                                                                              gtcaaTCCTGTCAGTCCTGACAAGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagagag         >  8:345469/4‑149 (MQ=255)
                                                                                                                                  aTCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTAc     <  2:11646/149‑1 (MQ=255)
                                                                                                                                   tCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTAcc    >  7:317712/1‑149 (MQ=255)
                                                                                                                                   tCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTAcc    >  1:40325/1‑149 (MQ=255)
                                                                                                                                     ttGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTAcccc  >  3:305702/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                            
AGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCCC  >  NC_002947/4547461‑4547742

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