Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 801,843 G→A 32.0% N43N (AAC→AAT dhaL ← dihydroxyacetone kinase, C‑terminal domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE801,8430GA32.0% 72.4 / 38.8 81N43N (AAC→AATdhaLdihydroxyacetone kinase, C‑terminal domain
Reads supporting (aligned to +/- strand):  ref base G (16/39);  new base A (9/17);  total (25/56)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAACCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGTAAGA  >  minE/801783‑801897
                                                            |                                                      
gAAACCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCata                                                  >  1:1044270/1‑67 (MQ=255)
gAAACCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCata                                                  >  1:692258/1‑67 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:3055151/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:496528/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:2246450/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:189799/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:33377/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               >  1:3282847/1‑69 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:1607565/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               >  1:3225136/1‑69 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:155148/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               >  1:1495967/1‑69 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:1349768/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               >  1:1335837/1‑69 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:810471/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               >  1:3189269/1‑69 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:2871361/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:1240578/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:3058413/69‑1 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               >  1:1118008/1‑69 (MQ=255)
 aaaCCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATAttt                                               <  1:296916/69‑1 (MQ=255)
         atCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTg                                       <  1:75613/69‑1 (MQ=255)
          tCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTg                                       <  1:3108239/68‑1 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:2756380/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:2928973/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:1036781/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:406567/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:426788/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:517286/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:605007/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:1274429/1‑67 (MQ=255)
           cTTTATCTGCGATAGCAGGGAGTTTTTCCACAACTTTGCTAAAGCCTCGATTCATATTTAGCCCGTg                                       >  1:656336/1‑67 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAg                                  <  1:2963315/67‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAg                                  <  1:3155693/67‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:2805308/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:1103830/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:114698/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:838937/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:1215110/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:127420/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:2722483/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:591906/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:3291983/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:1940841/68‑1 (MQ=255)
                tCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:2178116/68‑1 (MQ=255)
                         gCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAAtt                       >  1:3123347/1‑69 (MQ=255)
                         gCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAAt                        >  1:2100897/1‑68 (MQ=255)
                                 tttttCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATc               <  1:2870184/69‑1 (MQ=255)
                                 tttttCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATc               <  1:1858456/69‑1 (MQ=255)
                                 tttttCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATc               <  1:1325984/69‑1 (MQ=255)
                                   tttCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCa             >  1:305652/1‑69 (MQ=255)
                                   tttCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCa             >  1:3077090/1‑69 (MQ=255)
                                         ccacTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAAt                        >  1:2501877/1‑52 (MQ=255)
                                         ccacTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCg        >  1:1315207/1‑68 (MQ=255)
                                         ccacTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCg        >  1:2911026/1‑68 (MQ=255)
                                         ccacTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTAGCGATCCAGTcc         >  1:2862321/1‑67 (MQ=255)
                                          cacTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGt      <  1:1167752/69‑1 (MQ=255)
                                          cacTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGt      <  1:1677831/69‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGTaa    <  1:848233/69‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:896930/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:1293883/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:155932/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:59896/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:574073/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:906946/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:2077842/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:2145183/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:2493110/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:3244665/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:2717444/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:3207741/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:2828593/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:3176461/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:1019241/67‑1 (MQ=255)
                                            cTTTGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATccagtccagt      <  1:2904497/67‑1 (MQ=255)
                                             tttGCTAAAGCCTCGATTCATATTTAGCCCGTGGTCAGCGTCGCCAATCTCGCGATCCAGTCCTGt      <  1:880809/66‑1 (MQ=255)
                                              ttGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:2994650/38‑1 (MQ=255)
                                              ttGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGTAAGa  <  1:2043236/69‑1 (MQ=255)
                                              ttGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGTAAGa  <  1:2434993/69‑1 (MQ=255)
                                              ttGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGTAAGa  <  1:3151014/69‑1 (MQ=255)
                                                 cTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGc                                 <  1:156647/35‑1 (MQ=255)
                                                            |                                                      
GAAACCGATATCTTTATCTGCGATAGCAGGGAGTTTTTCCACCACTTTGCTAAAGCCTCGGTTCATATTTAGCCCGTGGTCAGCATCGCCAATTTCGCGATCCAGTCCGGTAAGA  >  minE/801783‑801897

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: