Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,009 A→T 14.9% Q40L (CAG→CTG)  cca → fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,0090AT14.9% 136.0 / 17.2 47Q40L (CAG→CTG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base T (2/5);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

TGTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACC  >  NC_000913/3201865‑3202106
                                                                                                                                                |                                                                                                 
tgtgTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACctgct                                                                                                >  4:1017986/1‑147 (MQ=255)
      tttGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgc                                                                                          <  3:1017986/148‑2 (MQ=255)
      tttGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgc                                                                                          <  1:589017/148‑2 (MQ=255)
      tttGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgc                                                                                          <  3:928629/148‑2 (MQ=255)
            gAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcag                                                                                                  >  2:255253/1‑134 (MQ=255)
            gAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcag                                                                                                  <  1:255253/134‑1 (MQ=255)
                    gatgGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgcc                                                                                         <  2:50814/135‑3 (MQ=255)
                    gatgGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgcc                                                                                         >  1:50814/1‑133 (MQ=255)
                      tgGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAgg                                                                                          >  1:241875/1‑132 (MQ=255)
                      tgGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAgg                                                                                          <  2:241875/132‑1 (MQ=255)
                            gAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                <  1:866441/120‑1 (MQ=255)
                            gAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                >  2:866441/1‑120 (MQ=255)
                              aGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                >  1:4276/1‑118 (MQ=255)
                              aGATTTATCTGGTCGGTGGTGCTGTTCGGGACGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                <  2:4276/118‑1 (MQ=255)
                                          tCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTc                                                                      >  2:845806/1‑132 (MQ=255)
                                          tCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGACTTTCCTGTGTTTc                                                                      <  1:845806/132‑1 (MQ=255)
                                             gtggtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTc                                                                               >  3:234037/1‑120 (MQ=255)
                                             gtggtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTc                                                                               <  4:234037/120‑1 (MQ=255)
                                                          ggATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAgg                                                                                              <  1:649354/92‑1 (MQ=255)
                                                          ggATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAgg                                                                                              >  2:649354/1‑92 (MQ=255)
                                                             tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtgtg                                                                          >  2:3699/1‑109 (MQ=255)
                                                             tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtgtg                                                                          <  1:3699/109‑1 (MQ=255)
                                                                 ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATc                                                                >  2:502494/1‑115 (MQ=255)
                                                                 ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATc                                                                <  1:502494/115‑1 (MQ=255)
                                                                     tAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCcg                                                                                       >  3:581270/1‑88 (MQ=255)
                                                                     tAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCcg                                                                                       <  4:581270/88‑1 (MQ=255)
                                                                         gCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagc                                                                                                 <  3:675680/74‑1 (MQ=255)
                                                                         gCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagc                                                                                                 >  4:675680/1‑74 (MQ=255)
                                                                          cTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTc                                                                      >  1:180816/1‑100 (MQ=255)
                                                                          cTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTc                                                                      <  2:180816/100‑1 (MQ=255)
                                                                                    aaGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAc                              <  1:279255/130‑1 (MQ=255)
                                                                                    aaGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAc                              >  2:279255/1‑130 (MQ=255)
                                                                                        cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 <  4:175712/91‑1 (MQ=255)
                                                                                        cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 <  2:755368/91‑1 (MQ=255)
                                                                                        cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 >  1:755368/1‑91 (MQ=255)
                                                                                        cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 >  3:175712/1‑91 (MQ=255)
                                                                                          gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                <  1:6613/138‑1 (MQ=255)
                                                                                          gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                >  2:6613/1‑138 (MQ=255)
                                                                                             aTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCt                                                                             >  1:507686/1‑74 (MQ=255)
                                                                                             aTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCt                                                                             <  2:507686/74‑1 (MQ=255)
                                                                                                   tggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACAcc  <  3:229650/143‑1 (MQ=255)
                                                                                                   tggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACAcc  >  4:229650/1‑143 (MQ=255)
                                                                                                    ggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                    <  1:313154/76‑1 (MQ=255)
                                                                                                    ggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                    >  2:313154/1‑76 (MQ=255)
                                                                                                                         ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagag                                               >  2:712958/1‑76 (MQ=255)
                                                                                                                         ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagag                                               <  1:712958/76‑1 (MQ=255)
                                                                                                                         ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttc            >  2:430921/1‑111 (MQ=255)
                                                                                                                         ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttc            <  1:430921/111‑1 (MQ=255)
                                                                                                                                  cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg          >  4:219782/1‑104 (MQ=255)
                                                                                                                                  cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg          <  3:219782/104‑1 (MQ=255)
                                                                                                                                     cgcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGa                      >  3:271498/1‑89 (MQ=255)
                                                                                                                                     cgcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGa                      <  4:271498/89‑1 (MQ=255)
                                                                                                                                        ggccTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCt                                       >  1:644357/4‑69 (MQ=37)
                                                                                                                                        ggccTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCt                                       <  2:644357/66‑1 (MQ=37)
                                                                                                                                          gCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg          >  3:590666/1‑96 (MQ=255)
                                                                                                                                          gCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg          <  4:590666/96‑1 (MQ=255)
                                                                                                                                                |                                                                                                 
TGTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACC  >  NC_000913/3201865‑3202106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: