Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,609 G→A 11.0% intergenic (‑264/‑80) 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,6090GA11.0% 142.4 / 6.4 55intergenic (‑264/‑80)pydP/pydBNCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base A (3/3);  total (29/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 = 8.91e-01

TCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCCCGCCGATGTCCTG  >  NC_002947/4547463‑4547754
                                                                                                                                                  |                                                                                                                                                 
tCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGa                                                                                                                                                 >  3:170899/1‑149 (MQ=255)
tCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGa                                                                                                                                                 <  6:245706/149‑1 (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)
                      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                                                 >  6:341032/1‑140 (MQ=255)
                                                                                                         ccAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTg                                                 <  5:341032/140‑1 (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)
                                                                                                                                               aCAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCCAGGAAAACCGCATGTCCCTGTTGATACGTGGCGCCACCGAGGTTACCAACGAAGACAGTTACCCCGCCGATGTCCtg  >  3:30953/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCCCGCCGATGTCCTG  >  NC_002947/4547463‑4547754

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