Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,568,563 A→C intergenic (‑44/‑435) mdtM ← / → yjiP multidrug efflux system protein/pseudogene, transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,568,5630AC98.2% 169.1 / NA 56intergenic (‑44/‑435)mdtM/yjiPmultidrug efflux system protein/pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (29/26);  total (30/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 = 1.00e+00

GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAT  >  NC_000913/4568419‑4568708
                                                                                                                                                |                                                                                                                                                 
ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCaa                                                                                                                                             <  1:444752/151‑1 (MQ=255)
 gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCAACGGGGGAAACGACATTCaa                                                                                                                                             >  2:908572/1‑150 (MQ=255)
  cTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGGGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAAAGACATTCAAt                                                                                                                                            >  2:785714/1‑150 (MQ=255)
   tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCCACGGGGGAAACGACATTCAATa                                                                                                                                           >  1:184883/1‑150 (MQ=255)
        cAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTTTCCGGGGAACCCCAACGGGGAAAACGACATTCAATACCGTc                                                                                                                                      >  2:107453/1‑150 (MQ=255)
        cAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCg                                                                                                                                     <  1:850532/151‑1 (MQ=255)
         aGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACATGGCATAACAAAAAACACCCTTTCCGGAGAACCCCCACGGGGGAAACGACATTCAATACCGTccg                                                                                                                                    >  1:615208/1‑149 (MQ=255)
            tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTa                                                                                                                                  <  1:643762/150‑1 (MQ=255)
             ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAAt                                                                                                                                <  1:922559/151‑1 (MQ=255)
                   aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGAGGCATAACAAACAAATCCTTGTCCGGAGAACCCCAACGGGGAAAAAGACATTTAATACCGTCGGGAATGACg                                                                                                                            >  2:8089/1‑149 (MQ=255)
                    cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGaaa                                                                                                                         <  2:857401/151‑1 (MQ=255)
                    cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAAAAAATCCTTGTCCGGAGAACCCCCACGGGGAAAACGACATTTAATACCCTCGGTAATgaagaa                                                                                                                          >  2:167078/1‑150 (MQ=255)
                           gcagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGt                                                                                                                  <  2:655521/151‑1 (MQ=255)
                             agcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTga                                                                                                                <  1:342886/151‑1 (MQ=255)
                                  aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAgg                                                                                                           <  2:104618/151‑1 (MQ=255)
                                  aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGAGGCAAAAAAAAAAAAACCTTGTTCGGAGAACACCAACGGGGGAAACGAAAATCAAAACCGTCGGTAATGACGAAAGGAAGGGGAGaa                                                                                                             >  1:675939/1‑149 (MQ=255)
                                    tCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACCTGGCATAACAAAAAACTCCTTGTCCGGGGAACCCCAACGGGGGAAAAGACATTTAATAACGTCGGGAAAGACGAAAAGCAGGGGAGAAggg                                                                                                          >  2:563210/1‑150 (MQ=255)
                                      aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCaa                                                                                                       <  2:13869/151‑1 (MQ=255)
                                      aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAAAAAACAACTCCTTTTCCGGAGAAACCCAAAGGGGGAAAAGACATTCAATAACCTCGGTAATGACGAAAGGCAGGGGAGAAggg                                                                                                          >  1:432112/1‑148 (MQ=255)
                                             aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACt                                                                                                <  2:429489/151‑1 (MQ=255)
                                             aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGGAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGGGAGAAGGGCAACGAGACt                                                                                                >  2:452649/1‑151 (MQ=255)
                                               cAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGGATAACAAACAACCCCCTTTCCGGGGAAACCCAACGGGGGAAAAGACAATCAATACCGTCGGGAATGAAGAAAAGCAGGGGGGAAGGGGAACGAGACTg                                                                                               >  1:720091/1‑150 (MQ=255)
                                                 aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCCACGGGGGAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGGGAGAAGGGCAACGAGACTGc                                                                                              >  1:188381/1‑149 (MQ=255)
                                                  cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaaa                                                                                           <  1:657523/151‑1 (MQ=255)
                                                   gCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAAACCCAACGGGGAAAACGACATTCAATACCGTCGGGAATGACGAAAGGCAGGGGAGAAGGGCAACGAGACTGCaaa                                                                                           >  2:158054/1‑150 (MQ=255)
                                                     cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGt                                                                                        <  1:533863/151‑1 (MQ=255)
                                                     cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGGAATGACGAAAGGCAGGGGAGAAGGGCAACGAGACTGCaaa                                                                                           >  2:860726/1‑148 (MQ=255)
                                                            aaaaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAAAAAAAAACTCCCTGTTCGGAGAACCCCAACGGGGAAAACGACATTCAATAACGTCGGGGATGACGAAAAGCAGGGGAGAAGGGCAACGAGACTGCAAAAGTACTGaa                                                                                  >  1:867303/1‑150 (MQ=255)
                                                             aaaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAAtt                                                                                <  1:515182/151‑1 (MQ=255)
                                                               aaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTc                                                                              <  2:351407/151‑1 (MQ=255)
                                                                aCAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAAAAAATCCTTGTCCGGGGAACCCCAACGGGGAAAAAGACATTCAATACCGTCGGGAATGACGAAAAGCAGGGGGGAAGGGCAACGAGAATGCAAAAGTAATGAAttt                                                                               >  1:807242/1‑149 (MQ=255)
                                                                 cAGCGTGGCGGCATGGCGGGGCAAAAAAAGTGGCATAACAAAAAACACCTTTTCTGGAGAAACCCCACGGGGGAAAAGACATTCAATTACGTCGGTAAAAATGGAAGTCGCGTGAGGAGGGGAAAGAGACCTCAAATGTACTGGAtttta                                                                             >  1:834776/1‑148 (MQ=255)
                                                                        gcggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCa                                                                     <  1:126638/151‑1 (MQ=255)
                                                                               ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCacac                                                              <  2:362220/151‑1 (MQ=255)
                                                                                gCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACt                                                             >  1:233446/1‑151 (MQ=255)
                                                                                                ggCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACacag                                             <  1:690900/151‑1 (MQ=255)
                                                                                                 gCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAacacac                                               >  1:855048/1‑148 (MQ=255)
                                                                                                 gCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACacaga                                            <  2:500200/151‑1 (MQ=255)
                                                                                                     aacaaacaACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACacagacag                                         <  2:785379/150‑1 (MQ=255)
                                                                                                     aacaaacaACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACacagac                                           >  2:477763/1‑148 (MQ=255)
                                                                                                     aacaaacaACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACacagac                                           >  1:657367/1‑148 (MQ=255)
                                                                                                     aacaaacaACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACAcagacaga                                        <  1:689010/151‑1 (MQ=255)
                                                                                                       caaacaaCTCCTTGTCCGGAGAACCCCAACGGGGCAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACacagacag                                         >  2:88044/1‑148 (MQ=255)
                                                                                                        aaacaaCTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGaa                                     <  1:341924/151‑1 (MQ=255)
                                                                                                          acaacTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAAt                                   >  1:258118/1‑151 (MQ=255)
                                                                                                          acaacTCCTTGTCCGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAAt                                   <  2:922618/151‑1 (MQ=255)
                                                                                                                      ccGGAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACtt                       >  2:186222/1‑151 (MQ=255)
                                                                                                                        ggAGAACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCg                     <  2:530602/151‑1 (MQ=255)
                                                                                                                          agaACCCCAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCggg                   >  1:772512/1‑151 (MQ=255)
                                                                                                                              ccccAACGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTc                 >  2:477934/1‑149 (MQ=255)
                                                                                                                                   aCGGGGAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGt          >  1:390248/1‑151 (MQ=255)
                                                                                                                                      gggAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGtat        >  2:272586/1‑150 (MQ=255)
                                                                                                                                       ggAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATaa      <  2:351840/151‑1 (MQ=255)
                                                                                                                                        gAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATaa      <  2:607713/150‑1 (MQ=255)
                                                                                                                                        gAAAACGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAAc     >  2:332803/1‑151 (MQ=255)
                                                                                                                                           aaCGACATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAt  <  2:671257/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAT  >  NC_000913/4568419‑4568708

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑