Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 122,068 C→T 22.9% R82R (CGC→CGT cueO → multicopper oxidase (laccase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE122,0680CT22.9% 68.1 / 13.5 52R82R (CGC→CGTcueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base C (9/31);  new base T (7/5);  total (16/36)
Fisher's exact test for biased strand distribution p-value = 3.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGC  >  minE/122015‑122129
                                                     |                                                             
aCTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACg                                                <  1:1488426/69‑1 (MQ=255)
aCTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACg                                                <  1:3201360/69‑1 (MQ=255)
aCTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACg                                                <  1:539650/69‑1 (MQ=255)
aCTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACg                                                <  1:3201453/69‑1 (MQ=255)
   aCCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGtt                                             <  1:2374823/69‑1 (MQ=255)
   aCCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGtt                                             <  1:812094/69‑1 (MQ=255)
   aCCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGtt                                             <  1:1645679/69‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGTGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:1977371/68‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:2682506/68‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:2159115/68‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:343665/68‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:420524/68‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:1061497/68‑1 (MQ=255)
          gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATc                                       <  1:1665883/68‑1 (MQ=255)
             ataACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTAca                                   >  1:3199157/1‑69 (MQ=255)
                aCGGCAATCTGCTGGGGCCGGCTGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTAcaacc                                >  1:214490/1‑69 (MQ=255)
                   gCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACg                                                <  1:2677733/50‑1 (MQ=255)
                                gCCGGCGGTGAAATTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAagag                  <  1:3199927/67‑1 (MQ=255)
                                gCCGGCGGTGAAATTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAagag                  <  1:2316910/67‑1 (MQ=255)
                                gCCGGCGGTGAAATTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAagag                  <  1:2789741/67‑1 (MQ=255)
                                gCCGGCGGTGAAATTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAagag                  <  1:2740404/67‑1 (MQ=255)
                                gCCGGCGGTGAAATTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAagag                  <  1:2563048/67‑1 (MQ=255)
                                       gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTg         <  1:2136818/69‑1 (MQ=255)
                                       gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTg         <  1:156982/69‑1 (MQ=255)
                                       gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTg         <  1:1301292/69‑1 (MQ=255)
                                       gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTg         <  1:1668392/69‑1 (MQ=255)
                                       gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTg         <  1:2465557/69‑1 (MQ=255)
                                       gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTg         >  1:2700743/1‑69 (MQ=255)
                                        tGAAATTACAGCGCTGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGc        <  1:1984351/69‑1 (MQ=255)
                                        tGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGc        <  1:1095163/69‑1 (MQ=255)
                                        tGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGc        <  1:440837/69‑1 (MQ=255)
                                           aaTTACAGCGCGGCAAAGCGGTAACGGTTTATATCTACAACCAACTGACGGAAGAGACAACGTTg         <  1:1136639/65‑1 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTTGCACt     >  1:2572282/1‑68 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTTGCACt     >  1:2314496/1‑68 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTTGCACt     >  1:2189357/1‑68 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTTGCACt     >  1:2065308/1‑68 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTTGCACt     >  1:195453/1‑68 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTTGCACt     >  1:583402/1‑68 (MQ=255)
                                            aTTACAGCGTGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACGACGTCGCACt     >  1:335631/1‑68 (MQ=255)
                                            aTTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGt           >  1:3220886/1‑62 (MQ=255)
                                            aTTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACt     >  1:651047/1‑68 (MQ=255)
                                              tACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGc  >  1:2454996/1‑69 (MQ=255)
                                              tACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGc  >  1:308800/1‑69 (MQ=255)
                                              tACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGc  >  1:2094815/1‑69 (MQ=255)
                                              tACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGc  >  1:330379/1‑69 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:899172/42‑1 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:757727/42‑1 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:1645975/42‑1 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:1679834/42‑1 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:1769399/42‑1 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:2455629/42‑1 (MQ=255)
                                               aCAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACt                            <  1:2724792/42‑1 (MQ=255)
                                                     |                                                             
ACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGC  >  minE/122015‑122129

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: