Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 123,017 G→A 68.0% G399S (GGC→AGC)  cueO → multicopper oxidase (laccase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE123,0170GA68.0% 87.0 / 55.4 106G399S (GGC→AGC) cueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (22/12);  new base A (38/34);  total (60/46)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 113 total aligned reads displayed.

GGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAT  >  minE/122959‑123084
                                                          |                                                                   
ggatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:3391515/69‑1 (MQ=255)
ggatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:1618223/69‑1 (MQ=255)
ggatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:3058714/69‑1 (MQ=255)
ggatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGAGGGAAGtt                                                           <  1:3075918/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACTGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:2393866/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTc                                                          >  1:3273021/1‑69 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:2925880/68‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:2755259/68‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:2443760/68‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGtt                                                           <  1:2176188/68‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:484395/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:3126417/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:2786958/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:249521/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:716559/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:2029661/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:1957384/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:1667478/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:2220926/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:1009285/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:806234/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:1467363/69‑1 (MQ=255)
 gatCACAGCCAGATGATGGGCCATATGGGGCACGACAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:1545529/69‑1 (MQ=255)
  atCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:485229/68‑1 (MQ=255)
  atCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:2854441/68‑1 (MQ=255)
  atCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:1755146/68‑1 (MQ=255)
              tgatgGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacca                                               >  1:1815653/1‑67 (MQ=255)
              tgatgGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCAt                                              >  1:1856592/1‑68 (MQ=255)
              tgatgGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCAt                                              >  1:2093885/1‑68 (MQ=255)
              tgatgGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCAt                                              >  1:1667383/1‑68 (MQ=255)
                      cATATGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:107865/48‑1 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCg                                                         >  1:1427091/1‑48 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTccacca                                               >  1:511775/1‑58 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTc                                                    >  1:309261/1‑53 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTc                                                    >  1:2640384/1‑53 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:2573985/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:3249406/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:912949/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:128899/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:2858804/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:2909124/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:3404958/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:1179103/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:2961381/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:1626854/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:1071865/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:951302/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:1472143/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:1573233/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:736397/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:827312/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:2411006/1‑59 (MQ=255)
                       atatGGGGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTCGATTTCCACCAt                                              >  1:133957/1‑59 (MQ=255)
                          tGGTGCACGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:497334/44‑1 (MQ=38)
                          tGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATc                                  >  1:1989777/1‑68 (MQ=255)
                          tGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCa                                 >  1:784100/1‑69 (MQ=255)
                            gggCACGGCAATATGAATCATATGAACAACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAc                               >  1:2102821/1‑69 (MQ=255)
                             ggCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACg                              >  1:3007289/1‑69 (MQ=255)
                             ggCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACg                              >  1:31671/1‑69 (MQ=255)
                             ggCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACg                              >  1:1111007/1‑69 (MQ=255)
                                 cGGCAATATGAATCATATGAACCACAGCGGGAAGTTc                                                          <  1:3244592/37‑1 (MQ=255)
                                   gCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATc                                  <  1:2455475/59‑1 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACg                              >  1:529014/1‑58 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACg                              >  1:487547/1‑58 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACg                              >  1:1702399/1‑58 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   <  1:640712/69‑1 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   <  1:1682858/69‑1 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   <  1:713942/69‑1 (MQ=255)
                                        atGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   <  1:25613/69‑1 (MQ=255)
                                         tGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   <  1:971746/68‑1 (MQ=255)
                                         tGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   <  1:2946037/68‑1 (MQ=255)
                                         tGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   >  1:1335535/1‑68 (MQ=255)
                                         tGAATCATATGAACCACGGCGGGAAGTTAGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGttt                   >  1:1846811/1‑68 (MQ=255)
                                           aaTCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGa                 <  1:73789/68‑1 (MQ=255)
                                           aaTCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGa                 <  1:404526/68‑1 (MQ=255)
                                           aaTCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGa                 <  1:1040621/68‑1 (MQ=255)
                                           aaTCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGa                 <  1:1412440/68‑1 (MQ=255)
                                           aaTCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGa                 <  1:2905632/68‑1 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGatat              >  1:1122908/1‑68 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGatat              >  1:1310097/1‑68 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:245123/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:1155207/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:1569198/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:192956/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:2358959/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:3316454/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:2473314/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:2792244/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:2921691/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:3313904/1‑69 (MQ=255)
                                              cATATGAACCACAGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:404494/1‑69 (MQ=255)
                                               atatGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATg             >  1:3192701/1‑68 (MQ=255)
                                                        aCGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGa   >  1:2965176/1‑69 (MQ=255)
                                                         cggcggGAAGTTCGATTTCCACCATGCCAACAAAATc                                  <  1:672034/37‑1 (MQ=255)
                                                         cggcggGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAt  <  1:2743310/69‑1 (MQ=255)
                                                         cggcggGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAt  <  1:1939467/69‑1 (MQ=255)
                                                         cggcggGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAt  <  1:739006/69‑1 (MQ=255)
                                                         cggcggGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAt  <  1:1692356/69‑1 (MQ=255)
                                                         cggcggGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAt  <  1:801129/69‑1 (MQ=255)
                                                         aggcggGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAt  >  1:598378/2‑69 (MQ=255)
                                                          |                                                                   
GGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACGGTCAGGCGTTTGATATGAACAAGCCGAT  >  minE/122959‑123084

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: