Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,558 A→C 8.1% V105V (GTT→GTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5580AC 8.1% 289.4 / 22.1 111V105V (GTT→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (63/39);  new base C (3/6);  total (66/45)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 114 total aligned reads displayed.

CGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312488‑312627
                                                                      |                                                                     
cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                                       <  1:3679459/71‑1 (MQ=255)
      ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                 >  1:3584377/1‑71 (MQ=255)
      ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                 >  1:2161056/1‑71 (MQ=255)
         cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               >  1:144054/1‑70 (MQ=255)
            gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACATGCGCGata                                                           >  1:4121072/1‑71 (MQ=255)
               gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATcc                                                        >  1:3016428/1‑71 (MQ=255)
                ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:3629191/1‑57 (MQ=255)
                ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                            <  1:882196/66‑1 (MQ=255)
                ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                       >  1:2171676/1‑71 (MQ=255)
                  gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                  <  1:510482/58‑1 (MQ=255)
                   aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                      <  1:1552607/69‑1 (MQ=255)
                   aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    <  1:3180559/71‑1 (MQ=255)
                    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:4026884/1‑70 (MQ=255)
                     cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                            >  1:1338591/1‑61 (MQ=255)
                     cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                            >  1:1272280/1‑61 (MQ=255)
                     cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                            >  1:2182012/1‑61 (MQ=255)
                        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:1240221/1‑49 (MQ=255)
                        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:1235457/1‑49 (MQ=255)
                        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                               >  1:92574/1‑71 (MQ=255)
                           gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:3028572/1‑69 (MQ=255)
                           gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:1210264/1‑69 (MQ=255)
                           gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:2360078/1‑69 (MQ=255)
                            cGGCATCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            <  1:737645/70‑2 (MQ=12)
                                 gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                       >  1:3535914/1‑69 (MQ=12)
                                 gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                       >  1:2110872/1‑69 (MQ=12)
                                 gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                       >  1:3531102/1‑69 (MQ=12)
                                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                                           <  1:1416572/50‑1 (MQ=255)
                                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:3570354/71‑1 (MQ=255)
                                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:2788957/71‑1 (MQ=255)
                                  caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                     <  1:3336376/71‑1 (MQ=255)
                                    aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:2351372/1‑71 (MQ=255)
                                    aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:2376829/1‑71 (MQ=255)
                                     cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTTAACCTGCGCGATATCCGCAt                                                    <  1:791388/53‑1 (MQ=255)
                                      gggaTCGACAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:3347239/68‑1 (MQ=255)
                                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                          <  1:2990753/62‑1 (MQ=255)
                                        cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 >  1:2012322/1‑69 (MQ=255)
                                        cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 >  1:3728487/1‑69 (MQ=255)
                                         aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              <  1:841145/71‑1 (MQ=255)
                                         aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              <  1:1578154/71‑1 (MQ=255)
                                         aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              <  1:2178774/71‑1 (MQ=255)
                                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                            >  1:1257401/1‑70 (MQ=255)
                                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                             <  1:3389683/69‑1 (MQ=255)
                                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                           >  1:3519374/1‑71 (MQ=255)
                                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                           >  1:3393173/1‑71 (MQ=255)
                                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                           >  1:1535932/1‑71 (MQ=255)
                                            gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                           >  1:3957681/1‑71 (MQ=255)
                                             ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:3038001/71‑1 (MQ=255)
                                             ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:1073045/71‑1 (MQ=255)
                                              cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                                          >  1:1442397/1‑38 (MQ=255)
                                              cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                         >  1:2646272/1‑71 (MQ=255)
                                              cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                         >  1:3511589/1‑71 (MQ=255)
                                              aaGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                         >  1:3665785/2‑71 (MQ=255)
                                               aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        >  1:2523077/1‑71 (MQ=255)
                                               aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAGAGGTCAGGGTGAt                        >  1:3466175/1‑71 (MQ=255)
                                               aGCTGCAGTTTGTTATGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        >  1:3491305/1‑71 (MQ=255)
                                                 cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                          <  1:734976/51‑1 (MQ=12)
                                                 cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                          <  1:1887770/51‑1 (MQ=12)
                                                 cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                          <  1:1655155/51‑1 (MQ=12)
                                                 cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                          <  1:3258854/51‑1 (MQ=12)
                                                 cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:1203974/67‑1 (MQ=255)
                                                 cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        >  1:3189015/1‑69 (MQ=255)
                                                 cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAAAGGTCAGGGTGATCt                      >  1:2820966/1‑71 (MQ=255)
                                                 cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        >  1:318973/1‑69 (MQ=255)
                                                  tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAtg                                                   <  1:3866301/41‑2 (MQ=255)
                                                  tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                     >  1:3588475/1‑71 (MQ=255)
                                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              <  1:3910271/60‑1 (MQ=255)
                                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                     <  1:2009352/69‑1 (MQ=255)
                                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                     <  1:3822765/69‑1 (MQ=255)
                                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                     <  1:3802011/69‑1 (MQ=255)
                                                      gTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACg                 >  1:3967590/1‑71 (MQ=255)
                                                      gTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACg                 >  1:1432303/1‑71 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                  <  1:2998188/53‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:1354870/54‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:2936538/54‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:633519/54‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:1632188/54‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:978370/61‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                <  1:3506615/71‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:3755505/61‑1 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                 >  1:1345308/1‑36 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:509066/1‑50 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                               >  1:3556662/1‑54 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        >  1:3636047/1‑61 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        >  1:674896/1‑61 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                   <  1:625446/66‑1 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:2535739/1‑71 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:1927374/1‑71 (MQ=255)
                                                          gTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                   <  1:1962980/65‑1 (MQ=255)
                                                           ttCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAcc            >  1:3602160/1‑71 (MQ=255)
                                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg          >  1:945818/1‑71 (MQ=255)
                                                              tGCACCTGAAGCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt         >  1:1220312/1‑71 (MQ=255)
                                                               gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        >  1:2641954/1‑71 (MQ=255)
                                                                cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           >  1:4083189/1‑67 (MQ=255)
                                                                   cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAg     >  1:1042034/1‑70 (MQ=255)
                                                                   cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt    >  1:1779245/1‑71 (MQ=255)
                                                                     gAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                     >  1:138925/1‑36 (MQ=255)
                                                                     gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa  >  1:3037525/1‑71 (MQ=255)
                                                                     gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa  >  1:2510550/1‑71 (MQ=255)
                                                                      aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        <  1:922801/64‑1 (MQ=255)
                                                                      aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc   <  1:4016983/69‑1 (MQ=255)
                                                                      |                                                                     
CGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312488‑312627

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