Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,693 C→G 8.8% P268R (CCG→CGG)  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,6930CG 8.8% 199.0 / 9.4 57P268R (CCG→CGG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base G (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGA  >  NC_000913/3202557‑3202833
                                                                                                                                        |                                                                                                                                            
cATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                                                                   <  4:769405/148‑1 (MQ=255)
cATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                                                                   <  3:562517/148‑1 (MQ=255)
       aaaTCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcatgg                                                                                                                            >  3:302872/1‑148 (MQ=255)
                cGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                   <  4:302872/148‑1 (MQ=255)
                       tCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtca                                                                                                                               >  1:103797/1‑129 (MQ=255)
                       tCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtca                                                                                                                               <  2:103797/129‑1 (MQ=255)
                         aTACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                   <  2:795460/139‑1 (MQ=255)
                         aTACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                   >  1:795460/1‑139 (MQ=255)
                           aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAg                                                                                                        <  4:312577/148‑1 (MQ=255)
                           aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAg                                                                                                        >  3:312577/1‑148 (MQ=255)
                              ttAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggt                                                                                                                                 >  4:559913/1‑120 (MQ=255)
                              ttAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggt                                                                                                                                 <  3:559913/120‑1 (MQ=255)
                              ttAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTgg                                                                                                     >  2:962878/1‑148 (MQ=255)
                               tAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCTGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                                                                           >  1:1055371/1‑109 (MQ=255)
                               tAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGa                                                                                                    <  1:962878/148‑1 (MQ=255)
                               tAATGACGCTCTCAATGGCGGAGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCACCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                                                                           <  2:1055371/109‑1 (MQ=255)
                                    aCGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                              >  1:739173/1‑133 (MQ=255)
                                    aCGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                              <  2:739173/133‑1 (MQ=255)
                                           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                                                                >  3:347990/1‑124 (MQ=255)
                                           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaac                                                                                               >  1:794575/1‑141 (MQ=255)
                                           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaac                                                                                               <  2:794575/141‑1 (MQ=255)
                                           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGCCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                                                                <  4:347990/124‑1 (MQ=255)
                                                             gTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAGAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcgg                                                                                                                       >  4:703942/1‑99 (MQ=255)
                                                             gTCCGCAGGTCGATGTCCGTTTCCCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcgg                                                                                                                       <  3:703942/99‑1 (MQ=255)
                                                               ccGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                                                                >  1:126972/1‑104 (MQ=255)
                                                               ccGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                                                                <  2:126972/104‑1 (MQ=255)
                                                                  cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                                                                 >  2:139389/1‑100 (MQ=255)
                                                                  cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                                                                 <  1:139389/100‑1 (MQ=255)
                                                                      tCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                            >  1:625172/1‑117 (MQ=255)
                                                                      tCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                            <  2:625172/117‑1 (MQ=255)
                                                                           gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGcc                                                                                                                          <  3:895419/82‑1 (MQ=255)
                                                                           gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGcc                                                                                                                          >  4:895419/1‑82 (MQ=255)
                                                                             ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                                                                     <  2:973220/85‑1 (MQ=255)
                                                                             ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                                                                     >  1:973220/1‑85 (MQ=255)
                                                                                 ttCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCaaa                                                                      <  4:790565/128‑1 (MQ=255)
                                                                                 ttCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCaaa                                                                      >  3:790565/1‑128 (MQ=255)
                                                                                     cAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                                                                <  2:690150/98‑1 (MQ=255)
                                                                                     cAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                                                                >  1:690150/1‑98 (MQ=255)
                                                                                           tATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                   >  2:600041/1‑73 (MQ=255)
                                                                                           tATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                   <  1:600041/73‑1 (MQ=255)
                                                                                                 aCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                               <  2:344119/71‑1 (MQ=255)
                                                                                                 aCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                               >  1:344119/1‑71 (MQ=255)
                                                                                                 aCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                              <  4:985794/72‑1 (MQ=255)
                                                                                                 aCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                              >  3:985794/1‑72 (MQ=255)
                                                                                                                       gcggccAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                            <  4:1093302/65‑1 (MQ=37)
                                                                                                                       gcggccAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                            >  3:1093302/4‑68 (MQ=37)
                                                                                                                       gcggccAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTg                                      >  1:106599/4‑122 (MQ=255)
                                                                                                                        cggccAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCa                               <  3:56630/126‑1 (MQ=255)
                                                                                                                        cggccAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCa                               >  4:56630/3‑128 (MQ=255)
                                                                                                                           ccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCaa        >  3:818133/1‑148 (MQ=255)
                                                                                                                             aGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                                                             <  3:480469/93‑1 (MQ=255)
                                                                                                                             aGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg                                                             >  4:480469/1‑93 (MQ=255)
                                                                                                                                aCTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATcc                    <  2:362735/131‑1 (MQ=255)
                                                                                                                                aCTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATcc                    >  1:362735/1‑131 (MQ=255)
                                                                                                                                 ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGa  <  4:818133/148‑1 (MQ=255)
                                                                                                                                   ctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgc                                                            >  4:1017443/1‑88 (MQ=255)
                                                                                                                                   ctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgc                                                            <  3:1017443/88‑1 (MQ=255)
                                                                                                                                   ctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTa                                                      <  4:405633/94‑1 (MQ=255)
                                                                                                                                   ctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTa                                                      >  3:405633/1‑94 (MQ=255)
                                                                                                                                        |                                                                                                                                            
CATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGA  >  NC_000913/3202557‑3202833

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