Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,587,031 A→G 24.9% S82S (AGT→AGC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,587,0310AG24.9% 188.1 / 59.2 156S82S (AGT→AGCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (30/87);  new base G (20/19);  total (50/106)
Fisher's exact test for biased strand distribution p-value = 5.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 158 total aligned reads displayed.

CGATATCAGGAGATGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCACCTGAGTGACCGA  >  minE/2586970‑2587099
                                                             |                                                                    
         gagaTGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:3091380/68‑1 (MQ=255)
         gagaTGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:2320705/68‑1 (MQ=255)
         gagaTGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:895570/68‑1 (MQ=255)
          agaTGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:1053603/69‑1 (MQ=255)
           gaTGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTAc                                                    >  1:1583862/1‑69 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:277079/68‑1 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:3037747/68‑1 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:298260/68‑1 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:2120225/68‑1 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:871118/68‑1 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:91712/68‑1 (MQ=255)
                 ccAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAg                                               <  1:1817259/68‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGTTAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:3383741/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGTCATGTa                                                     <  1:436306/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:339990/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:2689060/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:3068563/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:1852549/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:2521624/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:382441/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:2072937/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:457941/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:718097/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:922226/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTa                                                     <  1:2263492/60‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGc                                             <  1:94298/68‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGc                                             <  1:2108714/68‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGc                                             <  1:2568968/68‑1 (MQ=255)
                   aGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGc                                             <  1:682374/68‑1 (MQ=255)
                      ctcGAAGGTCAGTGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:2152677/55‑1 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       >  1:2637518/1‑55 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:2714482/55‑1 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:3097449/55‑1 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:3368878/55‑1 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       >  1:403965/1‑55 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       >  1:647017/1‑55 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       >  1:2276859/1‑55 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       >  1:101612/1‑55 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:1014044/55‑1 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       >  1:1370278/1‑55 (MQ=255)
                      ctcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATg                                                       <  1:1128787/55‑1 (MQ=255)
                       tcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCa                                        >  1:2020164/1‑69 (MQ=255)
                       tcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCa                                        >  1:2773887/1‑69 (MQ=255)
                       tcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCa                                        >  1:1559539/1‑69 (MQ=255)
                       tcGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCa                                        >  1:79716/1‑69 (MQ=255)
                              tCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGt                                                                 <  1:2706070/37‑1 (MQ=255)
                                agagTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTc                               <  1:628735/69‑1 (MQ=255)
                                agagTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTc                               <  1:2579829/69‑1 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:1278252/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:543466/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:524659/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:2134557/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:1441001/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:3332510/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:3088369/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:1659866/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:2947261/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:1056403/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:2243889/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAttt                         >  1:1996907/1‑65 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAtt                          >  1:291385/1‑64 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATAtt                          >  1:1540385/1‑64 (MQ=255)
                                          agGCATTACCTGCCGCATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGACCATTCATAttt                         >  1:484122/1‑65 (MQ=38)
                                           ggCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTc                    <  1:2134307/69‑1 (MQ=255)
                                           ggCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTc                    <  1:2282790/69‑1 (MQ=255)
                                           ggCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTc                    <  1:1171372/69‑1 (MQ=255)
                                           ggCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTc                    <  1:3256954/69‑1 (MQ=255)
                                           ggCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTc                    <  1:257788/69‑1 (MQ=255)
                                           ggCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTc                    <  1:1300861/69‑1 (MQ=255)
                                            gCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCTATCAAGCCCATTCATATTTTGCTCa                   <  1:574417/69‑1 (MQ=255)
                                            gCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCa                   <  1:2737655/69‑1 (MQ=255)
                                            gCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCa                   <  1:3139567/69‑1 (MQ=255)
                                               ttGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAgac                                    <  1:2412770/49‑3 (MQ=255)
                                               ttGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGcc                                    <  1:1729256/49‑1 (MQ=255)
                                               ttGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGcc                                    <  1:899360/49‑1 (MQ=255)
                                               ttGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGcc                                    <  1:1210016/49‑1 (MQ=255)
                                               ttGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGcc                                    <  1:545236/49‑1 (MQ=255)
                                               ttGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGcc                                    <  1:3246802/49‑1 (MQ=255)
                                                 gCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCa               >  1:526911/1‑68 (MQ=255)
                                                 gCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCAc              >  1:232936/1‑69 (MQ=255)
                                                 gCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCAc              >  1:1502761/1‑69 (MQ=255)
                                                 gCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCAc              >  1:1233101/1‑69 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:3140211/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1777399/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1789964/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1123396/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1808537/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1905538/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:2936883/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:72746/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1519354/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:3414748/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:2425607/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:2489664/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1471923/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1925533/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1322604/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:1333064/54‑1 (MQ=255)
                                                       cgcATCGCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTg                       <  1:2928648/54‑1 (MQ=255)
                                                             aCTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCACCTGAGTGAcc    >  1:2370103/1‑67 (MQ=255)
                                                             |                                                                    
CGATATCAGGAGATGTCCCAGTCTCGAAGGTCAGAGTGATAGAGGCATTGCCCGCCGCATCACTGGTTGAAGACATGTACATCAGGCCATCAAGCCCATTCATATTTTGCTCAATCACCTGAGTGACCGA  >  minE/2586970‑2587099

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