Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 122,440 C→T 24.0% A206A (GCC→GCT cueO → multicopper oxidase (laccase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE122,4400CT24.0% 33.4 / 5.7 29A206A (GCC→GCTcueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base C (9/13);  new base T (5/2);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCA  >  minE/122375‑122490
                                                                 |                                                  
gTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 <  1:413993/69‑1 (MQ=255)
gTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 <  1:1441017/69‑1 (MQ=255)
        tAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 >  1:3349423/1‑61 (MQ=255)
        tAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 >  1:2426090/1‑61 (MQ=255)
        tAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 <  1:1869157/61‑1 (MQ=255)
        tAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 <  1:2236370/61‑1 (MQ=255)
        tAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTg                                                 <  1:213195/61‑1 (MQ=255)
          aGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTg                                       <  1:883869/69‑1 (MQ=255)
          aGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTg                                       <  1:353830/69‑1 (MQ=255)
          aGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTg                                       <  1:2779892/69‑1 (MQ=255)
              aTTTAACGCTGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGa                                   >  1:1948117/1‑69 (MQ=255)
              aTTTAACGCTGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGa                                   >  1:792481/1‑69 (MQ=255)
              aTTTAACGCTGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGa                                   >  1:124062/1‑69 (MQ=255)
              aTTTAACGCTGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGa                                   >  1:2934983/1‑69 (MQ=255)
              aTTTAACGCTGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGa                                   >  1:2313563/1‑69 (MQ=255)
                ttAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATa                                 >  1:1584523/1‑69 (MQ=255)
                ttAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATa                                 >  1:2969962/1‑69 (MQ=255)
                                  ttgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCACCGGt                >  1:2719759/1‑68 (MQ=255)
                                  ttgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGt                >  1:1576074/1‑68 (MQ=255)
                                  ttgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGt                >  1:2164394/1‑68 (MQ=255)
                                  ttgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGt                >  1:678607/1‑68 (MQ=255)
                                  ttgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGt                >  1:842583/1‑68 (MQ=255)
                                   tgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTgg                                           <  1:160962/40‑1 (MQ=255)
                                   tgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTgg                                           <  1:2213750/40‑1 (MQ=255)
                                   tgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTgg                                           <  1:3104307/40‑1 (MQ=255)
                                   tgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTgg                                           <  1:2125439/40‑1 (MQ=255)
                                   tgattATCAACTGGATGTGATGACCGCCGCCGTGGGCTgg                                           <  1:1019753/40‑1 (MQ=255)
                                               ggATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCa  <  1:1873106/69‑1 (MQ=255)
                                                gATGTGATGACCGCCGCTGTGGGTTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCa  <  1:1551183/68‑1 (MQ=255)
                                                                 |                                                  
GTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCA  >  minE/122375‑122490

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: