Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,003 G→C 14.0% G38A (GGC→GCC)  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,0030GC14.0% 168.7 / 20.8 57G38A (GGC→GCC) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base C (3/5);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TTAGGAATTGTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACC  >  NC_000913/3201857‑3202106
                                                                                                                                                  |                                                                                                       
ttaGGAATTGTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGc                                                                                                        <  3:219349/148‑1 (MQ=255)
        tgtgTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACctgct                                                                                                >  4:1017986/1‑147 (MQ=255)
          tgtATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTAcctg                                                                                                  <  1:344841/144‑3 (MQ=255)
          tgtATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTAcctg                                                                                                  >  2:344841/1‑142 (MQ=255)
              tttGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgc                                                                                          <  3:1017986/148‑2 (MQ=255)
              tttGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgc                                                                                          <  3:928629/148‑2 (MQ=255)
              tttGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAgc                                                                                          <  1:589017/148‑2 (MQ=255)
                   cGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCggg                                                                                                         <  4:62105/128‑1 (MQ=255)
                   cGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCggg                                                                                                         >  3:62105/1‑128 (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                                                                                          <  2:241875/132‑1 (MQ=255)
                              tgGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAgg                                                                                          >  1:241875/1‑132 (MQ=255)
                                    gAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                >  2:866441/1‑120 (MQ=255)
                                    gAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                <  1:866441/120‑1 (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)
                                                                      gCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTa                                                                                                      >  1:868381/1‑80 (MQ=255)
                                                                      gCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTa                                                                                                      <  2:868381/80‑1 (MQ=255)
                                                                         ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATc                                                                <  1:502494/115‑1 (MQ=255)
                                                                         ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATc                                                                >  2:502494/1‑115 (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                                                                 >  1:755368/1‑91 (MQ=255)
                                                                                                cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 <  4:175712/91‑1 (MQ=255)
                                                                                                cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 <  2:755368/91‑1 (MQ=255)
                                                                                                cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                 >  3:175712/1‑91 (MQ=255)
                                                                                                  gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                >  2:6613/1‑138 (MQ=255)
                                                                                                  gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                <  1:6613/138‑1 (MQ=255)
                                                                                                     aTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCt                                                                             <  2:507686/74‑1 (MQ=255)
                                                                                                     aTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCt                                                                             >  1:507686/1‑74 (MQ=255)
                                                                                                           tggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACAcc  <  3:229650/143‑1 (MQ=255)
                                                                                                           tggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACAcc  >  4:229650/1‑143 (MQ=255)
                                                                                                            ggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                    >  2:313154/1‑76 (MQ=255)
                                                                                                            ggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                    <  1:313154/76‑1 (MQ=255)
                                                                                                                                 ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagag                                               <  1:712958/76‑1 (MQ=255)
                                                                                                                                 ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagag                                               >  2:712958/1‑76 (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)
                                                                                                                                                  gCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg          >  3:590666/1‑96 (MQ=255)
                                                                                                                                                  gCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg          <  4:590666/96‑1 (MQ=255)
                                                                                                                                                  |                                                                                                       
TTAGGAATTGTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACC  >  NC_000913/3201857‑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: