Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,138,452 T→G 6.9% L171R (CTG→CGG)  yadG → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,138,4520TG 6.9% 137.6 / 2.6 52L171R (CTG→CGG) yadGABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base G (2/2);  total (26/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 = 9.72e-01

TCAGCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCAT  >  NC_002947/2138306‑2138599
                                                                                                                                                  |                                                                                                                                                   
tCAGCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTgcg                                                                                                                                                  >  1:316337/1‑150 (MQ=255)
 cAGCTGGGCCTGTGGGACAAACGTGATTTGCAGTCGTGTTCTTTTTCCGGCTGCATGAAGCGGCGCCTTATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTTGACGATCCCACTGCCGGGGTGGATATCGAGCTgcgc                                                                                                                                                 <  2:585043/150‑1 (MQ=255)
          cTGTGGGACAAACGTGGTGTGATGTCGCGTTCGTTGGCCGGCTGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAta                                                                                                                                        <  2:336312/150‑2 (MQ=255)
          cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtg                                                                                                                                        >  1:213587/1‑150 (MQ=255)
               ggACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAg                                                                                                                                   >  1:1001095/1‑150 (MQ=255)
                  cAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtgtg                                                                                                                                      <  2:667508/144‑1 (MQ=255)
                  cAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtgtg                                                                                                                                      >  1:667508/1‑144 (MQ=255)
                  cAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCtt                                                                                                                                <  2:213587/150‑1 (MQ=255)
                   aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTc                                                                                                                               <  1:59935/150‑1 (MQ=255)
                         gATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACt                                                                                                                         <  2:90661/150‑1 (MQ=255)
                            gtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGa                                                                                                                                    >  2:919840/1‑136 (MQ=255)
                            gtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGa                                                                                                                                    <  1:919840/136‑1 (MQ=255)
                               cAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTg                                                                                                                   >  2:52854/1‑150 (MQ=255)
                                 gTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGaa                                                                                                                 <  2:1219719/150‑1 (MQ=255)
                                      gttcgttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGa                                                                                                            >  1:1047270/1‑150 (MQ=255)
                                              tCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATc                                                                                                    >  2:834541/1‑150 (MQ=255)
                                                       aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                           <  1:52854/150‑1 (MQ=255)
                                                           aGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTca                                                                                       <  1:375888/150‑1 (MQ=255)
                                                            gcggcgCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTcac                                                                                      <  2:1001095/150‑1 (MQ=255)
                                                                cgcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCaccacc                                                                                  <  2:82641/150‑1 (MQ=255)
                                                                cgcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCaccacc                                                                                  <  1:1195171/150‑1 (MQ=255)
                                                                 gcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGc                                                                                                                                  >  1:1006277/1‑101 (MQ=255)
                                                                 gcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGc                                                                                                                                  <  2:1006277/101‑1 (MQ=255)
                                                                 gcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCAccaccc                                                                                 >  2:376589/1‑150 (MQ=255)
                                                                 gcCTGATGATTGCCCGCGCATTGAGCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCAccaccc                                                                                 >  1:671397/1‑150 (MQ=255)
                                                                    tgatgaTTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAg                                                                                                          >  1:624189/1‑122 (MQ=255)
                                                                    tgatgaTTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAg                                                                                                          <  2:624189/122‑1 (MQ=255)
                                                                     gatgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTa                                                                             >  1:152919/1‑150 (MQ=255)
                                                                         atTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTg                                                                         <  1:376589/150‑1 (MQ=255)
                                                                             cccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGa                                                                                                                       <  1:644124/100‑1 (MQ=255)
                                                                             cccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGa                                                                                                                       >  2:644124/1‑100 (MQ=255)
                                                                              ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACTCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCCCCACCCACTATCTGGAcga                                                                    >  2:386806/1‑150 (MQ=255)
                                                                              ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCAAACTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCCTCACCATCATCCTCACCACCCACTATCTGGAaga                                                                    <  1:1265363/150‑1 (MQ=255)
                                                                                     ttGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGag                                                             <  2:772881/150‑1 (MQ=255)
                                                                                        aTCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGagcag                                                          <  2:1047270/150‑1 (MQ=255)
                                                                                         tCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGagcagc                                                         >  1:627181/1‑150 (MQ=255)
                                                                                           cACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtg                                                       >  2:25706/1‑150 (MQ=255)
                                                                                           cACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtg                                                       <  1:834541/150‑1 (MQ=255)
                                                                                                   gcgcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAgg                                                                  >  2:703926/1‑131 (MQ=255)
                                                                                                   gcgcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAgg                                                                  <  1:703926/131‑1 (MQ=255)
                                                                                                           tgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGca                                       >  2:718500/1‑150 (MQ=255)
                                                                                                                 tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                                 >  2:731458/1‑150 (MQ=255)
                                                                                                                   gACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGAc                               >  2:162472/1‑150 (MQ=255)
                                                                                                                       accccACCGCGGGGGTGGATATCGAGCGGCGGCGCTCGGTGTGGAGCGTCCTCACTGAGCTGACCCAGAAGGGCATCATCATCATCCTCACCCCCCACTATCTGGAAGAGGCGGAGCAGCTGGGCCGTAACATCGGCATCATCGACCATg                           <  1:386806/148‑1 (MQ=255)
                                                                                                                        aCCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATgg                          <  2:627181/150‑1 (MQ=255)
                                                                                                                                  ggggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCAtcgtc                <  1:25706/150‑1 (MQ=255)
                                                                                                                                     gTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCgag             >  2:294571/1‑150 (MQ=255)
                                                                                                                                       ggATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGaa           <  2:671397/150‑1 (MQ=255)
                                                                                                                                       ccacatCGAGCGGCGCAGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGCGGGTGGGTTGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGaa           >  2:675454/5‑150 (MQ=255)
                                                                                                                                       ccacatCGAGCGGCGCAGCTCGATATGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGaa           >  2:546994/5‑150 (MQ=255)
                                                                                                                                          tatCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAacac        >  1:1123094/1‑150 (MQ=255)
                                                                                                                                                gCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCAt  >  2:583013/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCAGCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCAT  >  NC_002947/2138306‑2138599

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