Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 931,499 C→T 5.4% A239A (GCG→GCA yddG ← predicted methyl viologen efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE931,4990CT 5.4% 291.1 / 4.2 93A239A (GCG→GCAyddGpredicted methyl viologen efflux pump
Reads supporting (aligned to +/- strand):  ref base C (42/45);  new base T (3/2);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

GCCGCAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACC  >  minE/931430‑931567
                                                                     |                                                                    
gccgcAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAgcg                                                                     >  1:4548283/1‑71 (MQ=255)
 ccgcAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAgcgc                                                                    >  1:3709472/1‑71 (MQ=255)
  cgcAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGTGGCGCTGAGCAGCACGGCTGCAAGCGCt                                                                   <  1:4088489/71‑1 (MQ=255)
  cgcAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCt                                                                   <  1:3046442/71‑1 (MQ=255)
  cgcAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCt                                                                   <  1:3331430/71‑1 (MQ=255)
   gcAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCt                                                                   >  1:2023815/1‑70 (MQ=255)
    cAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCt                                                                   >  1:1457810/1‑69 (MQ=255)
        ccATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGa                                                              <  1:929439/70‑1 (MQ=255)
           tCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAg                                                           <  1:4154415/70‑1 (MQ=255)
            cAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGt                                                          >  1:1652314/1‑70 (MQ=255)
             aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                         >  1:1093756/1‑70 (MQ=255)
             aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                         >  1:4579910/1‑70 (MQ=255)
             aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                         <  1:466854/70‑1 (MQ=255)
             aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                         >  1:3185868/1‑70 (MQ=255)
             aGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTa                                                         <  1:911771/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:3878668/1‑70 (MQ=255)
                        ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                              <  1:2548261/70‑1 (MQ=255)
                        ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                              <  1:3316984/70‑1 (MQ=255)
                        ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAAt                                             <  1:3594256/71‑1 (MQ=255)
                         cAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATa                                            >  1:1803735/1‑71 (MQ=255)
                         cAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATa                                            <  1:728425/71‑1 (MQ=255)
                         cAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATa                                            <  1:308792/71‑1 (MQ=255)
                          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGa                                                                 >  1:2484692/1‑49 (MQ=255)
                          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                           >  1:2294263/1‑71 (MQ=255)
                          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                           >  1:195681/1‑71 (MQ=255)
                          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                           >  1:3429517/1‑71 (MQ=255)
                            aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                          >  1:4494013/1‑70 (MQ=255)
                            aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                          >  1:2420407/1‑70 (MQ=255)
                            aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                          >  1:1308068/1‑70 (MQ=255)
                            aaagagaacgaCAGCGGGGCGCTGAGCAGCAGGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGa                                         <  1:3502674/68‑1 (MQ=255)
                             acgagaacgaCAGCGGGTCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGa                                         <  1:2310336/70‑1 (MQ=255)
                             acgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAg                                                                >  1:296966/1‑47 (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)
                                 gaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGaaa                                                            >  1:2585757/1‑47 (MQ=255)
                                 gaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCt                                     <  1:2762331/70‑1 (MQ=255)
                                   acgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGa                                                              <  1:1026630/43‑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:1078188/63‑1 (MQ=255)
                                                        gcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                     <  1:3616308/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:177766/1‑40 (MQ=255)
                                                            ggCTGCAAGTGCTGAGGAAAGGACGGGCGTAAAATACGaa                                        >  1:219954/1‑40 (MQ=255)
                                                            ggCTGCAAGTGCTGAGGAAAGGACGGGCGTAAAATACGaa                                        >  1:4248562/1‑40 (MQ=255)
                                                            ggCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc          >  1:705990/1‑70 (MQ=255)
                                                            ggCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc          >  1:279164/1‑70 (MQ=255)
                                                            ggCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc          >  1:2225048/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:4270215/1‑56 (MQ=255)
                                                                   aGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACAtt                 >  1:1049475/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)
                                                                     |                                                                    
GCCGCAGACCATCAGCGCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACC  >  minE/931430‑931567

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: