Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,836 G→A 31.9% V219V (GTG→GTA ygjG → putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,070,8360GA31.9% 44.6 / 23.7 53V219V (GTG→GTAygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (5/31);  new base A (6/11);  total (11/42)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGA  >  minE/2070779‑2070898
                                                         |                                                              
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:84539/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:2737541/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:2913807/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:1777287/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:1649182/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:319116/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:58380/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:835846/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATGTACCGTTTGGc                                                       <  1:1188394/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGTCATCTACCGTTTGGc                                                       <  1:3149774/67‑1 (MQ=255)
ggCGAAATCGACCTTCCGCAAACCTTTTATGCCACTGCTGCCAGGGTTTCGCCATGTACCGTTTGGc                                                       <  1:835924/67‑1 (MQ=255)
 gCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAat                                                    <  1:2492461/69‑1 (MQ=255)
                           tATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAat                                                    <  1:3265012/43‑1 (MQ=255)
                              gCCGTTACTGCCGGGCTTCCGTCATGTGCCTTTTGGCAATATCGAAGc                                            <  1:428991/48‑1 (MQ=255)
                                    aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCa                  >  1:2227414/1‑68 (MQ=255)
                                    aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCa                  >  1:1424714/1‑68 (MQ=255)
                                    aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCa                  >  1:129595/1‑68 (MQ=255)
                                    aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCa                  >  1:1655725/1‑68 (MQ=255)
                                    aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCa                  >  1:3134690/1‑68 (MQ=255)
                                          gggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaaa             <  1:1215032/67‑1 (MQ=255)
                                          gggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaaa             <  1:2679208/67‑1 (MQ=255)
                                          gggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaaa             <  1:2758346/67‑1 (MQ=255)
                                          gggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCaaaaaa             <  1:3027795/67‑1 (MQ=255)
                                             cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACg                        <  1:2314750/53‑1 (MQ=255)
                                             cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACg                        <  1:479324/53‑1 (MQ=255)
                                             cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACg                        <  1:1060656/53‑1 (MQ=255)
                                             cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACg                        <  1:214828/53‑1 (MQ=255)
                                             cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACg                        <  1:1781301/53‑1 (MQ=255)
                                                tcGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGCTCTTAACGa                       >  1:2145945/2‑51 (MQ=255)
                                                tcGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGCTCTTAACGa                       >  1:1297477/2‑51 (MQ=255)
                                                tcGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGCTCTTAACGa                       >  1:2059532/2‑51 (MQ=255)
                                                tcGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGCTCTTAACGa                       >  1:2928934/2‑51 (MQ=255)
                                                tcGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGCTCTTAACGa                       >  1:409526/2‑51 (MQ=255)
                                                tcGTCATGTACCGTTTGGCAATATCGAAGCCATGCTCACGGCTCTTAACGa                       >  1:380957/2‑51 (MQ=255)
                                                  gTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgat    <  1:1373926/68‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1585615/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:664557/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1334454/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:529832/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1411514/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:411499/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1416802/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:3245557/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:2199247/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1595535/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:3072550/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1805537/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1999101/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:2145900/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1015234/69‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatg   <  1:3135181/68‑1 (MQ=255)
                                                   tCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatg   <  1:3131352/68‑1 (MQ=255)
                                                    cATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGtgatga  <  1:1277115/68‑1 (MQ=255)
                                                         |                                                              
GGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGA  >  minE/2070779‑2070898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: