Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,460 A→G 16.9% intergenic (+154/+118) 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,4600AG16.9% 187.9 / 27.4 77intergenic (+154/+118)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 A (46/18);  new base G (7/6);  total (53/24)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CAAACGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTG  >  minE/2013390‑2013528
                                                                      |                                                                    
cAAACGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGa                                                                      <  1:622639/71‑1 (MQ=255)
  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: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: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:616682/1‑71 (MQ=255)
                                                  gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                    >  1:2823105/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:1517235/1‑70 (MQ=255)
                                                     tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                  >  1:911944/1‑70 (MQ=255)
                                                     tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                  >  1:1324369/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:2836133/1‑70 (MQ=255)
                                                              aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCaa         >  1:1435046/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)
                                                                      |                                                                    
CAAACGCATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTG  >  minE/2013390‑2013528

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: