Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 58,975 T→C 30.5% S656S (TCA→TCG hepA ← RNA polymerase‑associated helicase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE58,9750TC30.5% 41.4 / 21.2 49S656S (TCA→TCGhepARNA polymerase‑associated helicase protein
Reads supporting (aligned to +/- strand):  ref base T (14/20);  new base C (3/12);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GCTTCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGT  >  minE/58916‑59040
                                                           |                                                                 
gcTTCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGa                                                          <  1:1405553/69‑1 (MQ=255)
gcTTCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGa                                                          <  1:763282/69‑1 (MQ=255)
   tCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCa                                                            >  1:2536428/1‑64 (MQ=255)
   tCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCa                                                            >  1:418591/1‑64 (MQ=255)
   tCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCa                                                            >  1:498825/1‑64 (MQ=255)
   tCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCa                                                            >  1:1922102/1‑64 (MQ=255)
           cTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCa                                                <  1:130144/68‑1 (MQ=255)
                  cAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                        <  1:1123174/69‑1 (MQ=255)
                  cAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                        <  1:3000397/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:1493557/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:1073528/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:511756/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:649676/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:1865162/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:1818819/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:1729635/69‑1 (MQ=255)
                    gTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGta                                      <  1:2703613/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:1228756/69‑1 (MQ=255)
                       tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                        >  1:1962002/1‑64 (MQ=255)
                       tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTg                                        >  1:130940/1‑64 (MQ=255)
                       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)
                       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:3110144/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:3149819/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:3194565/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:310784/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:2857940/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:27762/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:229537/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:549355/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:1934483/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:83276/69‑1 (MQ=255)
                                                        cGGCGAGGCCAGATAGTTAATCAAATCGTTGTACACGCTATCGTAGATAGTGCGTCCGGTCGGGCAGgt  <  1:109174/69‑1 (MQ=255)
                                                           |                                                                 
GCTTCATGTTGCTCGCGGCAGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGT  >  minE/58916‑59040

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: