Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 470,313 C→G 7.2% I481M (ATC→ATG)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,3130CG 7.2% 173.5 / 3.1 53I481M (ATC→ATG) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (20/29);  new base G (2/2);  total (22/31)
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

GCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGC  >  NC_000913/470171‑470460
                                                                                                                                              |                                                                                                                                                   
gCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCAtt                                                                                                                                               >  8:205640/1‑149 (MQ=255)
 cccGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCAtt                                                                                                                                               <  6:126752/148‑1 (MQ=255)
  ccGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGc                                                                                                                                             >  2:396082/1‑149 (MQ=255)
        gCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGc                                                                                                                                       >  8:222598/1‑149 (MQ=255)
           aCCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGt                                                                                                                                     >  6:420257/1‑148 (MQ=255)
            cccAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGtta                                                                                                                                   <  6:452500/149‑1 (MQ=255)
             ccAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttat                                                                                                                                  <  1:197479/149‑1 (MQ=255)
             ccAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttat                                                                                                                                  <  7:146521/149‑1 (MQ=255)
             ccAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttat                                                                                                                                  <  2:369069/149‑1 (MQ=255)
               aGCAAGAGATTTAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttatta                                                                                                                                <  1:350558/149‑1 (MQ=255)
               aGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttatta                                                                                                                                <  6:413544/149‑1 (MQ=255)
               aGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttatta                                                                                                                                <  7:205640/149‑1 (MQ=255)
                gCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAg                                                                                                                               <  1:404001/149‑1 (MQ=255)
                 cAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGt                                                                                                                              <  4:172441/149‑1 (MQ=255)
                  aaGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGt                                                                                                                              >  1:23802/1‑148 (MQ=255)
                      gaTTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGGCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAAcg                                                                                                                         <  3:379535/149‑1 (MQ=255)
                           aGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGaa                                                                                                                    <  7:222598/149‑1 (MQ=255)
                             cATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAAt                                                                                                                  >  5:257244/1‑149 (MQ=255)
                                 tCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCt                                                                                                               <  6:257244/148‑1 (MQ=255)
                                   gcgcgGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCAt                                                                                                            <  1:396082/149‑1 (MQ=255)
                                       ggTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATGTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCtt                                                                                                        >  8:32694/1‑149 (MQ=255)
                                       gcttAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATGTCCATAGCTCGTGCGTTATTAGTCTACGCGGAAATCCTCATCCtt                                                                                                        >  8:39749/3‑149 (MQ=255)
                                         ttAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTtga                                                                                                      <  5:420257/149‑1 (MQ=255)
                                          tAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTtgat                                                                                                     >  3:130968/1‑149 (MQ=255)
                                              cAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGatgatg                                                                                                 <  1:342521/149‑1 (MQ=255)
                                              cAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTAATcc                                                                                                          <  4:383080/140‑1 (MQ=255)
                                              cAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTAATcc                                                                                                          >  3:383080/1‑140 (MQ=255)
                                                       tGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTcgg                                                                                        >  8:467370/1‑149 (MQ=255)
                                                           gATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCACCGCGGAAATCCTCATCCTAGATGATGCGCTTTCGGCGGt                                                                                    >  6:331918/1‑149 (MQ=255)
                                                           gATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGt                                                                                    <  7:82449/149‑1 (MQ=255)
                                                                    cGTCTACCGAAAGGTTACGAAACAGAGGTGGGCGAGCGCGTTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTACATTGATCGTGCGTTATTAGTCAACGCGGAAATCATCATCCTTGATGATGCGCTTTCGGCGGTggacggac                                                                            <  5:237630/148‑1 (MQ=255)
                                                                    cGTCTACCCCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGGGATGCTGTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGGGCGGTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTggacggac                                                                            <  5:331918/148‑1 (MQ=255)
                                                                       cTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGtt                                                                                                                                    >  2:356058/1‑89 (MQ=255)
                                                                       cTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGtt                                                                                                                                    <  1:356058/89‑1 (MQ=255)
                                                                        tACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttatta                                                                                                                                <  5:242334/92‑1 (MQ=255)
                                                                        tACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGttatta                                                                                                                                >  6:242334/1‑92 (MQ=255)
                                                                           cGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTg                                                                      <  2:249562/147‑1 (MQ=255)
                                                                           cGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCCGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTg                                                                      >  1:249562/1‑147 (MQ=255)
                                                                                ggTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCa                                                               <  2:140854/149‑1 (MQ=255)
                                                                                      gATACAGAGGTGGGCGAGCGCGGTGTGATGCTTACCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGTTGGACGGACGCACTGAGCACCAGATCCt                                                         <  1:115522/149‑1 (MQ=255)
                                                                                        tACAGAGGTGGGAGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGc                                                       <  4:130968/149‑1 (MQ=255)
                                                                                                     gAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGc                                                       <  7:12900/136‑1 (MQ=255)
                                                                                                     gAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGc                                                       >  8:12900/1‑136 (MQ=255)
                                                                                                           ggTGTTATGCTTACCCGAGTTCAAAAACAGCGTAACTCCATTCCTCGAGCGTTATTAGTCAACGCGCAAATACTCAACCTTGATGAAGAGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTgcgg                                    >  2:40415/1‑149 (MQ=255)
                                                                                                                 aTGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAgg                               >  5:118501/1‑148 (MQ=255)
                                                                                                                 aTGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAgg                               >  6:165010/1‑148 (MQ=255)
                                                                                                                   gCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTtga                                                                                                      >  3:330282/1‑75 (MQ=255)
                                                                                                                     tttCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGAt                                                            <  4:285924/115‑1 (MQ=255)
                                                                                                                     tttCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGAt                                                            >  3:285924/1‑115 (MQ=255)
                                                                                                                         cggcggGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACgg                      <  7:467370/149‑1 (MQ=255)
                                                                                                                              ggCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCGTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCTTAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTatt                  >  6:381031/1‑146 (MQ=255)
                                                                                                                                           tATGTCCATTGCTCGTGCGTTATTTGTCAACGCGGAAATGCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCGCATc    <  7:39749/149‑1 (MQ=255)
                                                                                                                                           tATGTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATc    <  7:32694/149‑1 (MQ=255)
                                                                                                                                             tctcCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGc  >  1:174757/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGC  >  NC_000913/470171‑470460

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