Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,721,377 G→A 18.2% G38G (GGG→GGA ydhJ → putative membrane fusion protein (MFP) of YdhJK efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,721,3770GA18.2% 166.6 / 29.3 77G38G (GGG→GGAydhJputative membrane fusion protein (MFP) of YdhJK efflux pump
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base A (6/8);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

ATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCT  >  NC_000913/1721279‑1721475
                                                                                                  |                                                                                                  
aTTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGa                                                                                                   >  1:1070700/1‑100 (MQ=255)
  tAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGaaa                                                                                                 <  2:349829/100‑1 (MQ=255)
  tAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGaaa                                                                                                 <  2:349764/100‑1 (MQ=255)
   aaGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGaaaa                                                                                                >  1:2068262/1‑100 (MQ=255)
      tATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATAc                                                                                             <  1:1903333/100‑1 (MQ=255)
      tATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATAc                                                                                             <  2:2535700/100‑1 (MQ=255)
       tttcctaaacaaTAATTGTAAGGGGATTTTCGGTTCTTGCCGGATGGTGTTTGTGGAATTATTACATGCGATCTCCGTGGAGACGGGACGGGaaaa                                                                                                <  1:123711/89‑1 (MQ=255)
          tctcAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCGGGATGGTGGCTGTGGAATTATTACAGGCAATCGTCGTGGTCACGCGACGGGAAAATAcgc                                                                                           >  2:380376/1‑98 (MQ=255)
          tctcAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACgcgc                                                                                         >  1:2362099/1‑100 (MQ=255)
          tctcAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCGGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACgcgc                                                                                         <  1:1862416/100‑1 (MQ=255)
            tcAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAg                                                                                       <  2:874226/100‑1 (MQ=255)
             caacaaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGa                                                                                      <  2:1884125/100‑1 (MQ=255)
                caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACa                                                                                   >  2:643388/1‑100 (MQ=255)
                caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACa                                                                                   <  2:2167897/100‑1 (MQ=255)
                caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACa                                                                                   <  1:766131/100‑1 (MQ=255)
                caaTCATTGTAACGGTAATTACGGTTCTCGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACa                                                                                   >  1:1353090/1‑100 (MQ=255)
                 aaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAg                                                                                  >  2:881333/1‑100 (MQ=255)
                  aTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAgg                                                                                 <  2:2212911/100‑1 (MQ=255)
                  aTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAgg                                                                                 <  2:1063732/100‑1 (MQ=255)
                     aTTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGttt                                                                              >  2:694727/1‑100 (MQ=255)
                     aTTGTAGCGGTAAGTGCGGTTCTTGACGGATGGTGGCTGTGGAACTAGTACATGCAGTCTGCGTGGACACGCGACGGGAAAATACGCTCATAACAGGttt                                                                              >  2:2583739/1‑100 (MQ=255)
                      ttGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTc                                                                             <  1:254217/100‑1 (MQ=255)
                      ttGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTc                                                                             <  2:577929/100‑1 (MQ=255)
                        gTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCACAGAACAGGTTTCTa                                                                           >  2:2549729/1‑100 (MQ=255)
                           gCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCa                                                                        >  2:1066607/1‑100 (MQ=255)
                             ggTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCAcc                                                                      >  2:781668/1‑100 (MQ=255)
                             ggTAGTTGCGGTACTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGCCACGCCACGGGGAAATACGCGCAGAACAGGTTTCTATCAcc                                                                      >  2:1702516/1‑100 (MQ=255)
                               tAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCAcccc                                                                    <  2:1070704/100‑1 (MQ=255)
                                 gtagaggTTTTTGACGGATGGTGGCGGTGGAATTATTACGTGCAAAATCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCcac                                                                  <  1:839607/95‑1 (MQ=255)
                                 gTTGCGGTTCTTGCCGGATGGTGGCTGTGGAAGTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCcac                                                                  <  1:582848/100‑1 (MQ=255)
                                  ttGCGGTTCTTGCCGGATGGTGTCTGTGGAATTATTACACGAAATCTCCGTGGACACGCGACGTGAAAATACGCGCAGAACAGTTTTCTATCAccc                                                                     <  1:360750/96‑1 (MQ=255)
                                   tGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAg                                                                >  1:1515166/1‑100 (MQ=255)
                                    gCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTGCGCGGACACGGGACGGGAAAATACCCGCAGACTCGGTTTCTATCACCCaaca                                                                 >  1:2296123/1‑98 (MQ=255)
                                    gCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAgg                                                               >  2:1481802/1‑100 (MQ=255)
                                      ggTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGgtg                                                             >  1:2438538/1‑100 (MQ=255)
                                      ggTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGgtg                                                             >  2:2404059/1‑100 (MQ=255)
                                        ttCTTGCCGGATGGTGGCTGTGAAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTc                                                           <  1:1371391/100‑1 (MQ=255)
                                            tGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGa                                                       <  2:1307534/100‑1 (MQ=255)
                                            tGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGa                                                       <  2:1347110/100‑1 (MQ=255)
                                             gCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGAc                                                      >  1:815572/1‑100 (MQ=255)
                                             gCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGAc                                                      >  2:747091/1‑100 (MQ=255)
                                             gCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGAc                                                      <  1:2457452/100‑1 (MQ=255)
                                                 gATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTAt                                                  <  2:1023057/100‑1 (MQ=255)
                                                    ggtggCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTAttgt                                               >  1:929493/1‑100 (MQ=255)
                                                     gtggCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTAttgtt                                              >  1:1374457/1‑100 (MQ=255)
                                                       ggCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGa                                            >  1:577752/1‑100 (MQ=255)
                                                        gCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAg                                           <  1:1781903/100‑1 (MQ=255)
                                                        gCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAg                                           <  1:1091331/100‑1 (MQ=255)
                                                           gtgGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTg                                        >  1:715818/1‑100 (MQ=255)
                                                           gtgGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTg                                        >  1:2045782/1‑100 (MQ=255)
                                                             ggAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGaa                                      <  1:616646/100‑1 (MQ=255)
                                                             ggAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGaa                                      >  2:724706/1‑100 (MQ=255)
                                                              gAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAat                                     <  1:486410/100‑1 (MQ=255)
                                                              gAATTATTACATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAat                                     <  1:2458193/100‑1 (MQ=255)
                                                                attattACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAatat                                   <  2:2180238/100‑1 (MQ=255)
                                                                 ttattaCATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGATGTATTGTTGAGCTGAAtata                                  <  2:1329864/100‑1 (MQ=255)
                                                                   attaCATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATaaa                                <  2:243765/100‑1 (MQ=255)
                                                                    ttaCATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCACACAGGTGTCGGGCCGTATTGTTGAGCTGAATATaaaa                               >  2:1509172/1‑100 (MQ=255)
                                                                     taCATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAg                              <  1:2457400/100‑1 (MQ=255)
                                                                      aCATGCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGa                             >  2:2282541/1‑100 (MQ=255)
                                                                         tGCAATCTCCGTGGACACGCGACGGGAAAATACGCACAGAACAGGTTTCTATCACCCCAAAGGTGTCGGGACGTATTGGTGAGCTGCATATTAAAGAc                            >  1:999795/1‑98 (MQ=255)
                                                                          gCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAAc                         <  2:1564347/100‑1 (MQ=255)
                                                                          gCAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAAc                         >  2:1518149/1‑100 (MQ=255)
                                                                           cAATCTCCGTGGACACGCGACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAAcc                        <  2:1708247/100‑1 (MQ=255)
                                                                                tcCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTg                   <  2:493080/100‑1 (MQ=255)
                                                                                        cacGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAAcgc           <  2:2139461/100‑1 (MQ=255)
                                                                                        cacGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAAcgc           >  1:1491584/1‑100 (MQ=255)
                                                                                          cgcgACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGttgtctttt         <  2:695379/100‑8 (MQ=255)
                                                                                          cgcgACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCgg         >  2:2436747/1‑100 (MQ=255)
                                                                                          cgcgACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGCACATATTGTTGAGCAGAATATAAACAACAACCAGCAGGTTACCGCgg         >  1:1931839/1‑100 (MQ=255)
                                                                                          cgcgACGGGAAAATACGCGCAGAACAGGATTCTGTCACCCCACAGGTGTCGGGACGTATTGTTGCGCTGAATATAAAAGACACCCAGCTGGTTAAc             >  1:1696148/1‑96 (MQ=255)
                                                                                           gcgACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCggg        <  1:2404044/100‑1 (MQ=255)
                                                                                            cgACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCgggg       <  1:2566082/100‑1 (MQ=255)
                                                                                            cgACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCgggg       <  1:718517/100‑1 (MQ=255)
                                                                                            cgACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCgggg       >  2:1928875/1‑100 (MQ=255)
                                                                                            cgACGGAAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCgggg       <  2:1445042/100‑1 (MQ=255)
                                                                                              aCGGGAAAATACGCGCCGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGa     <  1:1514570/100‑1 (MQ=255)
                                                                                                 ggAAAATAGGCGCAGATCAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGAtc   <  1:1562351/99‑1 (MQ=255)
                                                                                                 ggAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGAtct  <  1:2486834/100‑1 (MQ=255)
                                                                                                 ggAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGAtct  <  2:2062528/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
ATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCT  >  NC_000913/1721279‑1721475

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: