Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 59,005 T→C 29.3% V646V (GTA→GTG hepA ← RNA polymerase‑associated helicase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE59,0050TC29.3% 45.4 / 20.9 51V646V (GTA→GTGhepARNA polymerase‑associated helicase protein
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base C (3/12);  total (12/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGT  >  minE/58939‑59066
                                                                  |                                                             
tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                             >  1:927985/1‑69 (MQ=255)
tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                             >  1:3198619/1‑69 (MQ=255)
tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                             >  1:2223996/1‑69 (MQ=255)
tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                             >  1:2045849/1‑69 (MQ=255)
    aTCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCt                                                          <  1:1592710/68‑1 (MQ=255)
         aTCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTa                                                    <  1:3124187/69‑1 (MQ=255)
           cGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTaaa                                                  <  1:2629632/69‑1 (MQ=255)
              cAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATa                                                >  1:2702309/1‑68 (MQ=255)
              cAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAg                                               >  1:1183622/1‑69 (MQ=255)
              cAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAg                                               >  1:1407218/1‑69 (MQ=255)
                      tCGGTCTCATCCGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCg                                       >  1:2601154/1‑69 (MQ=255)
                      tCGGTCTCATCCGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCg                                       >  1:2336275/1‑69 (MQ=255)
                      tCGGTCTCAACCGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCg                                       >  1:1594538/1‑69 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTTTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:318766/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:3110144/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:310784/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:2857940/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:3149819/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:27762/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:3194565/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:549355/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:83276/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:109174/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:229537/69‑1 (MQ=255)
                                 cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                            <  1:1934483/69‑1 (MQ=255)
                                        gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCa                     <  1:3400213/69‑1 (MQ=255)
                                        gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCa                     <  1:72567/69‑1 (MQ=255)
                                        gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCa                     <  1:2023424/69‑1 (MQ=255)
                                              tAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCa               >  1:1212815/1‑69 (MQ=255)
                                               aGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCa               >  1:2756985/1‑68 (MQ=255)
                                                gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc             <  1:2964840/69‑1 (MQ=255)
                                                gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc             <  1:2874743/69‑1 (MQ=255)
                                                gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc             <  1:2362958/69‑1 (MQ=255)
                                                gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc             <  1:2646646/69‑1 (MQ=255)
                                                gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc             <  1:1052019/69‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1402997/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1205613/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:78367/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:446500/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:3426324/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:3376809/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1251512/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1377088/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1578900/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1883653/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:1916157/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:3029549/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:26313/64‑1 (MQ=255)
                                                       aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCAGGCAGGTGTGCTCAAATGCATCCa           <  1:888128/64‑1 (MQ=255)
                                                        gATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa           <  1:536534/63‑1 (MQ=255)
                                                             ttGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGt  <  1:2504419/67‑1 (MQ=255)
                                                                  |                                                             
TTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGT  >  minE/58939‑59066

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: