Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,827,357 G→A 52.0% A298V (GCA→GTA)  dcuA ← C4‑dicarboxylate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,827,3570GA52.0% ‑2.1 / 110.5 75A298V (GCA→GTA) dcuAC4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base G (18/18);  new base A (14/25);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GTCGGCAGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACAGCAGAGCAGAAGCAAAGA  >  minE/2827288‑2827423
                                                                     |                                                                  
gTCGGCAGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGc                                                                   <  1:3940276/71‑1 (MQ=255)
  cGGCAGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGc                                                                   >  1:3525425/1‑69 (MQ=255)
  cGGCAGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGc                                                                   >  1:2561874/1‑69 (MQ=255)
     cAGAATGAACAGCCCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCagag                                                              >  1:49350/1‑71 (MQ=255)
      aGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGc                                                                   <  1:936149/65‑1 (MQ=255)
      aGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGc                                                                   <  1:3517165/65‑1 (MQ=255)
        aaTGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCa                                                           >  1:2013617/1‑71 (MQ=255)
          tGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATc                                                         <  1:3249699/71‑1 (MQ=255)
            aaCAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCAt                                                          >  1:1115424/1‑68 (MQ=255)
            aaCAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCgg                                                       >  1:2816025/1‑71 (MQ=255)
              cagaccagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCa                                                     >  1:477650/1‑71 (MQ=255)
               agaccagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCAt                                                    <  1:2382624/71‑1 (MQ=255)
                gacgagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCAt                                                    <  1:568435/70‑1 (MQ=255)
                gaccagacACCGCAGCGAAAGAAGCAACAGGGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATc                                                   <  1:2838293/71‑1 (MQ=255)
                gaccagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATc                                                   <  1:3505826/71‑1 (MQ=255)
                gaccagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCAt                                                    >  1:460566/1‑70 (MQ=255)
                gaccagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCAt                                                    >  1:465492/1‑70 (MQ=255)
                   cagacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGt                                                >  1:148141/1‑71 (MQ=255)
                     gacACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGc                                              <  1:3360491/71‑1 (MQ=255)
                      acacCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGc                                                             <  1:1789858/55‑1 (MQ=255)
                      acacCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCt                                             <  1:1265864/71‑1 (MQ=255)
                       cacCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCtt                                            <  1:1263857/71‑1 (MQ=255)
                        acCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCa                                                           >  1:817863/1‑55 (MQ=255)
                           tcAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCATTACCAGAGc                                                             <  1:3978942/49‑1 (MQ=38)
                           gCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCa                                                           >  1:27663/1‑52 (MQ=255)
                           gCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCa                                                  >  1:1664361/1‑61 (MQ=255)
                            cAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCg                                       <  1:760327/71‑1 (MQ=255)
                              gCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCtt                                            <  1:1086509/64‑1 (MQ=255)
                               cGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGtt                                    >  1:2285654/1‑71 (MQ=255)
                                gAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTg                                   <  1:1187818/71‑1 (MQ=255)
                                gAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTg                                   >  1:326627/1‑71 (MQ=255)
                                 aaaGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTg                                   >  1:660755/1‑70 (MQ=255)
                                 aaaGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATc                                                         >  1:436301/1‑48 (MQ=255)
                                 aaaGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTg                                   >  1:2639429/1‑70 (MQ=255)
                                 aaaGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGc                                  <  1:373301/71‑1 (MQ=255)
                                  aagaagCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCttt                                           >  1:1996985/1‑61 (MQ=255)
                                   agaagCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCa                                                     >  1:3530644/1‑50 (MQ=255)
                                   agaagCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAg                                <  1:2740174/71‑1 (MQ=255)
                                   agaagCAACAGCGGTCAGCGCTGAAACGTTCAGTGCCAGAGCCATCGGCa                                                     >  1:3525875/1‑50 (MQ=255)
                                    gaagCAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGc                               <  1:1418006/71‑1 (MQ=255)
                                       gCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTg                            <  1:2547253/71‑1 (MQ=255)
                                        cAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTga                           <  1:744151/71‑1 (MQ=255)
                                        cAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGt                                     <  1:1454339/61‑1 (MQ=255)
                                        cAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAg                                <  1:1754437/66‑1 (MQ=255)
                                        cAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAg                                <  1:3402266/66‑1 (MQ=255)
                                        cAACAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTga                           <  1:740465/71‑1 (MQ=255)
                                         aaCAGCGGTCAGCTGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTga                           <  1:834654/70‑1 (MQ=255)
                                          aCAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCAt                                                    <  1:1542013/44‑1 (MQ=255)
                                          aCAGCGGTCAGCGGTGAAACGTTCAGTACCAGAGCCATCGGCATCAGt                                                <  1:1304057/48‑1 (MQ=255)
                                           cagcggtcagcggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGagag                        <  1:54505/71‑1 (MQ=255)
                                           cagcggtcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGagag                        <  1:3355769/71‑1 (MQ=255)
                                           cagcggtcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGagag                        <  1:1412708/71‑1 (MQ=255)
                                           cagcggtcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGagag                        <  1:2322638/71‑1 (MQ=255)
                                            agcggtcagcggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGt                                                >  1:1778448/1‑46 (MQ=255)
                                            agcggtcagcggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGt                       <  1:1303808/71‑1 (MQ=255)
                                            agcggtcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGagag                        <  1:1358983/70‑1 (MQ=255)
                                             gcggtcagcgttGAAACGTTCAGTACCAGAGCCATCGGCATCa                                                  >  1:1542326/1‑43 (MQ=38)
                                             gcggtcagcggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTa                      >  1:1509770/1‑71 (MQ=255)
                                               ggtcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCttt                                           >  1:467041/1‑48 (MQ=255)
                                                 tcagcggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGc                                        <  1:4055964/49‑1 (MQ=255)
                                                 tcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCtt                                            <  1:3824400/45‑1 (MQ=255)
                                                 tcagcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTa                      >  1:1813675/1‑67 (MQ=255)
                                                   agcggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTga                           >  1:119673/1‑60 (MQ=255)
                                                    gcggtGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTg                                         <  1:1283334/45‑1 (MQ=255)
                                                      ggtGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACagcagag              <  1:3173733/70‑1 (MQ=255)
                                                        tGAAACGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACagcagagcag           <  1:3861806/71‑1 (MQ=255)
                                                         gAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCa                                 >  1:1617009/1‑48 (MQ=255)
                                                         gAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCa                                 >  1:3518380/1‑48 (MQ=255)
                                                          aaaCGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTg                            >  1:3466354/1‑52 (MQ=255)
                                                          aaaCGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACagcagagcag           >  1:2531554/1‑69 (MQ=255)
                                                            aCGTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCt                             <  1:3552552/49‑1 (MQ=255)
                                                              gTTCAGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCt                             <  1:3064720/47‑1 (MQ=255)
                                                                 cAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACAGCAGAGCAGAAGCaatg   >  1:1166505/1‑68 (MQ=255)
                                                                  aGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACAGCAGAGCAGAAGCAaaga  <  1:3419423/70‑1 (MQ=255)
                                                                  aGTACCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACAGCAGAGCAGAAGCAaaga  <  1:1404452/70‑1 (MQ=255)
                                                                     |                                                                  
GTCGGCAGAATGAACAGACCAGACACCGCAGCGAAAGAAGCAACAGCGGTCAGCGGTGAAACGTTCAGTGCCAGAGCCATCGGCATCAGTGCTTTTGCGGTTGCAGCCTGAGAGTACAGCAGAGCAGAAGCAAAGA  >  minE/2827288‑2827423

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: