Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 931,511 T→G 6.4% V235V (GTA→GTC yddG ← predicted methyl viologen efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE931,5110TG 6.4% 229.3 / 11.4 94V235V (GTA→GTCyddGpredicted methyl viologen efflux pump
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base G (4/2);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATG  >  minE/931442‑931578
                                                                     |                                                                   
cAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGt                                                                     >  1:1652314/1‑70 (MQ=255)
 aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    >  1:3185868/1‑70 (MQ=255)
 aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    >  1:4579910/1‑70 (MQ=255)
 aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    <  1:911771/70‑1 (MQ=255)
 aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    >  1:1093756/1‑70 (MQ=255)
 aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    <  1:466854/70‑1 (MQ=255)
    gcgcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    >  1:381432/1‑67 (MQ=255)
    gcgcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAgg                                                                >  1:2855293/1‑71 (MQ=255)
    gcgcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAgg                                                                >  1:3609003/1‑71 (MQ=255)
      gcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCg                                                              >  1:187950/1‑71 (MQ=255)
       ccTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCg                                                              <  1:4199814/70‑1 (MQ=255)
        cTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCGCGGCTGCAAGCGCTGAGGAAAGTACAGGCGt                                                             >  1:707070/1‑70 (MQ=255)
        cTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGt                                                                     <  1:2949006/62‑1 (MQ=255)
        cTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGt                                                                     <  1:3839302/62‑1 (MQ=255)
        cTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTa                                                            <  1:3738026/71‑1 (MQ=255)
            ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                         <  1:3316984/70‑1 (MQ=255)
            ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                         <  1:2548261/70‑1 (MQ=255)
            ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                         >  1:3878668/1‑70 (MQ=255)
            ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAAt                                                        <  1:3594256/71‑1 (MQ=255)
             cAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATa                                                       <  1:728425/71‑1 (MQ=255)
             cAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATa                                                       >  1:1803735/1‑71 (MQ=255)
             cAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATa                                                       <  1:308792/71‑1 (MQ=255)
              agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                                      >  1:195681/1‑71 (MQ=255)
              agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                                      >  1:3429517/1‑71 (MQ=255)
              agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                                      >  1:2294263/1‑71 (MQ=255)
                aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                                     >  1:4494013/1‑70 (MQ=255)
                aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                                     >  1:1308068/1‑70 (MQ=255)
                aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                                     >  1:2420407/1‑70 (MQ=255)
                aaagagaacgaCAGCGGGGCGCTGAGCAGCAGGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGa                                                    <  1:3502674/68‑1 (MQ=255)
                 acgagaacgaCAGCGGGTCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGa                                                    <  1:2310336/70‑1 (MQ=255)
                 acgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                                    <  1:3623272/54‑1 (MQ=255)
                 acgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGaa                                                   >  1:2541770/1‑71 (MQ=255)
                 acgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGa                                                    <  1:2315796/70‑1 (MQ=255)
                   gagaacgaCAGCGGGGCGCTGAGTAACACGGCTGCAAGTGCTGAGGAAAGGACGGGCGTaaaa                                                         <  1:2456821/63‑1 (MQ=255)
                   gagaacgaCAGCGGGGCGCTGAGTAACACGGCTGCAAGTGCTGAGGAAAGGACGGGCGTaaaa                                                         <  1:2003876/63‑1 (MQ=255)
                   gagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAAcc                                                 <  1:2072365/71‑1 (MQ=255)
                    agaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCt                                                <  1:1997467/71‑1 (MQ=255)
                     gaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCt                                                <  1:2762331/70‑1 (MQ=255)
                        cgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGaa                                                   >  1:1115834/1‑64 (MQ=255)
                         gaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaa                                                           <  1:3408388/55‑1 (MQ=255)
                             gCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                                     <  1:139828/57‑1 (MQ=255)
                              cGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATaa                                      <  1:299672/71‑1 (MQ=255)
                               ggggCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAAt                                     >  1:2723826/1‑71 (MQ=255)
                               ggggCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAAt                                     >  1:1854054/1‑71 (MQ=255)
                                    gcTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                                <  1:2848767/71‑1 (MQ=255)
                                       gAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                         <  1:1391681/43‑1 (MQ=255)
                                        agcagcACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                         >  1:4393610/1‑42 (MQ=255)
                                          cagcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaa                                                           <  1:2964033/38‑1 (MQ=255)
                                           agcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACAt                             <  1:1153927/67‑1 (MQ=255)
                                            gcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                                <  1:2570046/63‑1 (MQ=255)
                                            gcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                                <  1:3616308/63‑1 (MQ=255)
                                            gcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                                <  1:1078188/63‑1 (MQ=255)
                                              aCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGt                                  >  1:3383940/1‑59 (MQ=255)
                                              aCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGt                                  >  1:1039202/1‑59 (MQ=255)
                                              aCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCacg                      >  1:3632096/1‑69 (MQ=255)
                                              aCGCCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATg                      <  1:938935/71‑1 (MQ=255)
                                               cGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                         >  1:1552703/1‑35 (MQ=255)
                                               cGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                     <  1:3228335/71‑1 (MQ=255)
                                               cGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                     <  1:2423196/71‑1 (MQ=255)
                                               cGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                     <  1:4318387/71‑1 (MQ=255)
                                                ggCTGCAAGTGCTGAGGAAAGGACGGGCGTAAAATACGaa                                                   >  1:4248562/1‑40 (MQ=255)
                                                ggCTGCAAGTGCTGAGGAAAGGACGGGCGTAAAATACGaa                                                   >  1:219954/1‑40 (MQ=255)
                                                ggCTGCAAGTGCTGAGGAAAGGACGGGCGTAAAATACGaa                                                   >  1:177766/1‑40 (MQ=255)
                                                ggCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                     >  1:279164/1‑70 (MQ=255)
                                                ggCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                     >  1:2225048/1‑70 (MQ=255)
                                                ggCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                     >  1:705990/1‑70 (MQ=255)
                                                ggCTGCAAGAGCTGAGGAAAGGACGGGCGTAAAATACGaa                                                   >  1:4128742/1‑40 (MQ=255)
                                                 gCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGa                                                    >  1:4405791/1‑38 (MQ=255)
                                                 gCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTTACATTGCCATGCaa                   <  1:2762955/71‑1 (MQ=255)
                                                 gCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCaa                   <  1:3859284/71‑1 (MQ=255)
                                                 gCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCaa                   <  1:3796210/71‑1 (MQ=255)
                                                     cAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATaa                                      >  1:3377123/1‑48 (MQ=255)
                                                     cAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCaa                   <  1:2984253/67‑1 (MQ=255)
                                                     cAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAatat                <  1:279496/70‑1 (MQ=255)
                                                      aaGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCAt                                        <  1:3045214/45‑1 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTa                                               >  1:640460/1‑37 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTa                                               >  1:1029124/1‑37 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACAtt                            >  1:1049475/1‑56 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACAtt                            >  1:4270215/1‑56 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAata                 >  1:3148085/1‑67 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAata                 >  1:1725198/1‑67 (MQ=255)
                                                       aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATAcc             <  1:1945295/71‑1 (MQ=255)
                                                        gcgcTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATAcc             >  1:1259480/1‑70 (MQ=255)
                                                           cTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATaa                                      >  1:892426/1‑42 (MQ=255)
                                                           cTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGt                                  <  1:3030806/46‑1 (MQ=255)
                                                           cTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGt                                  <  1:4323627/46‑1 (MQ=255)
                                                           cTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACa         <  1:2007879/71‑1 (MQ=255)
                                                            tGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACAt        >  1:2932774/1‑71 (MQ=255)
                                                             gAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGt                                  <  1:3257032/44‑1 (MQ=255)
                                                              aGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCTACATTc      >  1:3932677/1‑71 (MQ=255)
                                                              aGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACAtt       <  1:1202858/70‑1 (MQ=255)
                                                              aGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTc      >  1:3770902/1‑71 (MQ=255)
                                                               ggAAAGTACAGGCGTAAAATACGAACCTACCGCCTTAATGGTGAc                               <  1:1134642/45‑1 (MQ=255)
                                                               ggAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGAc                               <  1:1036843/45‑1 (MQ=255)
                                                               ggAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTcc     <  1:30771/71‑1 (MQ=255)
                                                                gAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATg                                    >  1:3991591/1‑39 (MQ=255)
                                                                  aaGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATg  >  1:1446881/1‑71 (MQ=255)
                                                                     tACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCa    >  1:4061286/1‑66 (MQ=255)
                                                                     |                                                                   
CAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATG  >  minE/931442‑931578

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