Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,617,423 C→T R94C (CGC→TGC)  marR → MarR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,617,4230CT98.1% 285.3 / NA 103R94C (CGC→TGC) marRMarR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (45/56);  minor base G (0/2);  total (45/58)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCA  >  NC_000913/1617351‑1617493
                                                                        |                                                                      
gCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGtgcg                                                                     <  1:3412210/76‑1 (MQ=255)
  aCTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTgcggc                                                                   <  2:3658674/76‑1 (MQ=255)
  aCTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTgcggc                                                                   <  2:1373863/76‑1 (MQ=255)
  aCTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTgcggc                                                                   >  2:1931717/1‑76 (MQ=255)
   cTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTgcggcg                                                                  >  1:284327/1‑76 (MQ=255)
    tGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGt                                                                 >  1:3064268/1‑76 (MQ=255)
     gACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGtc                                                                >  2:5020847/1‑75 (MQ=255)
              gCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTaaaa                                                       >  1:604714/1‑76 (MQ=255)
              gCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTaaaa                                                       <  2:4140907/76‑1 (MQ=255)
              gCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTaaaa                                                       >  1:3320941/1‑76 (MQ=255)
               cTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAAc                                                      <  2:1765969/76‑1 (MQ=255)
                tGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACt                                                     >  1:5830965/1‑76 (MQ=255)
                 gggtCGACTGGTCTGTAAAGGCTGGGTGGAAAGTTTGCCGAACCCGCATGACAAGTGCGGCGTACTGGTAAAACtt                                                    <  2:1938815/73‑1 (MQ=255)
                 ggATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACtt                                                    <  2:5169699/76‑1 (MQ=255)
                  gATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTa                                                   <  2:1425279/76‑1 (MQ=255)
                     cGCCTGGTGTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAGCTTtcca                                                <  2:384062/76‑1 (MQ=255)
                     cGCCTGGTGTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTgcca                                                <  2:410569/76‑1 (MQ=255)
                     cGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTacca                                                <  2:3803389/76‑1 (MQ=255)
                      gCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTaccac                                               >  1:2975846/1‑76 (MQ=255)
                      gCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTaccac                                               >  2:2484852/1‑76 (MQ=255)
                       gcTGGGCTGGAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTaccacc                                              <  2:5430455/75‑1 (MQ=255)
                       ccTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTaccacc                                              >  1:4986935/1‑76 (MQ=255)
                       ccTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTaccacc                                              <  2:2682126/76‑1 (MQ=255)
                         tGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTCCCACcgg                                            <  1:1235788/76‑1 (MQ=255)
                         tGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcgg                                            >  1:3871523/1‑76 (MQ=255)
                         tGGTCTGTAAAAGCTGGGTGGAAAGGTTGCCGAAAACGAATGACAAGTGAGGCGTACTAGGAAAACTTACCACcgg                                            >  1:4513858/1‑76 (MQ=255)
                          ggTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcggc                                           >  1:364480/1‑76 (MQ=255)
                          ggTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcggc                                           <  2:5305330/76‑1 (MQ=255)
                          ggTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcggc                                           >  2:3865244/1‑76 (MQ=255)
                           gTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcggcg                                          >  1:2559864/1‑76 (MQ=255)
                           gTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcggcg                                          <  2:4785103/76‑1 (MQ=255)
                            tCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACcggcgg                                         >  1:337053/1‑76 (MQ=255)
                             cTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCggcggc                                        <  1:4024716/76‑1 (MQ=255)
                             cTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCggcggc                                        >  2:4294951/1‑76 (MQ=255)
                             cTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCggcggc                                        >  2:3612662/1‑76 (MQ=255)
                             cTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCggcggc                                        >  2:5374162/1‑76 (MQ=255)
                             cTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGGACGGGTAAAACTTACCACCggcggc                                        <  1:4061197/76‑1 (MQ=255)
                               gTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgc                                      <  2:3925903/76‑1 (MQ=255)
                               gTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgc                                      <  1:1563251/76‑1 (MQ=255)
                                tAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcg                                     >  1:1646733/1‑76 (MQ=255)
                                tAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcg                                     <  2:3660687/76‑1 (MQ=255)
                                 naaGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcacgg                                    >  2:2631961/2‑76 (MQ=255)
                                 aaaGGCTGGGTGGAAAGGTTTCCGAACCCGAATGACAAGTGCGGCGTACTGGTTAAACTTTCCACCGgcggcgcgg                                    <  2:3771731/76‑1 (MQ=255)
                                 aaaGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcgg                                    <  2:3741111/76‑1 (MQ=255)
                                 aaaGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcgg                                    >  2:5256030/1‑76 (MQ=255)
                                  aaGTCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcggc                                   <  2:3797555/76‑1 (MQ=255)
                                  aaGGCTGGGTGGACAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcggc                                   <  2:3950571/76‑1 (MQ=255)
                                  aaGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcggc                                   >  1:2342291/1‑76 (MQ=255)
                                  aaGGCTGGGTGGAAAGGTTGCCGAACCCGAAGGACAAGTGCGGCGTACTGGTAAAACTTACCACCGgcggcgcggc                                   <  2:408487/76‑1 (MQ=255)
                                   aGGCTGGGTAGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCa                                  <  1:4389804/76‑1 (MQ=255)
                                   aGGCTAGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCa                                  >  1:2932700/1‑76 (MQ=255)
                                    ggCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCaa                                 <  1:5602076/76‑1 (MQ=255)
                                    ggCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCaa                                 <  1:3058025/76‑1 (MQ=255)
                                     gCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTATTGGTAAAACTTACCACCGGCGGCGCGGCAat                                <  2:3977311/76‑1 (MQ=255)
                                     gCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAat                                <  1:3393180/76‑1 (MQ=255)
                                     gCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAat                                <  2:4238638/76‑1 (MQ=255)
                                     gCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGAGCGGCGat                                >  1:4184584/1‑76 (MQ=255)
                                       tGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTCCCACCGGCGGCGCGGCCatat                              >  2:409473/1‑76 (MQ=255)
                                       tGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAatat                              >  2:1725027/1‑76 (MQ=255)
                                       tGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAatat                              <  1:436376/76‑1 (MQ=255)
                                       tGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAatat                              <  2:2165479/76‑1 (MQ=255)
                                       tGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAatat                              >  1:4741172/1‑76 (MQ=255)
                                       tGGGTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAatat                              <  1:4823327/76‑1 (MQ=255)
                                        gggTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATAtg                             <  1:52869/76‑1 (MQ=255)
                                          gTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATAtgtg                           <  2:5038751/76‑1 (MQ=255)
                                          gTGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATAtgtg                           >  1:3428961/1‑76 (MQ=255)
                                          gTGGAAAGCCTGCCGAACCCGAATGACAAGTGCGGCCGACTGGTAAAACTTACCACCGGCGGCGCGGCAATAtgtg                           >  2:117026/1‑76 (MQ=255)
                                           tGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGa                          >  2:1218735/1‑76 (MQ=255)
                                           tGGAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGa                          >  1:3579916/1‑76 (MQ=255)
                                           gggAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGa                          <  2:504087/75‑1 (MQ=255)
                                            ggAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGaa                         <  1:1160732/76‑1 (MQ=255)
                                             gAAAGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAAc                        >  2:4327581/1‑76 (MQ=255)
                                              aaaGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACa                       <  2:5560027/76‑1 (MQ=255)
                                              aaaGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACa                       >  2:1515708/1‑76 (MQ=255)
                                              aaaGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACa                       <  2:2622081/76‑1 (MQ=255)
                                               aaGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACaa                      >  1:491738/1‑76 (MQ=255)
                                               aaGGTTCGCGAACCCGAATGACGAGGGCGGCGTACTGGTAAAACTTCCCACCGGCGGCGCGGCAATATGTGAACaa                      <  2:5658925/76‑1 (MQ=255)
                                                aGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAt                     >  1:1565973/1‑76 (MQ=255)
                                                aGGTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAt                     >  1:1604565/1‑76 (MQ=255)
                                                aGGTGGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAt                     <  2:2015503/76‑1 (MQ=255)
                                                  gTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGc                   <  1:4294952/76‑1 (MQ=255)
                                                  gTTGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGc                   <  1:5285653/76‑1 (MQ=255)
                                                   ttGCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGcc                  >  2:5573051/1‑76 (MQ=255)
                                                   tagCCGATCCCGAATGACAAGTGCGGCGTACTGGTAACACGTCCCGCCGGCGGCGCGGCAATATGTGAACTATGcc                  <  1:2563790/74‑1 (MQ=255)
                                                     gCCGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCAt                <  2:528781/76‑1 (MQ=255)
                                                      ccGAACCCGAATGACAAGGGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATc               <  2:2560880/76‑1 (MQ=255)
                                                      cagAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCCATATGTGAAAAATGCCATc               >  1:412618/3‑76 (MQ=255)
                                                       cGAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCa              >  1:4108727/1‑76 (MQ=255)
                                                        gAACCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCaa             <  1:5256031/76‑1 (MQ=255)
                                                          aCCCGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAAtt           <  2:2001810/76‑1 (MQ=255)
                                                            ccGAATGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAg         >  2:5349694/1‑76 (MQ=255)
                                                              gAAGGACAAGTGCGGCGGACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGtt       <  2:1414367/76‑1 (MQ=255)
                                                                aTGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTgg     >  1:1803056/1‑76 (MQ=255)
                                                                aTGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTgg     >  1:1573703/1‑76 (MQ=255)
                                                                aTGACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGCTgc     >  2:5718735/1‑75 (MQ=255)
                                                                 aggcAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGc    <  2:5654835/73‑1 (MQ=255)
                                                                  gACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGcc   <  1:3332316/76‑1 (MQ=255)
                                                                  gACAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGcc   <  2:4316341/76‑1 (MQ=255)
                                                                   aCAAGTGCGGCGTACTGGTAAACCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCa  <  2:2018176/76‑1 (MQ=255)
                                                                   aCAAGTGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCa  <  2:2467515/76‑1 (MQ=255)
                                                                        |                                                                      
GCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCA  >  NC_000913/1617351‑1617493

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