Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 60,568 G→A 46.0% R125R (CGC→CGT hepA ← RNA polymerase‑associated helicase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE60,5680GA46.0% 11.2 / 65.9 50R125R (CGC→CGThepARNA polymerase‑associated helicase protein
Reads supporting (aligned to +/- strand):  ref base G (19/8);  new base A (11/12);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCT  >  minE/60505‑60633
                                                               |                                                                 
gCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCa                                                            <  1:2880633/71‑1 (MQ=255)
  cGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATa                                                          >  1:154489/1‑71 (MQ=255)
  cGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATa                                                          >  1:727923/1‑71 (MQ=255)
   ccTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATAc                                                         >  1:2857088/2‑71 (MQ=255)
       tACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTc                                                              <  1:2942707/62‑1 (MQ=255)
        aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGGGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                     >  1:2064756/1‑70 (MQ=255)
        aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGt                                                      >  1:2910628/1‑69 (MQ=255)
        aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                     >  1:3101906/1‑70 (MQ=255)
        aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                     >  1:2994315/1‑70 (MQ=255)
        aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                     >  1:1462688/1‑70 (MQ=255)
        aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                     >  1:2470097/1‑70 (MQ=255)
         cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCa                                                    <  1:3500529/70‑1 (MQ=255)
         cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCa                                                    <  1:1312672/70‑1 (MQ=255)
         cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCaa                                                   <  1:1119029/71‑1 (MQ=255)
         cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCaa                                                   <  1:2407305/71‑1 (MQ=255)
         cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCaa                                                   >  1:1324110/1‑71 (MQ=255)
              tACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGc                                              >  1:3535222/1‑71 (MQ=255)
              tACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGc                                              >  1:1657259/1‑71 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGccc                                            >  1:3794265/1‑71 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGccc                                            >  1:3809005/1‑71 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGccc                                            >  1:861999/1‑71 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGcac                                            >  1:801746/1‑69 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGccc                                            >  1:2088791/1‑71 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGccc                                            >  1:362504/1‑71 (MQ=255)
                cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGccc                                            >  1:2038499/1‑71 (MQ=255)
                   aaCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTc                                                              <  1:823220/50‑1 (MQ=255)
                    aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCa                                        <  1:132161/71‑1 (MQ=255)
                    aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCa                                        <  1:3443404/71‑1 (MQ=255)
                     cTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCaa                                                   <  1:3534315/59‑1 (MQ=255)
                     cTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCaa                                                   <  1:2266357/59‑1 (MQ=255)
                          cGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGa                                  >  1:3127165/1‑71 (MQ=255)
                          cGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGa                                  >  1:2092263/1‑71 (MQ=255)
                               gAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACgg                               <  1:1736853/69‑1 (MQ=255)
                                    tttACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCg                        <  1:1636181/71‑1 (MQ=255)
                                    tttACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCg                        <  1:677243/71‑1 (MQ=255)
                                     ttACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCgg                       >  1:640213/1‑71 (MQ=255)
                                     ttACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCgg                       >  1:1935350/1‑71 (MQ=255)
                                            cgcgATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGt                              <  1:3603185/57‑1 (MQ=255)
                                              cgATAACGCAGCGCAAAACGGTCCATACGGTCAAt                                                  >  1:335879/1‑35 (MQ=255)
                                              cgATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGaa              >  1:118023/1‑71 (MQ=255)
                                                   aCGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGaa              >  1:659271/1‑66 (MQ=255)
                                                   aCGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGaa              >  1:1268821/1‑66 (MQ=255)
                                                   aCGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGaa              >  1:917587/1‑66 (MQ=255)
                                                   aCGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAacagaa         >  1:151164/1‑68 (MQ=255)
                                                   aCGCAGCGCAAAGCGGACCATAAGCTCAATCTGCCCGGCAAACAGACGGt                              >  1:3674938/1‑50 (MQ=255)
                                                          gcAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCt  <  1:2801426/71‑1 (MQ=255)
                                                          gcAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCt  <  1:1278601/71‑1 (MQ=255)
                                                           cAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGc                  <  1:3378488/54‑1 (MQ=255)
                                                           cAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGc                  <  1:251402/54‑1 (MQ=255)
                                                           cAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGc                  <  1:1982749/54‑1 (MQ=255)
                                                               |                                                                 
GCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCT  >  minE/60505‑60633

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: