Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,193 G→A 7.8% intergenic (+16/+52) ynfM → / ← PP_5073 efflux transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,1930GA 7.8% 172.8 / 2.8 60intergenic (+16/+52)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (3/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGGCGGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTT  >  NC_002947/5790054‑5790337
                                                                                                                                           |                                                                                                                                                
ggcgggcgGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGTCACAcc                                                                                                                                                <  1:1129741/142‑1 (MQ=255)
         ggTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCACGCACCCACAACCtgtg                                                                                                                               <  2:723950/150‑1 (MQ=255)
          gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                              >  2:762599/1‑150 (MQ=255)
          gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                              <  1:717658/150‑1 (MQ=255)
                   cACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtg                                                                                                                     <  2:1218303/150‑1 (MQ=255)
                   cACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGATGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTAATGCATTCGTGGGAACACCCGCACCCACAACCTGTGAGAGCGGgtg                                                                                                                     >  2:904741/1‑150 (MQ=255)
                    aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgcgt                                                                                                                    <  1:527848/150‑1 (MQ=255)
                             gCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                        <  2:278953/137‑1 (MQ=255)
                             gCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                        >  1:278953/1‑137 (MQ=255)
                                  aaCGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGc                                                                                                      <  1:1159892/150‑1 (MQ=255)
                                     gggATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCat                                                                                                   >  1:698389/1‑150 (MQ=255)
                                      ggATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                  <  1:279718/150‑1 (MQ=255)
                                      ggATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                  <  1:881010/150‑1 (MQ=255)
                                       gATCGGGCTGTTCATTGGCAGCTTGCTTTCTGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatca                                                                                                 <  1:904741/150‑1 (MQ=255)
                                                 ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtga                                                                                       >  1:1247725/1‑150 (MQ=255)
                                                      tGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtc                                                                                  <  1:424105/150‑1 (MQ=255)
                                                        gCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTGGATGAtctc                                                                                >  1:740669/1‑150 (MQ=255)
                                                            cTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAAACTGTGGGAGc                                                                                                                          <  2:672990/104‑1 (MQ=255)
                                                            cTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAAACTGTGGGAGc                                                                                                                          >  1:672990/1‑104 (MQ=255)
                                                               gCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGc                                                                                                                                                 >  1:802788/1‑78 (MQ=255)
                                                               gCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGc                                                                                                                                                 <  2:802788/78‑1 (MQ=255)
                                                               gCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAata                                                                         <  1:1166010/150‑1 (MQ=255)
                                                                  ggCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcc                                                                      <  1:63450/150‑1 (MQ=255)
                                                                   gCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcca                                                                     >  1:135250/1‑150 (MQ=255)
                                                                    cGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccac                                                                    <  2:785049/150‑1 (MQ=255)
                                                                     ggTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                         <  2:425480/128‑1 (MQ=255)
                                                                     ggTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                         >  1:425480/1‑128 (MQ=255)
                                                                          ccTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCca                                                                                                                                        <  1:1151394/76‑1 (MQ=255)
                                                                          ccTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCca                                                                                                                                        >  2:1151394/1‑76 (MQ=255)
                                                                          ccTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacatccac                                                              >  2:1143673/1‑150 (MQ=255)
                                                                           cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGcgtgcc                                                                                                                 >  2:781666/1‑96 (MQ=17)
                                                                           cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatca                                                                                            >  2:1135662/1‑119 (MQ=255)
                                                                           cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatca                                                                                            <  1:1135662/119‑1 (MQ=255)
                                                                           cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                             >  2:1106493/1‑150 (MQ=255)
                                                                           cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                             >  1:149270/1‑150 (MQ=255)
                                                                              gATGGTGGCGTTGCAGCTGAGCAAGTTGCCACCAAAAACGGTGTGATGCATTCGTGTGCACGCCCGCACCCACAACCTGTGGCAGCGGGTGTGACCACGCAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTcc                                                          <  2:516524/150‑1 (MQ=255)
                                                                               aTGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCa                                                                                                               <  2:186444/96‑1 (MQ=255)
                                                                               aTGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCa                                                                                                               >  1:186444/1‑96 (MQ=255)
                                                                                     gCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCAAAATCCACGTCCCGGCtca                                                   >  2:873640/1‑150 (MQ=255)
                                                                                              cTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCt                                          <  2:1247725/150‑1 (MQ=255)
                                                                                                 aGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAag                                                                                                          >  1:171005/1‑83 (MQ=255)
                                                                                                 aGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAag                                                                                                          <  2:171005/83‑1 (MQ=255)
                                                                                                   cAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTc                                      <  1:615196/149‑1 (MQ=255)
                                                                                                   cAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTc                                      >  2:615196/1‑149 (MQ=255)
                                                                                                           cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAcc                             >  1:1100312/1‑150 (MQ=255)
                                                                                                              ccAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGt                          <  2:135250/150‑1 (MQ=255)
                                                                                                              ccAAAAAGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGt                          >  2:851489/1‑150 (MQ=255)
                                                                                                                aacaCGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGCGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCTCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCa                        <  1:361034/147‑1 (MQ=255)
                                                                                                                aaaaCGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgatg                                                                                     <  2:728450/89‑1 (MQ=255)
                                                                                                                aaaaCGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgatg                                                                                     >  1:728450/1‑89 (MQ=255)
                                                                                                                  aaCGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGc                      <  2:698389/150‑1 (MQ=255)
                                                                                                                          gATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGAGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCtt              <  1:873640/150‑1 (MQ=255)
                                                                                                                             gCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                 >  2:40389/1‑80 (MQ=255)
                                                                                                                             gCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                 <  1:40389/80‑1 (MQ=255)
                                                                                                                              tcttCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTc          >  2:724869/3‑150 (MQ=25)
                                                                                                                               aTTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCt         >  1:1018523/1‑150 (MQ=255)
                                                                                                                               aTTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCt         >  2:99325/1‑150 (MQ=255)
                                                                                                                                    tGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTAACCTTGGTCtt        <  1:256824/146‑1 (MQ=25)
                                                                                                                                     gggCATGCCCGCACCCACAACCTGTGGGAGCGGGTGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGt   <  1:367067/150‑1 (MQ=255)
                                                                                                                                      ggCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGtt  <  1:770547/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                
GGCGGGCGGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTT  >  NC_002947/5790054‑5790337

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