Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,575 C→G 6.7% A100P (GCG→CCG) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5750CG 6.7% 378.4 / 10.3 119A100P (GCG→CCG) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (59/50);  new base G (3/5);  total (64/55)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

AACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAC  >  minE/312507‑312645
                                                                    |                                                                      
aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                                        <  1:1552607/69‑1 (MQ=255)
aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                      <  1:3180559/71‑1 (MQ=255)
 aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                      >  1:4026884/1‑70 (MQ=255)
    tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                   >  1:35390/1‑70 (MQ=255)
     cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                                 >  1:92574/1‑71 (MQ=255)
        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                >  1:2360078/1‑69 (MQ=255)
         cGGCATCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                                              <  1:737645/70‑2 (MQ=12)
         cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                              <  1:1260947/70‑1 (MQ=255)
           gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                           <  1:304415/71‑1 (MQ=255)
              gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                         >  1:2110872/1‑69 (MQ=12)
              gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                         >  1:3535914/1‑69 (MQ=12)
              gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                         >  1:3531102/1‑69 (MQ=12)
              gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                        <  1:2788957/71‑1 (MQ=255)
               caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                       <  1:3336376/71‑1 (MQ=255)
                 aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                     >  1:2376829/1‑71 (MQ=255)
                   gggaTCGACAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                   <  1:3347239/68‑1 (MQ=255)
                   ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                                            <  1:2990753/62‑1 (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:3445840/1‑71 (MQ=255)
                         gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                             >  1:3519374/1‑71 (MQ=255)
                         gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                             >  1:1535932/1‑71 (MQ=255)
                         gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                             >  1:3957681/1‑71 (MQ=255)
                         gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                             >  1:2122566/1‑71 (MQ=255)
                          ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                            <  1:3038001/71‑1 (MQ=255)
                          ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                            <  1:1073045/71‑1 (MQ=255)
                           cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                           >  1:2646272/1‑71 (MQ=255)
                           aaGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                           >  1:3665785/2‑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: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)
                               tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAtg                                                                     <  1:3866301/41‑2 (MQ=255)
                               tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                       >  1:3588475/1‑71 (MQ=255)
                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                <  1:3910271/60‑1 (MQ=255)
                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                       <  1:3802011/69‑1 (MQ=255)
                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                       <  1:2009352/69‑1 (MQ=255)
                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                       <  1:3822765/69‑1 (MQ=255)
                                  aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                    <  1:2507194/71‑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:2936538/54‑1 (MQ=255)
                                    tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                   <  1:1354870/54‑1 (MQ=255)
                                    tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                   <  1:1632188/54‑1 (MQ=255)
                                    tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                   <  1:633519/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)
                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                >  1:1927374/1‑71 (MQ=255)
                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                >  1:2535739/1‑71 (MQ=255)
                                       gTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                     <  1:1962980/65‑1 (MQ=255)
                                        ttCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAcc                              >  1:3602160/1‑71 (MQ=255)
                                          cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                        >  1:507193/1‑43 (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:1760959/1‑71 (MQ=255)
                                                cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                      >  1:1779245/1‑71 (MQ=255)
                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGGACCCGTGGAGt                      >  1:765945/1‑70 (MQ=255)
                                                  gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                    >  1:3037525/1‑71 (MQ=255)
                                                  gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                    >  1:2510550/1‑71 (MQ=255)
                                                      cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                >  1:1798092/1‑71 (MQ=255)
                                                      cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                <  1:2621067/71‑1 (MQ=255)
                                                       tGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCAc                   <  1:2003998/67‑1 (MQ=255)
                                                        gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                     >  1:1874241/1‑48 (MQ=255)
                                                        gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                     >  1:4011056/1‑48 (MQ=255)
                                                        gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                >  1:813574/1‑69 (MQ=255)
                                                         cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGATTCACTGTTag             <  1:2926808/71‑1 (MQ=255)
                                                         cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa              <  1:1194559/70‑1 (MQ=255)
                                                         cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa              <  1:234697/70‑1 (MQ=255)
                                                         cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa              <  1:932816/70‑1 (MQ=255)
                                                         cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGGGATCTGCACGGTACCCGTGGAGTCACTGTTag             <  1:2374867/71‑1 (MQ=255)
                                                           cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                     <  1:1127175/61‑1 (MQ=255)
                                                             atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag             <  1:2601741/67‑1 (MQ=255)
                                                              tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                                  <  1:2038381/45‑1 (MQ=255)
                                                               atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                            >  1:3374315/1‑50 (MQ=255)
                                                               atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag             <  1:299290/65‑1 (MQ=255)
                                                                tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                          >  1:2764503/1‑35 (MQ=255)
                                                                tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                          >  1:90525/1‑35 (MQ=255)
                                                                tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGa                        <  1:842902/53‑1 (MQ=255)
                                                                tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa       <  1:2803921/70‑1 (MQ=255)
                                                                  cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg     >  1:598377/1‑70 (MQ=255)
                                                                  cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg     <  1:894640/70‑1 (MQ=255)
                                                                  cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg     >  1:3645555/1‑70 (MQ=255)
                                                                   gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa   <  1:1712299/71‑1 (MQ=255)
                                                                    cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc  <  1:4011723/71‑1 (MQ=255)
                                                                    |                                                                      
AACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAC  >  minE/312507‑312645

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: