Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 931,514 A→G 6.7% P234P (CCT→CCC yddG ← predicted methyl viologen efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE931,5140AG 6.7% 242.3 / 7.7 89P234P (CCT→CCCyddGpredicted methyl viologen efflux pump
Reads supporting (aligned to +/- strand):  ref base A (41/42);  new base G (4/2);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATGCAGCA  >  minE/931446‑931583
                                                                    |                                                                     
gcgcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAgg                                                                     >  1:3609003/1‑71 (MQ=255)
gcgcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAgg                                                                     >  1:2855293/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)
    cTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTa                                                                 <  1:3738026/71‑1 (MQ=255)
        ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                              <  1:3316984/70‑1 (MQ=255)
        ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                              >  1:3878668/1‑70 (MQ=255)
        ccAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTaaaa                                                              <  1:2548261/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)
          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                                           >  1:2294263/1‑71 (MQ=255)
          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                                           >  1:195681/1‑71 (MQ=255)
          agaacgagaacgACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATAc                                                           >  1:3429517/1‑71 (MQ=255)
            aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                                          >  1:2420407/1‑70 (MQ=255)
            aacgagaacgaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACg                                                          >  1:4494013/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)
             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:1078188/63‑1 (MQ=255)
                                        gcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                                     <  1:3616308/63‑1 (MQ=255)
                                        gcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGa                                     <  1:2570046/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:4318387/71‑1 (MQ=255)
                                           cGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                          <  1:3228335/71‑1 (MQ=255)
                                           cGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGc                          <  1:2423196/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: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:3796210/71‑1 (MQ=255)
                                             gCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCaa                        <  1:3859284/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)
                                                                  aCAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATgcag    >  1:2081873/1‑70 (MQ=255)
                                                                    aGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATgcagca  >  1:2872592/1‑70 (MQ=255)
                                                                    |                                                                     
GCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATGCAGCA  >  minE/931446‑931583

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: