Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,553,793 C→T 8.1% R154C (CGC→TGC)  PP_3138 → VirK domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,553,7930CT 8.1% 139.4 / 2.4 49R154C (CGC→TGC) PP_3138VirK domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (27/18);  new base T (2/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

AACCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGC  >  NC_002947/3553646‑3553935
                                                                                                                                                   |                                                                                                                                               
aaCCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCg                                                                                                                                                >  4:14771/1‑149 (MQ=255)
 aCCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGc                                                                                                                                               <  3:94705/149‑1 (MQ=255)
    tGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑c                                                                                                                                           >  3:60211/1‑149 (MQ=255)
                aCCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGaaa                                                                                                                                 >  3:86396/1‑147 (MQ=255)
                aCCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGaaa                                                                                                                                 <  4:86396/147‑1 (MQ=255)
                aCCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAAcc                                                                                                                               >  6:67016/1‑149 (MQ=255)
                  cTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAAcctg                                                                                                                             >  8:149671/1‑149 (MQ=255)
                           cACGCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGCTCAAGCGCGACGGCGAGCTGGGGCTGAACCTGTTCCAGCCCGACCTGCGGGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCAt                                                                                                                    <  4:264532/149‑1 (MQ=255)
                            aCTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                   >  2:196505/1‑149 (MQ=255)
                                            ccTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCT‑GGCTGCGGCCCAGGGTGaa                                                                                                                                  >  7:27062/1‑118 (MQ=255)
                                            ccTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCT‑GGCTGCGGCCCAGGGTGaa                                                                                                                                  <  8:27062/118‑1 (MQ=255)
                                             cTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCAt                                                                                                  >  7:26951/1‑149 (MQ=255)
                                               ggTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGGGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATcc                                                                                                <  7:121109/149‑1 (MQ=255)
                                                gTGCTGGTCCGGCCAATCTGGTTCAAGGGCGAAGGCGAGCGGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATGCCTTGCCTTCACGCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGGAGGGCATcc                                                                                                <  8:123661/148‑1 (MQ=255)
                                                gTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAggg                                                                                                                                      <  7:74816/110‑1 (MQ=255)
                                                gTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAggg                                                                                                                                      >  8:74816/1‑110 (MQ=255)
                                                      gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAgg                                                                                         >  4:154472/1‑149 (MQ=255)
                                                      gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAg                                                                                          <  3:14771/148‑1 (MQ=255)
                                                       aCCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTTGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGt                                                                                        <  7:149671/149‑1 (MQ=255)
                                                           gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAACCCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATc                                                                                     >  4:83358/1‑148 (MQ=255)
                                                                tCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTTTCgctgtc                                                                               >  6:203313/1‑144 (MQ=255)
                                                                           cgcgAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCt                                                                    <  1:196505/149‑1 (MQ=255)
                                                                            gcgAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTg                                                                   <  3:276654/149‑1 (MQ=255)
                                                                                  gCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATc                                                             >  1:244500/1‑149 (MQ=255)
                                                                                      gCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGCTGAAACCAGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTAcc                                                         >  2:250847/1‑149 (MQ=255)
                                                                                                   cTGTTCCAGGCCGACCCGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacct                                                  >  1:236361/1‑143 (MQ=255)
                                                                                                   cTGTTCCAGGCCGACCCGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacct                                                  <  2:236361/143‑1 (MQ=255)
                                                                                                        ccAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACt                                       >  5:173533/1‑149 (MQ=255)
                                                                                                          aGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGAc                                                                                  >  4:281757/1‑104 (MQ=255)
                                                                                                          aGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGAc                                                                                  <  3:281757/104‑1 (MQ=255)
                                                                                                           ggCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAg                                                                                          <  5:176895/95‑1 (MQ=255)
                                                                                                           ggCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAg                                                                                          >  6:176895/1‑95 (MQ=255)
                                                                                                            gCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACaaa                                                                                           >  5:54253/1‑93 (MQ=255)
                                                                                                            gCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACaaa                                                                                           <  6:54253/93‑1 (MQ=255)
                                                                                                            gCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCg                                    >  4:242546/1‑148 (MQ=255)
                                                                                                               gACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAgg                                >  5:190020/1‑149 (MQ=255)
                                                                                                                           gCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACg                    >  1:277831/1‑149 (MQ=255)
                                                                                                                              tCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAg                 >  4:97629/1‑149 (MQ=255)
                                                                                                                                ccTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGcc               >  4:74711/1‑149 (MQ=255)
                                                                                                                                   aggtttCACCCT‑GGCTGCGGCCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGt            >  8:79620/5‑149 (MQ=255)
                                                                                                                                     gtttCACCCT‑GGCTGCGGCCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTg          <  7:79620/147‑1 (MQ=255)
                                                                                                                                       ttCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCCGCCTGTTGCt        >  2:257681/1‑149 (MQ=255)
                                                                                                                                          aCCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCg      >  2:140569/1‑148 (MQ=255)
                                                                                                                                           cccTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGaa    <  8:26951/149‑1 (MQ=255)
                                                                                                                                           cccTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGaa    <  2:244500/149‑1 (MQ=255)
                                                                                                                                            ccTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAg   >  7:295957/1‑149 (MQ=255)
                                                                                                                                             cTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGt            >  4:139749/1‑139 (MQ=255)
                                                                                                                                             cTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGt            <  3:139749/139‑1 (MQ=255)
                                                                                                                                             cTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGc  <  4:60211/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AACCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGC  >  NC_002947/3553646‑3553935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑