Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,076,180 C→T 29.0% I453I (ATC→ATT 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,1800CT29.0% 81.9 / 32.8 79I453I (ATC→ATTydhSconserved hypothetical protein with FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (23/33);  new base T (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 = 6.39e-01

GACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGA  >  minE/1076112‑1076245
                                                                    |                                                                 
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:206815/69‑1 (MQ=255)
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:878741/69‑1 (MQ=255)
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:2513551/69‑1 (MQ=255)
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:740025/69‑1 (MQ=255)
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:161529/69‑1 (MQ=255)
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:2688827/69‑1 (MQ=255)
gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:15467/69‑1 (MQ=255)
                 tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                   >  1:2224657/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                   >  1:1641089/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                   >  1:2598357/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                   >  1:1224745/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGTGGGCAGCGTCCGCTTAAAGTGAAAGATATTCCTTTCCCTGGACTAc                                                   >  1:1561542/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:1001054/52‑1 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATc                                                                   <  1:1197061/52‑1 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                   >  1:2737625/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                   >  1:2779412/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                   >  1:2933413/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                   >  1:858842/1‑68 (MQ=255)
                 tGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTAc                                                   >  1:911511/1‑68 (MQ=255)
                      tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:2736850/69‑1 (MQ=255)
                      tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:968432/69‑1 (MQ=255)
                      tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:1640326/69‑1 (MQ=255)
                       ttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                            <  1:472269/69‑1 (MQ=255)
                       ttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                            <  1:1061679/69‑1 (MQ=255)
                       ttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                            <  1:2154327/69‑1 (MQ=255)
                        tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:1891152/67‑1 (MQ=255)
                        tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:1203737/67‑1 (MQ=255)
                        tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                            >  1:1235240/1‑68 (MQ=255)
                        tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                            >  1:2784898/1‑68 (MQ=255)
                        tttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCa                                            >  1:1060416/1‑68 (MQ=255)
                              gATGCCGGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:3266705/1‑68 (MQ=255)
                              gATGCCCTCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:787229/1‑68 (MQ=255)
                              gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:1078686/1‑68 (MQ=255)
                              gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:3039718/1‑68 (MQ=255)
                              gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:2170766/1‑68 (MQ=255)
                              gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:33372/1‑68 (MQ=255)
                              gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:3350062/1‑68 (MQ=255)
                              gATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                      >  1:2437516/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:1963861/69‑1 (MQ=255)
                                            aCGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGa                       <  1:2285654/69‑1 (MQ=255)
                                            aCGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGa                       <  1:1813946/69‑1 (MQ=255)
                                               tCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAAt                    <  1:30388/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:3248732/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:2557028/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:2797619/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:296943/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:1737017/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:2182327/69‑1 (MQ=255)
                                                ccGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:2098149/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:1256279/39‑1 (MQ=255)
                                                    ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:2209602/39‑1 (MQ=255)
                                                    ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:3396276/39‑1 (MQ=255)
                                                    ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:645017/39‑1 (MQ=255)
                                                    ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:1777116/39‑1 (MQ=255)
                                                    ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:556286/39‑1 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGGTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa              >  1:3371641/1‑69 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGc                                             <  1:3305092/38‑1 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa              >  1:2720623/1‑69 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa              >  1:1434324/1‑69 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa              >  1:2578371/1‑69 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa              >  1:904231/1‑69 (MQ=255)
                                                     tAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGa              >  1:1048341/1‑69 (MQ=255)
                                                        aGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATc                   <  1:2102108/61‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:2837528/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:850307/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:891151/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:603438/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:1006432/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:434242/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:410742/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:3435040/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:3232017/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:27734/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:2705195/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:2531267/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:1677124/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:2132796/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:2120397/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:1943505/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:1936989/69‑1 (MQ=255)
                                                                 tatTCCTTTCCCTGGACTACGCGAACAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa  <  1:1827687/69‑1 (MQ=255)
                                                                 tatCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGggagga                       <  1:423809/48‑1 (MQ=255)
                                                                    |                                                                 
GACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGA  >  minE/1076112‑1076245

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: