Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,480 T→G 17.1% intergenic (+174/+98) yghA → / ← exbD predicted glutathionylspermidine synthase, with NAD(P)‑binding Rossmann‑fold domain/membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,013,4800TG17.1% 164.2 / 34.4 76intergenic (+174/+98)yghA/exbDpredicted glutathionylspermidine synthase, with NAD(P)‑binding Rossmann‑fold domain/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base T (46/17);  new base G (7/6);  total (53/23)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGAT  >  minE/2013412‑2013550
                                                                    |                                                                      
gTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGcgtcg                                                                      >  1:983786/1‑71 (MQ=255)
 ttGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGcgtcgt                                                                     <  1:3222267/71‑1 (MQ=255)
  tGCAATAGATTTTGTAGGCCTGATAAGCGTAGAGCATCAGGCATTGAGCACCTAATGCCGGATGcgtcgt                                                                     >  1:172582/1‑70 (MQ=255)
   gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTg                                                                    <  1:90943/70‑1 (MQ=255)
    cAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGa                                                                   <  1:3608013/70‑1 (MQ=255)
      aTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAAc                                                                 <  1:1029503/70‑1 (MQ=255)
          aTTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCtt                                                            <  1:1374020/71‑1 (MQ=255)
           ttttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                                           <  1:3422601/71‑1 (MQ=255)
            tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                                           <  1:2858905/70‑1 (MQ=255)
            tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                                           <  1:3774097/70‑1 (MQ=255)
            tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                                           >  1:2200219/1‑70 (MQ=255)
            tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTAt                                                          >  1:3264642/1‑71 (MQ=255)
             ctGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATcc                                                        >  1:3593501/2‑71 (MQ=14)
              tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACg                                                                >  1:1459170/1‑63 (MQ=255)
              tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATc                                                         <  1:408657/70‑1 (MQ=255)
               gTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCg                                                       >  1:83353/1‑71 (MQ=255)
                tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATc                                                         <  1:2163813/67‑1 (MQ=14)
                tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATc                                                         <  1:1368354/67‑1 (MQ=14)
                tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATcc                                                        <  1:2892251/68‑1 (MQ=14)
                tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATc                                                         <  1:3696935/67‑1 (MQ=14)
                tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATc                                                         <  1:1182418/67‑1 (MQ=14)
                 aGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATc                                                         >  1:309970/1‑66 (MQ=14)
                 aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGcgtc                                                                       >  1:353167/1‑53 (MQ=255)
                 aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGAACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGc                                                     >  1:517085/1‑71 (MQ=255)
                  ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                    >  1:2894970/1‑70 (MQ=14)
                  ggCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGc                                                     >  1:3596814/1‑70 (MQ=255)
                   gCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGcc                                                    >  1:1356117/1‑70 (MQ=255)
                    ccTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                  >  1:1120287/1‑70 (MQ=14)
                    ccTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                 <  1:3237446/71‑1 (MQ=14)
                     cTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                 >  1:2008421/1‑70 (MQ=11)
                     cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTa                                                  >  1:2139927/1‑70 (MQ=255)
                      tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTAc                                                 >  1:917597/1‑70 (MQ=255)
                      tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTAc                                                 >  1:103821/1‑70 (MQ=255)
                       gATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTAcgag                                              >  1:1688264/1‑68 (MQ=2)
                       gATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACa                                                >  1:1284305/1‑70 (MQ=255)
                       gATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACa                                                >  1:2907219/1‑70 (MQ=255)
                        aTAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCtt                                                            >  1:1902853/1‑56 (MQ=2)
                        aTAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACaaa                                              >  1:2412091/1‑71 (MQ=255)
                         tAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAAc                                             >  1:1491971/1‑71 (MQ=255)
                           aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACaaa                                              >  1:2414028/1‑68 (MQ=255)
                           aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACg                                            >  1:429728/1‑70 (MQ=255)
                            gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                                          >  1:263891/1‑71 (MQ=255)
                            gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                                          >  1:2823105/1‑71 (MQ=255)
                            gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                                          >  1:2867528/1‑71 (MQ=255)
                            gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                                          >  1:616682/1‑71 (MQ=255)
                             cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCg                                                       <  1:2575131/57‑1 (MQ=255)
                             cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                                          >  1:884712/1‑70 (MQ=255)
                              gTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCAc                                         >  1:3109714/1‑70 (MQ=255)
                               tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                                        >  1:1517235/1‑70 (MQ=255)
                               tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                                        >  1:1324369/1‑70 (MQ=255)
                               tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                                        >  1:911944/1‑70 (MQ=255)
                               tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGa                                       <  1:2053553/71‑1 (MQ=255)
                                  cgcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATc                                                         >  1:3552393/1‑50 (MQ=255)
                                   gcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTa                                                  >  1:36356/1‑56 (MQ=255)
                                   gcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGAt                                   >  1:2871749/1‑71 (MQ=255)
                                      tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTc                                 >  1:1952030/1‑70 (MQ=255)
                                      tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCa                                <  1:2082838/71‑1 (MQ=255)
                                        aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCaa                               >  1:1435046/1‑70 (MQ=255)
                                        aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCaa                               >  1:3271401/1‑70 (MQ=255)
                                        aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCaa                               >  1:2836133/1‑70 (MQ=255)
                                            aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACACACGCa                                          >  1:692061/1‑55 (MQ=255)
                                            aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGt                          >  1:2046440/1‑71 (MQ=255)
                                             ttGAGCACCTAATGCCGGATGCGTCGTGAACGCCtt                                                            <  1:2279280/36‑1 (MQ=255)
                                              tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTg                        >  1:1093660/1‑71 (MQ=255)
                                              tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAAAGCACGAAGATTCAAc                              >  1:1468535/1‑65 (MQ=255)
                                                 gCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCaa                     >  1:506707/1‑71 (MQ=255)
                                                  cACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCaa                     >  1:868365/1‑70 (MQ=255)
                                                   aCCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCaa                     >  1:1642281/1‑69 (MQ=255)
                                                      tAATGACGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAt                >  1:3340311/1‑71 (MQ=255)
                                                        aTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGa                                    >  1:3070173/1‑49 (MQ=255)
                                                         tGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAtttt             >  1:3020079/1‑71 (MQ=255)
                                                         tGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAtttt             >  1:3117297/1‑71 (MQ=255)
                                                           ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTg            <  1:2283825/70‑1 (MQ=255)
                                                           ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTg            <  1:1134529/70‑1 (MQ=255)
                                                           ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGt           <  1:914513/71‑1 (MQ=255)
                                                            cGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTg            >  1:1243747/1‑69 (MQ=255)
                                                              gATGCGTCGTTAACGCCTTATCCGGCCTACAAACTCACGAAGATTCAACATGTTGCAATAGATTTTGTAgg        <  1:570896/71‑1 (MQ=255)
                                                               aTGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGc       >  1:2857214/1‑71 (MQ=255)
                                                               aTGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGc       <  1:1795362/71‑1 (MQ=255)
                                                                tGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGcc      >  1:1204901/1‑71 (MQ=255)
                                                                  cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAg                                     <  1:2476160/38‑1 (MQ=255)
                                                                    tcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGAt                                   >  1:2930259/1‑38 (MQ=255)
                                                                    tcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCtgat  <  1:3010807/71‑1 (MQ=255)
                                                                    |                                                                      
GTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGAT  >  minE/2013412‑2013550

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: