Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,642,037 G→A 25.5% V558I (GTT→ATT)  ggt → gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,642,0370GA25.5% 99.4 / 28.2 82V558I (GTT→ATT) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base G (42/19);  new base A (6/15);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 1.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

TGCTGGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCGGGG  >  minE/2641972‑2642102
                                                                 |                                                                 
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:988605/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:1609943/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:1484871/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:1956107/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:2092451/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:2129970/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:2188204/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:2229296/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:2294595/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:433742/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:910483/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:1170978/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCCATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:3319161/69‑1 (MQ=255)
tgctgGAAGCAAAAGGTCAGAAAGTGGCGATGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:1621999/69‑1 (MQ=255)
tgctgGAAGCAAAAGCTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTg                                                                <  1:1778047/69‑1 (MQ=255)
     gAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTGGGCCg                                                           >  1:2559400/1‑69 (MQ=255)
     gAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTGGGCCg                                                           >  1:2650433/1‑69 (MQ=255)
     gAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTGGGCCg                                                           >  1:2964487/1‑69 (MQ=255)
     gAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTGGGCCg                                                           >  1:2410770/1‑69 (MQ=255)
     gAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTGGGCCg                                                           >  1:1905693/1‑69 (MQ=255)
     gAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGATTGGGCCg                                                           >  1:1208514/1‑69 (MQ=255)
            aaGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGGTTGGGCCGGACGGTg                                                    <  1:3360148/69‑1 (MQ=255)
            aaGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGGTTGGGCCGGACGGTg                                                    <  1:1414058/69‑1 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGTGTCGGATGGTGa                                                   >  1:1942167/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:1439919/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:376881/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:3235241/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:1338811/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:316081/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:2072458/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:281193/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:936181/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:1259824/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:1646940/1‑69 (MQ=255)
             aGGTCAGAAAGTGGCGCTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGa                                                   >  1:1808599/1‑69 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:3103206/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:2905366/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:2734872/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:1298445/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:2534023/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:2920714/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:2913809/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:240136/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:639045/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:680447/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:1367870/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:827306/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:845707/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:2070673/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:938635/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:970480/1‑55 (MQ=255)
                            gcTGAATGAGGCGATGGGCAGCACGCAAAGCATTATGGTTGGGTCGGATGGTGAg                                                  >  1:1923160/1‑55 (MQ=255)
                                            ggCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTg                    >  1:1191284/1‑69 (MQ=255)
                                            ggCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTg                    >  1:2863551/1‑69 (MQ=255)
                                            ggCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTg                    >  1:2023896/1‑69 (MQ=255)
                                            ggCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTg                    >  1:1283925/1‑69 (MQ=255)
                                            ggCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTg                    >  1:2519690/1‑69 (MQ=255)
                                                 tACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGgatga               >  1:616323/1‑69 (MQ=255)
                                                   cgcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGATGAttt            <  1:2763443/67‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:725240/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:2856636/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:1478029/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:1968420/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:589609/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:2390944/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:1024403/68‑1 (MQ=255)
                                                    gcaAAGCATTATGGTTGGGTCGGATGGTGAGTTGTAAGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:369707/68‑1 (MQ=255)
                                                     caAAGCATTATGGTTGGGTCGGATGGTGAGTTGTACGGCGCATCCGACTCGCGCTTGGTGGATGAttt            <  1:1248063/68‑1 (MQ=255)
                                                     caAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:3342268/68‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGGATCCGACCCGCGCTCGGTGGATGATTTAACggggggt  >  1:2174843/1‑68 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:1444547/59‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:1084877/59‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:2880348/59‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:2430339/59‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:3145200/59‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGAttt            <  1:3143956/59‑1 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCgggg  >  1:2301417/1‑69 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCgggg  >  1:663346/1‑69 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCgggg  >  1:2194005/1‑69 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCgggg  >  1:2366551/1‑69 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCgggg  >  1:3041266/1‑69 (MQ=255)
                                                              aTGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCgggg  >  1:1708927/1‑69 (MQ=255)
                                                                 |                                                                 
TGCTGGAAGCAAAAGGTCAGAAAGTGGCGCTGAAAGAGGCGATGGGCAGTACACAAAGCATTATGGTTGGGCCGGACGGTGAGTTGTACGGCGCATCCGACCCGCGCTCGGTGGATGATTTAACGGCGGGG  >  minE/2641972‑2642102

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: