Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,468 A→T 14.3% intergenic (+162/+110) 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,4680AT14.3% 169.4 / 18.7 73intergenic (+162/+110)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 (45/15);  new base T (6/3);  total (52/21)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAT  >  minE/2013405‑2013536
                                                               |                                                                    
tCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGa                                                               >  1:716606/1‑71 (MQ=255)
   aCATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGc                                                            <  1:482222/71‑1 (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:2858905/70‑1 (MQ=255)
                   tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                             >  1:2200219/1‑70 (MQ=255)
                   tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                             <  1:3774097/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: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: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:2867528/1‑71 (MQ=255)
                                   gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                            >  1:616682/1‑71 (MQ=255)
                                   gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                            >  1:263891/1‑71 (MQ=255)
                                   gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCa                            >  1:2823105/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: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)
                                                               |                                                                    
TCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAT  >  minE/2013405‑2013536

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: