Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 60,568 G→A 43.1% 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,5680GA43.1% 25.8 / 62.1 51R125R (CGC→CGThepARNA polymerase‑associated helicase protein
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (8/14);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCAA  >  minE/60501‑60638
                                                                   |                                                                      
gcAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCgg                                                                     <  1:1837831/71‑1 (MQ=255)
           tACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATa                                                               <  1:1817016/66‑1 (MQ=255)
           tACGGCATACGGAACTGTT‑‑CTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATa                                                               <  1:4051416/64‑1 (MQ=255)
            aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                          >  1:823181/1‑70 (MQ=255)
            aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                          >  1:3307323/1‑70 (MQ=255)
            aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                          >  1:2260752/1‑70 (MQ=255)
            aCGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTc                                                          >  1:1189538/1‑70 (MQ=255)
             cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCa                                                         <  1:485460/70‑1 (MQ=255)
             cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCa                                                         <  1:3290342/70‑1 (MQ=255)
             cGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCa                                                         <  1:1993761/70‑1 (MQ=255)
              ggCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAAt                                                       <  1:4064867/71‑1 (MQ=255)
                 aTACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACgg                                                                     <  1:412212/54‑1 (MQ=255)
                    cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGccc                                                 >  1:1346844/1‑71 (MQ=255)
                    cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGccc                                                 >  1:672926/1‑71 (MQ=255)
                    cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGccc                                                 >  1:1871186/1‑71 (MQ=255)
                    cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGccc                                                 >  1:2175627/1‑71 (MQ=255)
                    cGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGccc                                                 >  1:462705/1‑71 (MQ=255)
                      gAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCgg                                               >  1:190824/1‑71 (MQ=255)
                      gAACTGTTCGCTGGAATATTTACGCGCGCGATAACGAAGCGCAAAACGGTCCATa                                                               <  1:2598813/55‑1 (MQ=255)
                        aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCa                                             <  1:1702231/71‑1 (MQ=255)
                        aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCa                                             <  1:1148324/71‑1 (MQ=255)
                        aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCa                                             <  1:183316/71‑1 (MQ=255)
                        aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCa                                             <  1:2948651/71‑1 (MQ=255)
                        aCTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCa                                             <  1:795378/71‑1 (MQ=255)
                          tGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCaaa                                           >  1:1567812/1‑71 (MQ=255)
                          tGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCaaa                                           >  1:2206139/1‑71 (MQ=255)
                          tGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCaaa                                           >  1:1494109/1‑71 (MQ=255)
                           gTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAAt                                                       >  1:1997040/1‑58 (MQ=255)
                           gTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCaaa                                           >  1:4220589/1‑70 (MQ=255)
                           gTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAAAGGTCCATACGGTCAATCTGCCCGGCaaa                                           >  1:375968/1‑70 (MQ=255)
                           gTTCGCTGAAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAAt                                                       >  1:2123092/1‑58 (MQ=255)
                            ttCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCaaa                                           <  1:2314138/69‑1 (MQ=255)
                            ttCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCaaa                                           <  1:1165828/69‑1 (MQ=255)
                            ttCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACa                                         <  1:2113698/71‑1 (MQ=255)
                             tCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAg                                        >  1:1165335/1‑71 (MQ=255)
                                  ggAATATTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAAt                                                       >  1:3782147/1‑51 (MQ=255)
                                   gAATATTTACGCGCGCGATAACGCAGCGCAAAACgg                                                                     <  1:3544628/36‑1 (MQ=255)
                                    aaTATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTcc                                 <  1:2874476/71‑1 (MQ=255)
                                     atatTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTcc                                 <  1:3045068/70‑1 (MQ=255)
                                     atatTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTcat                                >  1:93995/1‑69 (MQ=255)
                                     atatTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCt                                >  1:2761254/1‑71 (MQ=255)
                                     atatTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCt                                >  1:2684015/1‑71 (MQ=255)
                                     atatTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCt                                >  1:1209822/1‑71 (MQ=255)
                                      tatTTACGCGCGCGATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTg                               <  1:3159682/71‑1 (MQ=255)
                                          tACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCa                                                         >  1:2702437/1‑41 (MQ=255)
                                              cgcgcgATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGc                       <  1:3786494/71‑1 (MQ=255)
                                                   gATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCgg                            <  1:2237566/61‑1 (MQ=255)
                                                   gATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTg                     <  1:3983988/68‑1 (MQ=255)
                                                   gATAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAac                  <  1:557407/71‑1 (MQ=255)
                                                   gATAACGCAGCGCAAAACGGTCCATACGGGCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGa                    <  1:355604/69‑1 (MQ=255)
                                                     tAACGCAGCGCAAAACGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAaca                 >  1:3743280/1‑70 (MQ=255)
                                                                   gCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCaa  <  1:2739569/71‑1 (MQ=255)
                                                                   |                                                                      
GCAGGCCGCTGTACGGCATACGGAACTGTTCGCTGGAATATTTACGCGCGCGATAACGCAGCGCAAAGCGGTCCATACGGTCAATCTGCCCGGCAAACAGACGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCAA  >  minE/60501‑60638

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: