Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,412,668 G→A 31.8% M21I (ATG→ATA aslA → acrylsulfatase‑like enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,412,6680GA31.8% 59.4 / 32.6 69M21I (ATG→ATAaslAacrylsulfatase‑like enzyme
Reads supporting (aligned to +/- strand):  ref base G (8/39);  new base A (12/10);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 3.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

TTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACG  >  minE/2412602‑2412726
                                                                  |                                                          
ttGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATg                                                            >  1:1362441/1‑67 (MQ=255)
          ttttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:1980313/69‑1 (MQ=255)
           tttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:1429789/68‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:2371832/67‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:959744/67‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:890427/67‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:584432/67‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:3087348/67‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:1898195/67‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:1822228/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:2312476/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              >  1:3262029/1‑69 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:355039/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:465792/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              >  1:1073180/1‑69 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:1909196/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:1767687/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:1755522/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:275135/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:605748/69‑1 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              >  1:2786528/1‑69 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              >  1:2790555/1‑69 (MQ=255)
            ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCGCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:2538550/69‑1 (MQ=255)
             tCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:1264356/68‑1 (MQ=255)
             tCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:2092804/68‑1 (MQ=255)
                ttttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                                                <  1:2599382/63‑1 (MQ=255)
                  ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAcccc                                        >  1:1818228/1‑69 (MQ=255)
                  ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAcccc                                        >  1:3100164/1‑69 (MQ=255)
                  ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAccc                                         >  1:1521256/1‑68 (MQ=255)
                  ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAcc                                          >  1:326179/1‑67 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:3205990/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:873795/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:2539798/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:2195057/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:660999/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:1550988/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:1678766/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:2257438/68‑1 (MQ=255)
                    cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:1734001/68‑1 (MQ=255)
                     aCCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg                                       <  1:1907313/67‑1 (MQ=255)
                       ccAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCa                                    <  1:570526/68‑1 (MQ=255)
                        cAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTACTCTCATGGCCAGCgcagca                                                <  1:2856553/55‑1 (MQ=255)
                            cGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGcc                              <  1:596100/69‑1 (MQ=255)
                                       tgCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:2595533/68‑1 (MQ=255)
                                       tgCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:2571554/68‑1 (MQ=255)
                                       tgCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:1288156/68‑1 (MQ=255)
                                       tgCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:194377/68‑1 (MQ=255)
                                       tgCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:2428542/68‑1 (MQ=255)
                                       tgCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:2836132/68‑1 (MQ=255)
                                        gCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                                              <  1:1391478/41‑1 (MQ=255)
                                                     ctctTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGAt     <  1:109176/69‑1 (MQ=255)
                                                        ttCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACcgcacgca                    <  1:1078471/51‑1 (MQ=255)
                                                          ccTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACg  <  1:1098227/67‑1 (MQ=255)
                                                          ccTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACg  <  1:895225/67‑1 (MQ=255)
                                                          ccTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACg  <  1:1373760/67‑1 (MQ=255)
                                                          ccTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACg  <  1:847897/67‑1 (MQ=255)
                                                          ccTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACg  <  1:2479877/67‑1 (MQ=255)
                                                          ccTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACg  <  1:2572331/67‑1 (MQ=255)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:2327553/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:2167918/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:1964561/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:1945796/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:2824330/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:1785677/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:1621358/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:1336678/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:3223864/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:277152/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgca                        >  1:1164192/1‑40 (MQ=38)
                                                               cATAGCCAGCGCAGCAGATACCCCGTCAACCGCCACcgc                         >  1:893928/1‑39 (MQ=38)
                                                                  |                                                          
TTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCCACCGCACGCAAAGGCTTTGCCGGATACG  >  minE/2412602‑2412726

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