Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 58,995 G→A 22.6% L650L (CTG→TTG)  hepA ← RNA polymerase‑associated helicase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE58,9950GA22.6% 89.1 / 18.1 66L650L (CTG→TTG) hepARNA polymerase‑associated helicase protein
Reads supporting (aligned to +/- strand):  ref base G (12/39);  new base A (3/12);  total (15/51)
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.42e-01

CAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCA  >  minE/58934‑59057
                                                             |                                                              
cAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                                         <  1:1123174/69‑1 (MQ=255)
cAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                                         <  1:3000397/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1073528/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1729635/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1818819/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1865162/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1493557/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:149080/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1468782/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1418154/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:2703613/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:1228756/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:511756/69‑1 (MQ=255)
  gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                                       <  1:649676/69‑1 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                                         >  1:1962002/1‑64 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                                         >  1:130940/1‑64 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                    >  1:2223996/1‑69 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                    >  1:2045849/1‑69 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                    >  1:3198619/1‑69 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATAc                                                    >  1:927985/1‑69 (MQ=255)
     tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGACAGATAGTTAATCAGATCGTTg                                                         >  1:3246476/1‑64 (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:1407218/1‑69 (MQ=255)
                   cAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAg                                      >  1:1183622/1‑69 (MQ=255)
                           tCGGTCTCATCCGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCg                              >  1:2336275/1‑69 (MQ=255)
                           tCGGTCTCATCCGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCg                              >  1:2601154/1‑69 (MQ=255)
                           tCGGTCTCAACCGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCg                              >  1:1594538/1‑69 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTTTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:318766/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:83276/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:3110144/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:310784/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:3149819/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:2857940/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:27762/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:3194565/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:229537/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:549355/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:109174/69‑1 (MQ=255)
                                      cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt                   <  1:1934483/69‑1 (MQ=255)
                                             gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCa            <  1:72567/69‑1 (MQ=255)
                                             gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCa            <  1:3400213/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:1052019/69‑1 (MQ=255)
                                                     gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc    <  1:2964840/69‑1 (MQ=255)
                                                     gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc    <  1:2874743/69‑1 (MQ=255)
                                                     gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc    <  1:2646646/69‑1 (MQ=255)
                                                     gTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATc    <  1:2362958/69‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1578900/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:3029549/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1883653/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:78367/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1916157/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1402997/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1377088/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:446500/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:3426324/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:26313/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:3376809/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1205613/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:1251512/64‑1 (MQ=255)
                                                            aGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCAGGCAGGTGTGCTCAAATGCATCCa  <  1:888128/64‑1 (MQ=255)
                                                             gATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCa  <  1:536534/63‑1 (MQ=255)
                                                             |                                                              
CAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCA  >  minE/58934‑59057

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: