Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,461 G→T 12.3% intergenic (+155/+117) 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,4610GT12.3% 209.3 / 11.8 76intergenic (+155/+117)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 G (45/18);  new base T (6/3);  total (52/24)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

AACGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAA  >  minE/2013392‑2013531
                                                                     |                                                                      
aaCGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGc                                                                       <  1:2578610/71‑1 (MQ=255)
 aCGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAg                                                                        >  1:871683/1‑69 (MQ=255)
 aCGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCa                                                                      >  1:1520295/1‑71 (MQ=255)
     aTGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCt                                                                   <  1:48596/70‑1 (MQ=255)
             tCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGa                                                          >  1:716606/1‑71 (MQ=255)
              cAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGc                                                                       <  1:3698880/57‑1 (MQ=255)
                aCATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCAc                                                                     >  1:3058957/1‑57 (MQ=255)
                aCATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCAc                                                                     >  1:3752914/1‑57 (MQ=255)
                aCATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGc                                                       <  1:482222/71‑1 (MQ=255)
                aCATGTAGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCAc                                                                     >  1:510251/1‑57 (MQ=255)
                 cATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGcg                                                      <  1:1552784/71‑1 (MQ=255)
                  aTGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGa                                                          >  1:343130/1‑66 (MQ=255)
                    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:3774097/70‑1 (MQ=255)
                                tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                        >  1:2200219/1‑70 (MQ=255)
                                tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                        <  1:2858905/70‑1 (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:1368354/67‑1 (MQ=14)
                                    tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATc                                      <  1:2163813/67‑1 (MQ=14)
                                    tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATcc                                     <  1:2892251/68‑1 (MQ=14)
                                    tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATc                                      <  1:1182418/67‑1 (MQ=14)
                                    tAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATc                                      <  1:3696935/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:103821/1‑70 (MQ=255)
                                          tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTAc                              >  1:917597/1‑70 (MQ=255)
                                           gATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTAcgag                           >  1:1688264/1‑68 (MQ=2)
                                           gATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACa                             >  1:2907219/1‑70 (MQ=255)
                                           gATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACa                             >  1:1284305/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:616682/1‑71 (MQ=255)
                                                gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                       >  1:2867528/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:1324369/1‑70 (MQ=255)
                                                   tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                     >  1:1517235/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:2836133/1‑70 (MQ=255)
                                                            aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCaa            >  1:3271401/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)
                                                                     |                                                                      
AACGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAA  >  minE/2013392‑2013531

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: