Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,852,078 T→G T467P (ACG→CCG)  yidK ← solute:sodium symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,852,0780TG96.5% 285.1 / 1.2 87T467P (ACG→CCG) yidKsolute:sodium symporter family transporter
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (32/52);  total (33/54)
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.66e-01

CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851991‑3852162
                                                                                       |                                                                                    
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    >  2:977971/1‑90 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:99779/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:859643/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:807425/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:68986/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:1423615/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:556374/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    >  1:528542/1‑90 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:142501/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:146979/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:2450/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:234461/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    >  1:167988/1‑90 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:425706/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:1186377/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    >  1:1033173/1‑90 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:1385278/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:89529/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:719393/90‑1 (MQ=255)
  gCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAg                                                                                  <  2:248943/90‑1 (MQ=255)
  gCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAg                                                                                  <  2:5076/90‑1 (MQ=255)
  gCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAg                                                                                  <  2:1261919/90‑1 (MQ=255)
     aaCGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGcc                                                                               <  2:42739/90‑1 (MQ=255)
       cGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAg                                                                             <  2:125166/90‑1 (MQ=255)
       cGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAg                                                                             <  2:670250/90‑1 (MQ=255)
        gCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGc                                                                            <  1:1521341/90‑1 (MQ=255)
                   ggTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGaa                                                                 >  2:213448/1‑90 (MQ=255)
                     tGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAAt                                                               >  2:660964/1‑90 (MQ=255)
                     tGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAAt                                                               >  1:1392957/1‑90 (MQ=255)
                      gAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    >  1:539588/1‑68 (MQ=255)
                      gAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:539588/68‑1 (MQ=255)
                      gAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATg                                                              <  2:136230/90‑1 (MQ=255)
                      gAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATg                                                              >  2:1311398/1‑90 (MQ=255)
                           gCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGt                                                         <  1:922054/90‑1 (MQ=255)
                            cGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTc                                                        <  1:586359/90‑1 (MQ=255)
                             gTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGt                                                         >  2:1100375/1‑88 (MQ=255)
                             gTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGt                                                         <  1:1100375/88‑1 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    >  1:685605/1‑90 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    <  2:1137185/90‑1 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    <  2:424621/90‑1 (MQ=255)
                                   cgcgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTc                                                 <  2:396636/90‑1 (MQ=255)
                                   cgcgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTc                                                 <  2:117175/90‑1 (MQ=255)
                                     cgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCAc                                               >  1:213004/1‑90 (MQ=255)
                                     cgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCAc                                               >  1:1024679/1‑90 (MQ=255)
                                         ttGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGa                                           <  1:660964/90‑1 (MQ=255)
                                          tGATAAAACCGATCCCCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  2:1048661/90‑1 (MQ=255)
                                          tGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  2:371577/90‑1 (MQ=255)
                                          tGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACAGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  2:1325629/90‑1 (MQ=255)
                                               aaaCCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTg                                     >  2:782758/1‑90 (MQ=255)
                                                 aCCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGAt                                   >  2:1063412/1‑90 (MQ=255)
                                                 aCCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGAt                                   <  2:1238573/90‑1 (MQ=255)
                                                 aCCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGAt                                   >  1:1346612/1‑90 (MQ=255)
                                                  ccGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATg                                  <  1:62834/90‑1 (MQ=255)
                                                  ccGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATg                                  <  1:133773/90‑1 (MQ=255)
                                                     aTCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTg                               >  1:1356148/1‑90 (MQ=255)
                                                     aTCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTg                               >  1:739352/1‑90 (MQ=255)
                                                       cACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGAt                             <  1:1223396/90‑1 (MQ=255)
                                                           aGCATCACGACCACGTTGATGCAGAACGGACAGGcc                                                                               >  2:878793/1‑36 (MQ=255)
                                                           aGCATCACGACCACGTTGATGCAGAACGGACAGGcc                                                                               <  1:878793/36‑1 (MQ=255)
                                                           aGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGa                                           >  1:142982/1‑72 (MQ=255)
                                                           aGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGa                                           <  2:142982/72‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAAGGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCACATAGTTGATGGTGATGTAGc                        <  1:1436930/90‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAg                                                                             <  2:1430637/37‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAg                                                                             >  1:1430637/1‑37 (MQ=39)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        <  1:213448/90‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        <  1:1351049/90‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        <  2:258790/90‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGACGTAGc                        >  2:550431/1‑90 (MQ=255)
                                                                 aCGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGAtt                   <  1:589102/90‑1 (MQ=255)
                                                                 aCGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGAtt                   <  2:1397975/90‑1 (MQ=255)
                                                                 aCGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGAtt                   >  2:1220479/1‑90 (MQ=255)
                                                                   gACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTAt                 <  1:1396247/90‑1 (MQ=255)
                                                                   gACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTAt                 >  1:1193320/1‑90 (MQ=255)
                                                                      cACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGcc              >  2:215653/1‑90 (MQ=255)
                                                                      cACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGcc              >  2:1436440/1‑90 (MQ=255)
                                                                       aCGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCa             >  2:32580/1‑90 (MQ=255)
                                                                       aCGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCa             >  1:1437827/1‑90 (MQ=255)
                                                                        cGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        <  1:88337/78‑1 (MQ=255)
                                                                        cGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        >  2:88337/1‑78 (MQ=255)
                                                                        cGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCaa            >  2:499300/1‑90 (MQ=255)
                                                                        cGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCaa            <  1:278096/90‑1 (MQ=255)
                                                                        cGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCaa            <  1:446865/90‑1 (MQ=255)
                                                                           tGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATcc         <  2:1356148/90‑1 (MQ=255)
                                                                           tGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATcc         <  2:1028644/90‑1 (MQ=255)
                                                                            gATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATccc        >  2:184247/1‑90 (MQ=255)
                                                                            gATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATccc        >  1:1025867/1‑90 (MQ=255)
                                                                                  gAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGc  >  2:1115437/1‑90 (MQ=255)
                                                                                       |                                                                                    
CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851991‑3852162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: