Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,553,796 A→G 8.5% S155G (AGC→GGC)  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,7960AG 8.5% 129.6 / 3.0 47S155G (AGC→GGC) PP_3138VirK domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (26/17);  new base G (2/2);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGC  >  NC_002947/3553650‑3553935
                                                                                                                                                  |                                                                                                                                            
tGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑c                                                                                                                                           >  3:60211/1‑149 (MQ=255)
            aCCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGaaa                                                                                                                                 <  4:86396/147‑1 (MQ=255)
            aCCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGaaa                                                                                                                                 >  3:86396/1‑147 (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    <  2:244500/149‑1 (MQ=255)
                                                                                                                                       cccTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGaa    <  8:26951/149‑1 (MQ=255)
                                                                                                                                        ccTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAg   >  7:295957/1‑149 (MQ=255)
                                                                                                                                         cTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGt            <  3:139749/139‑1 (MQ=255)
                                                                                                                                         cTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGt            >  4:139749/1‑139 (MQ=255)
                                                                                                                                         cTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGc  <  4:60211/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAG‑CCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGC  >  NC_002947/3553650‑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: ‑