Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,160,092 T→G S569R (AGT→AGG mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,160,0920TG98.0% 171.2 / ‑4.1 50S569R (AGT→AGGmdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/29);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

GCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCACCCTCAAGCCACGCG  >  NC_000913/2159951‑2160234
                                                                                                                                             |                                                                                                                                              
gCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGttt                                                                                                                                       <  2:655138/151‑1 (MQ=255)
   cAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGTACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAg                                                                                                                                    <  1:547048/151‑1 (MQ=255)
   cAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAg                                                                                                                                    <  2:887483/151‑1 (MQ=255)
    aaTCATGCCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGttt                                                                                                                                       <  1:356270/147‑1 (MQ=255)
     aTCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGc                                                                                                                                  <  2:182245/151‑1 (MQ=255)
          aCCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATgc                                                                                                                             >  2:439650/1‑151 (MQ=255)
          aCCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATgc                                                                                                                             >  1:302688/1‑151 (MQ=255)
                cTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGt                                                                                                                        >  1:531341/1‑150 (MQ=255)
                cTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGt                                                                                                                        >  1:692691/1‑150 (MQ=255)
                 tggtggGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTa                                                                                                                       >  1:286297/1‑150 (MQ=255)
                     gggCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGtt                                                                                                                   <  2:755474/150‑1 (MQ=255)
                      ggCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGc                                                                                                                 <  1:634260/151‑1 (MQ=255)
                        cGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTg                                                                                                                  >  1:518192/1‑148 (MQ=255)
                         gtggtgCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAgg                                                                                                              <  1:692954/151‑1 (MQ=255)
                           ggtgCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGAt                                                                                                            <  1:474570/151‑1 (MQ=255)
                                 gctTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATg                                                                                                      <  2:224522/151‑1 (MQ=255)
                                   ttGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTAGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGa                                                                                                     >  1:311144/1‑150 (MQ=255)
                                        aCCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAAtt                                                                                                >  2:92383/1‑150 (MQ=255)
                                        aCCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAAtt                                                                                                >  1:790658/1‑150 (MQ=255)
                                              gcgcTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTg                                                                                         <  1:553912/151‑1 (MQ=255)
                                               cgcTGAATATCTGGCTGTATATCTCGATCCCGAAATCCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGa                                                                                        <  1:933595/151‑1 (MQ=255)
                                                gcTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGAc                                                                                       >  1:936478/1‑151 (MQ=255)
                                                  tGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACg                                                                                      >  1:545106/1‑150 (MQ=255)
                                                  tGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACg                                                                                      >  1:299038/1‑150 (MQ=255)
                                                     atatCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATcc                                                                                  <  1:268040/151‑1 (MQ=255)
                                                          tGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGc                                                                               >  1:227021/1‑149 (MQ=255)
                                                                 atatCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATaa                                                                      <  1:7530/151‑1 (MQ=255)
                                                                 atatCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATa                                                                       <  1:582462/150‑1 (MQ=255)
                                                                   atCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATg                                                                    >  1:404709/1‑151 (MQ=255)
                                                                           cccGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGc                                                            <  1:701013/151‑1 (MQ=255)
                                                                            ccGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCt                                                           >  2:695807/1‑151 (MQ=255)
                                                                             cGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCtt                                                          >  1:224770/1‑151 (MQ=255)
                                                                                      tcttcCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCCGGCGGATCAGAGGATTTTGTTTCAGGCGATGCGCGGTCAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCgg                                                 >  2:398270/1‑151 (MQ=255)
                                                                                      tcttcCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGc                                                   <  2:861418/149‑1 (MQ=255)
                                                                                           ccGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGtttg                                             >  1:151248/1‑148 (MQ=255)
                                                                                           ccGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTcgc                                            <  1:781151/151‑1 (MQ=255)
                                                                                                gCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAg                                         >  2:678790/1‑149 (MQ=255)
                                                                                                  aGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGaa                                     <  2:33551/151‑1 (MQ=255)
                                                                                                   ggACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAAc                                    <  1:409704/151‑1 (MQ=255)
                                                                                                      cacTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGc                                 <  1:194563/151‑1 (MQ=255)
                                                                                                         tGGCGTGTTGGTGGGCGGGAGTCAGGCGGATCTGAGGCTTTCGTTTCAGGCGATGCGCGCTAAGTTGCTGGATTTCATGAAAATTATCGGTGACGATCCGGCAGTGGATAATGTCGCCGGCTTTACAGGCGGTTCGCGAGTGAACAGc                                 <  1:334267/148‑1 (MQ=255)
                                                                                                          ggCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCggg                              <  2:587403/150‑1 (MQ=255)
                                                                                                          ggCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCggg                              <  2:97545/150‑1 (MQ=255)
                                                                                                                   aTGGGCGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTa                    <  1:883386/151‑1 (MQ=255)
                                                                                                                        cGGGATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCAcc               <  1:630315/151‑1 (MQ=255)
                                                                                                                          ggATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCACCCt             <  1:251875/151‑1 (MQ=255)
                                                                                                                           gATTCAGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCACCCTc            <  2:928612/151‑1 (MQ=255)
                                                                                                                                aGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCACCCTCAAGcc       <  1:349922/151‑1 (MQ=255)
                                                                                                                                aGGCGGATCAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGGGGACAGCGGGAAGATGTTTATTAACCTCAAGcc       >  1:586374/1‑151 (MQ=255)
                                                                                                                                        cAGAGGATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCACCCTCAAGCCAcgcg  >  2:233834/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGGCTTTACAGGCGGTTCGCGAGTGAACAGCGGGATGATGTTTATCACCCTCAAGCCACGCG  >  NC_000913/2159951‑2160234

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