Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,466 Δ1 bp 18.1% intergenic (+160/+112) 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,4660T.18.1% 182.5 / 43.4 72intergenic (+160/+112)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 (44/15);  new base . (7/6);  total (51/21)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAT  >  minE/2013397‑2013536
                                                                     |                                                                      
aTGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCt                                                                        <  1:48596/70‑1 (MQ=255)
        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:3774097/70‑1 (MQ=255)
                           tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTa                                             <  1:2858905/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: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:917597/1‑70 (MQ=255)
                                     tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTAc                                   >  1:103821/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: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:911944/1‑70 (MQ=255)
                                              tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACg                          >  1:1517235/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)
                                                                     |                                                                      
ATGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAT  >  minE/2013397‑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: