Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,481,282 G→A 18.6% A362T (GCT→ACT)  glpA → sn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,481,2820GA18.6% 156.0 / 26.3 70A362T (GCT→ACT) glpAsn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (20/37);  new base A (4/9);  total (24/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TGCTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGG  >  minE/1481212‑1481351
                                                                      |                                                                     
tGCTTAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATgg                                                                       <  1:651783/71‑1 (MQ=255)
tGCTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATgg                                                                       >  1:1282994/1‑71 (MQ=255)
tGCTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATgg                                                                       >  1:2060961/1‑71 (MQ=255)
tGCTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATgg                                                                       <  1:1486816/71‑1 (MQ=255)
 gCTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGc                                                                      >  1:3354093/1‑71 (MQ=255)
  cTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCt                                                                     <  1:479811/71‑1 (MQ=255)
  cTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCt                                                                     <  1:3628309/71‑1 (MQ=255)
     aaCGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGa                                                                   <  1:1648144/70‑1 (MQ=255)
     aaCGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGa                                                                   <  1:2504484/70‑1 (MQ=255)
      aCGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAAt                                                                 >  1:1517170/1‑71 (MQ=255)
      aCGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAAt                                                                 <  1:3302272/71‑1 (MQ=255)
      aCGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGACTGAAt                                                                 <  1:235552/71‑1 (MQ=255)
       cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATgg                                                                       >  1:126977/1‑64 (MQ=255)
       cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAAt                                                                 <  1:1652288/70‑1 (MQ=255)
       cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATg                                                                <  1:978081/71‑1 (MQ=255)
       cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATg                                                                >  1:2022635/1‑71 (MQ=255)
       cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATg                                                                <  1:3316705/71‑1 (MQ=255)
       cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGACTGAAt                                                                 <  1:1911579/70‑1 (MQ=255)
        gcgATGGTCTGGACTGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGGTGGCTGAATgg                                                               <  1:2047663/71‑1 (MQ=255)
        gcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                               <  1:444290/71‑1 (MQ=255)
          gATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTTaa                                                                  <  1:1194643/66‑1 (MQ=255)
          gATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                               >  1:390884/1‑69 (MQ=255)
           aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATg                                                                <  1:3080013/67‑1 (MQ=255)
           aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATg                                                                <  1:2070762/67‑1 (MQ=255)
               tCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGATGATGGCTGAATg                                                                <  1:2301032/63‑1 (MQ=255)
                cTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTa                                                           >  1:3378366/1‑67 (MQ=255)
                   gACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGAcgc                                                    <  1:1279282/71‑1 (MQ=255)
                   gACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGAcg                                                     <  1:3624479/70‑1 (MQ=255)
                   gACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGAc                                                      <  1:786672/69‑1 (MQ=255)
                     cGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCgg                                                  <  1:3742607/71‑1 (MQ=255)
                       gATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGACTGAATGGGCTACCGACGCGGTa                                                <  1:3205296/71‑1 (MQ=255)
                        aTTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGACTGAATGGGCTACCGACGCGGTAt                                               <  1:2588667/71‑1 (MQ=255)
                          ttATCACCATCACCGGTGGCAATCTGATGACCTATCGGCTGATGGCTGAAt                                                                 <  1:2372487/51‑1 (MQ=255)
                            atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCg                                           >  1:1756713/1‑71 (MQ=255)
                             tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCg                                           <  1:3599200/70‑1 (MQ=255)
                               accatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCa                                         <  1:1672119/70‑1 (MQ=255)
                               accatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCa                                         >  1:2159596/1‑70 (MQ=255)
                                   tcaccGGTGTCAAACTGATGACCTATCGGCTGATGACTGAATGGGc                                                             <  1:387746/46‑1 (MQ=255)
                                   tcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTg                                    >  1:2284082/1‑71 (MQ=255)
                                   tcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTg                                    >  1:1296741/1‑71 (MQ=255)
                                   tcaccGGTGGCAAACTGATGACCTATCGGCTGATGACTGAATGGGCTAc                                                          >  1:2518953/1‑49 (MQ=255)
                                   tcaccGGTGGCAAACTGATGACCTATCGGCTGATGACTGAATGGGCTAc                                                          >  1:2080953/1‑49 (MQ=255)
                                     accGGTGGCAAACTGATGACCTATCGGCTGATGACTGAATGGGCTACCg                                                        >  1:710125/1‑49 (MQ=255)
                                      ccGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTggg                                  >  1:1606624/1‑70 (MQ=255)
                                       cGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCaaa                                       >  1:3017095/1‑64 (MQ=255)
                                       cGGTGGCAAACTGATGACCTATCGGCTGATGACTGAATGGTCTACCGACGCgg                                                  <  1:3453425/53‑1 (MQ=255)
                                         gTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCaaa                                       >  1:1974832/1‑62 (MQ=255)
                                         gTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCaa                               >  1:1970090/1‑70 (MQ=255)
                                          tGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAaca                             <  1:2494430/71‑1 (MQ=255)
                                            gCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCa                                         <  1:2274322/57‑1 (MQ=255)
                                             cAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACAcgc                          <  1:3259374/71‑1 (MQ=255)
                                             cAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACAcgc                          <  1:1699960/71‑1 (MQ=255)
                                              aaaCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAAc                                      <  1:3106214/58‑1 (MQ=255)
                                               aaCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCTGTATGCCGCAAACTgg                                   >  1:356323/1‑60 (MQ=255)
                                               aaCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAAc                                      <  1:1151073/57‑1 (MQ=255)
                                               aaCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                        <  1:1846483/71‑1 (MQ=255)
                                                aCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTAt                                               <  1:1792730/47‑1 (MQ=255)
                                                aCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGcc                                            >  1:2830702/1‑50 (MQ=255)
                                                aCTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGcc                       >  1:1419106/1‑71 (MQ=255)
                                                 cTGATGACCTATCGGCTGATGGCAGAATGGGCTACCGACGCGGTATGCCGCAAACTg                                    >  1:2871797/1‑57 (MQ=255)
                                                   gatgaCCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCaa                                        <  1:1903581/51‑1 (MQ=255)
                                                   gatgaCCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCaa                                        <  1:3454870/51‑1 (MQ=255)
                                                    atgaCCTATCGGCTGATGACTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTg                   <  1:757107/71‑1 (MQ=255)
                                                     tgaCCTATCGGCTGATGACTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTg                   <  1:752495/70‑1 (MQ=255)
                                                          tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGc                                 <  1:2901113/51‑1 (MQ=255)
                                                           aTCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACt            >  1:654000/1‑71 (MQ=255)
                                                            tCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACt            <  1:2610388/70‑1 (MQ=255)
                                                            tCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACt            <  1:78971/70‑1 (MQ=255)
                                                                cTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGa       <  1:3911351/71‑1 (MQ=255)
                                                                cTGATGACTGAATGGGCTACCGACGCGGTATGCCGCAAACTggg                                  >  1:729889/1‑44 (MQ=255)
                                                                 tGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTAc                >  1:3918869/1‑61 (MQ=255)
                                                                   aTGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATc     <  1:1821353/70‑1 (MQ=255)
                                                                   aTGACTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATc     <  1:1750305/70‑1 (MQ=255)
                                                                    tGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAacac                            <  1:3111914/46‑1 (MQ=255)
                                                                    tGGCTGAATGCGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTg   >  1:2234755/1‑71 (MQ=255)
                                                                      gCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCt                    >  1:3574455/1‑52 (MQ=255)
                                                                      gCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTgg  <  1:2399675/70‑1 (MQ=255)
                                                                      |                                                                     
TGCTGAACGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGG  >  minE/1481212‑1481351

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: