Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,076,192 G→A 29.0% G457G (GGG→GGA ydhS → conserved hypothetical protein with FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,076,1920GA29.0% 80.9 / 33.2 79G457G (GGG→GGAydhSconserved hypothetical protein with FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base G (23/33);  new base A (5/18);  total (28/51)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGC  >  minE/1076129‑1076259
                                                               |                                                                   
tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                                 >  1:1224745/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                                 >  1:2598357/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                                 >  1:2224657/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                                 >  1:1641089/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                                 >  1:1561542/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                                 >  1:911511/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                                 >  1:858842/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                                 >  1:2933413/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                                 >  1:2779412/1‑68 (MQ=255)
tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                                 >  1:2737625/1‑68 (MQ=255)
     tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:968432/69‑1 (MQ=255)
     tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:2736850/69‑1 (MQ=255)
     tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:1640326/69‑1 (MQ=255)
      ttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                                          <  1:1061679/69‑1 (MQ=255)
      ttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                                          <  1:2154327/69‑1 (MQ=255)
      ttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                                          <  1:472269/69‑1 (MQ=255)
       tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:1891152/67‑1 (MQ=255)
       tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:1203737/67‑1 (MQ=255)
       tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                                          >  1:2784898/1‑68 (MQ=255)
       tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                                          >  1:1060416/1‑68 (MQ=255)
       tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                                          >  1:1235240/1‑68 (MQ=255)
             gATGCCGGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:3266705/1‑68 (MQ=255)
             gATGCCCTCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:787229/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:2437516/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:2170766/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:3350062/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:33372/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:1078686/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                    >  1:3039718/1‑68 (MQ=255)
             gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCCAGCAATTACa                                                    >  1:1600531/1‑68 (MQ=255)
                   cgcgGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAAc                                              <  1:704549/68‑1 (MQ=255)
                   cgcgGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAAc                                              <  1:3102356/68‑1 (MQ=255)
                           aCGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGa                                     <  1:2285654/69‑1 (MQ=255)
                           aCGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGa                                     <  1:1813946/69‑1 (MQ=255)
                           aCGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGa                                     <  1:1963861/69‑1 (MQ=255)
                              tCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAAt                                  <  1:30388/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:2098149/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:296943/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:1737017/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:2797619/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:3248732/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:2557028/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:2182327/69‑1 (MQ=255)
                               ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:444412/69‑1 (MQ=255)
                                  cTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTAAAGAAAACAGGGGATGAAATc                                 <  1:2947486/66‑1 (MQ=255)
                                   ttaaagttaaagATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:2146518/39‑1 (MQ=255)
                                   ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:556286/39‑1 (MQ=255)
                                   ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:645017/39‑1 (MQ=255)
                                   ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:2209602/39‑1 (MQ=255)
                                   ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:3396276/39‑1 (MQ=255)
                                   ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:1777116/39‑1 (MQ=255)
                                   ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:1256279/39‑1 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGGTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa                            >  1:3371641/1‑69 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                                           <  1:3305092/38‑1 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa                            >  1:904231/1‑69 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa                            >  1:1048341/1‑69 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa                            >  1:1434324/1‑69 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa                            >  1:2720623/1‑69 (MQ=255)
                                    tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa                            >  1:2578371/1‑69 (MQ=255)
                                       aGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                                 <  1:2102108/61‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:434242/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:1006432/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:1677124/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:1827687/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:891151/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:1936989/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:850307/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:1943505/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:2120397/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:603438/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:2132796/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:410742/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:3435040/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:2531267/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:2705195/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:27734/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:3232017/69‑1 (MQ=255)
                                                tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                <  1:2837528/69‑1 (MQ=255)
                                                tatCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGggagga                                     <  1:423809/48‑1 (MQ=255)
                                                              ggCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGc  >  1:2349506/1‑69 (MQ=255)
                                                              ggCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGc  >  1:2150121/1‑69 (MQ=255)
                                                              ggCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGc  >  1:2914637/1‑69 (MQ=255)
                                                              ggCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGc  >  1:2538725/1‑69 (MQ=255)
                                                              ggCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGc  >  1:1212244/1‑69 (MQ=255)
                                                               |                                                                   
TGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGC  >  minE/1076129‑1076259

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: