Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 121,875 C→T 12.1% A18V (GCT→GTT)  cueO → multicopper oxidase (laccase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE121,8750CT12.1% 210.5 / 16.1 74A18V (GCT→GTT) cueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base C (23/42);  new base T (5/4);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCAC  >  minE/121805‑121944
                                                                      |                                                                     
aaGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGc                                                                       >  1:590621/1‑71 (MQ=255)
 aGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGc                                                                       <  1:1614044/70‑1 (MQ=255)
 aGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCt                                                                      <  1:1863462/71‑1 (MQ=255)
  gAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCtt                                                                     >  1:1489593/1‑71 (MQ=255)
   aaTAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt                                                                    <  1:84301/71‑1 (MQ=255)
     tAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGc                                                                  >  1:3417050/1‑71 (MQ=255)
     tAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCCTCGGCTTTGc                                                                  <  1:3414074/71‑1 (MQ=255)
       aGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGcc                                                                 <  1:1747676/70‑1 (MQ=255)
         gAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGTTTTGCCGCt                                                              <  1:1994481/71‑1 (MQ=255)
         gAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCt                                                              <  1:1491594/71‑1 (MQ=255)
           aaTAACTATGCAACGTCGTGATTTCTTAAAATATTGCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                                                            <  1:2424829/71‑1 (MQ=255)
           aaTAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                                                            >  1:1727743/1‑71 (MQ=255)
           aaTAACTACGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                                                            >  1:148638/1‑71 (MQ=255)
            aTAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           >  1:156789/1‑71 (MQ=255)
            aTAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           <  1:1982421/71‑1 (MQ=255)
            aTAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           <  1:3906097/71‑1 (MQ=255)
             tAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           <  1:1645349/70‑1 (MQ=255)
                  aTGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           <  1:2183505/65‑1 (MQ=255)
                   tGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGcc                                                      <  1:3843463/69‑1 (MQ=255)
                       aCGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt                                                                    <  1:184344/51‑1 (MQ=255)
                       aCGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                                                            <  1:1404052/59‑1 (MQ=255)
                       aCGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                                                            <  1:1805622/59‑1 (MQ=255)
                       aCGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                                                            <  1:2114768/59‑1 (MQ=255)
                        cgtcgtGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGt                                               <  1:3824551/71‑1 (MQ=255)
                        cgtcgtGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGt                                               <  1:3555955/71‑1 (MQ=255)
                        cgtcgtGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGt                                               <  1:1889771/71‑1 (MQ=255)
                         gtcgtGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa                                              >  1:33125/1‑71 (MQ=255)
                            gtGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTAtt                                            <  1:3565871/70‑1 (MQ=255)
                              gATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCg                                                                <  1:3482592/48‑1 (MQ=255)
                              gATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGc                                                               <  1:2129737/49‑1 (MQ=255)
                              gATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGc                                         >  1:1617110/1‑71 (MQ=255)
                              gATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGc                                         >  1:463178/1‑71 (MQ=255)
                               aTTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCg                                        >  1:2802018/1‑71 (MQ=255)
                                 ttCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           >  1:224514/1‑50 (MQ=255)
                                 ttCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           >  1:2303756/1‑50 (MQ=255)
                                   cTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                                                           >  1:1870221/1‑48 (MQ=255)
                                     tAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGaa                                  >  1:3174185/1‑71 (MQ=255)
                                      aaaaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAg                                    <  1:2075761/68‑1 (MQ=255)
                                      aaaaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGaa                                  <  1:3211729/70‑1 (MQ=255)
                                       aaaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAAc                                 <  1:3428258/70‑1 (MQ=255)
                                       aaaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAAc                                 <  1:3686485/70‑1 (MQ=255)
                                        aaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGc                               <  1:952732/71‑1 (MQ=255)
                                        aaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGc                               <  1:965033/71‑1 (MQ=255)
                                          tatTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg                                                          <  1:1038330/42‑1 (MQ=255)
                                          tatTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCa                                     >  1:1412000/1‑63 (MQ=255)
                                          tatTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGccc                             <  1:3144369/71‑1 (MQ=255)
                                           atTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg                                                <  1:1051607/51‑1 (MQ=255)
                                            ttCCGTCGCGCTGGGTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAg                                    <  1:1932405/62‑1 (MQ=255)
                                            ttCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCa                            <  1:1365535/70‑1 (MQ=255)
                                            ttCCGTCACGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCaa                           >  1:4064649/1‑71 (MQ=255)
                                               cGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGt                        <  1:113711/71‑1 (MQ=255)
                                                 tCGCGCTGGGTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTa                      <  1:3709621/71‑1 (MQ=255)
                                                 tCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa                                              <  1:3420762/47‑1 (MQ=255)
                                                   gcgcTGGGTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTAcc                    >  1:1007277/1‑71 (MQ=255)
                                                   gcgcTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTAcc                    <  1:3165303/71‑1 (MQ=255)
                                                      cTGGGTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGaa                                  >  1:441068/1‑54 (MQ=255)
                                                      cTGGGTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGaa                                  >  1:1314121/1‑54 (MQ=255)
                                                      cTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCAGCAGAACGCCCAACGTTACCGAt                 <  1:1265310/71‑1 (MQ=255)
                                                       tGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa                                              <  1:126781/41‑1 (MQ=255)
                                                         ggTGTGGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGt                        >  1:2534492/1‑61 (MQ=255)
                                                          gtgtGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGaa                                  <  1:1994296/50‑1 (MQ=255)
                                                          gtgtGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCt             >  1:3639307/1‑71 (MQ=255)
                                                             tgGCTTCGGTTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGa           <  1:917214/70‑1 (MQ=255)
                                                             tgGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGAt          >  1:2927940/1‑71 (MQ=255)
                                                               gCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGCACg                                >  1:142938/1‑47 (MQ=255)
                                                               gCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGAttt        >  1:290740/1‑71 (MQ=255)
                                                               gCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGAttt        >  1:2824085/1‑71 (MQ=255)
                                                               gCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGAtt         <  1:3445035/70‑1 (MQ=255)
                                                               gCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGAtt         <  1:178068/70‑1 (MQ=255)
                                                                cTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTg       <  1:3286603/71‑1 (MQ=255)
                                                                 ttCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGc      >  1:3517488/1‑71 (MQ=255)
                                                                   cGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCt     <  1:262364/70‑1 (MQ=255)
                                                                    ggCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTc    <  1:3793248/70‑1 (MQ=255)
                                                                    ggCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCa   >  1:3352028/1‑71 (MQ=255)
                                                                    ggCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCa   >  1:3469811/1‑71 (MQ=255)
                                                                     gCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCAc  >  1:1697881/1‑71 (MQ=255)
                                                                      |                                                                     
AAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCAC  >  minE/121805‑121944

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: