Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,399,566 T→G 8.9% R108R (CGA→CGC) ‡ PP_3878 ← minor capsid protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,399,5660TG 8.9% 158.8 / 8.3 56R108R (CGA→CGC) ‡PP_3878minor capsid protein C
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base G (2/3);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

AGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTC  >  NC_002947/4399423‑4399706
                                                                                                                                               |                                                                                                                                            
aGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAg                                                                                                                                         >  4:278362/1‑149 (MQ=255)
aGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAg                                                                                                                                         <  3:247828/149‑1 (MQ=255)
aGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAg                                                                                                                                         >  6:228873/1‑149 (MQ=255)
   ccGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGccc                                                                                                                                      <  5:8801/149‑1 (MQ=255)
   ccGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGccc                                                                                                                                      <  8:167088/149‑1 (MQ=255)
         agtGGCTTGATCTGGGTCGCGGCCCGAATGGCCGCCGCCAGTTCAAAAGCGCCGCCGGCGCTGCCGCCGGGGCGGTCAATGTCGAGCACAATGTGTTCGACCATGGGGGCAGGGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATa                                                                                                                                <  8:82979/147‑1 (MQ=255)
         aTTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCAt                                                                                                                                 <  4:237215/148‑1 (MQ=255)
         aTTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTCCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCAt                                                                                                                                 <  7:71010/148‑1 (MQ=255)
          ttGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCtgaataa                                                                                                                               <  1:116885/149‑7 (MQ=255)
             gTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTg                                                                                                                            >  2:171307/1‑149 (MQ=255)
             gTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTg                                                                                                                            >  2:265319/1‑149 (MQ=255)
              tttGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCt                                                                                                                                     <  4:80110/139‑1 (MQ=255)
              tttGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCt                                                                                                                                     >  3:80110/1‑139 (MQ=255)
                tGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCg                                                                                                                                            >  2:214164/1‑130 (MQ=255)
                tGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCg                                                                                                                                            <  1:214164/130‑1 (MQ=255)
                       ggTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATg                                                                                                                      <  8:181809/145‑1 (MQ=255)
                       ggTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATg                                                                                                                      <  5:207453/145‑1 (MQ=255)
                           gcgGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCtgaa                                                                                                                                  <  6:230709/129‑4 (MQ=255)
                           gcgGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCtgaa                                                                                                                                  >  5:230709/1‑126 (MQ=255)
                                 cGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTAGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGtt                                                                                                        >  8:27934/1‑149 (MQ=255)
                                   aaTGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTc                                                                                                       >  8:11769/1‑148 (MQ=255)
                                         cgccgcCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGc                                                                                                >  1:159967/1‑149 (MQ=255)
                                                  ttCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGAc                                                                                       <  1:92951/149‑1 (MQ=255)
                                                     aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTca                                                                                    <  7:11769/149‑1 (MQ=255)
                                                     aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTca                                                                                    >  7:59479/1‑149 (MQ=255)
                                                     aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTc                                                                                     <  1:161228/148‑1 (MQ=255)
                                                     aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTc                                                                                     <  1:171307/148‑1 (MQ=255)
                                                        aGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACaa                                                                                 <  1:265319/149‑1 (MQ=255)
                                                                 ggCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCAt                                                                        >  5:34714/1‑149 (MQ=255)
                                                                      tGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACg                                                                   >  3:320379/1‑149 (MQ=255)
                                                                        ccgccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACg                                                                   <  6:236676/147‑1 (MQ=255)
                                                                           ccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGaa                                                                              <  5:263884/133‑1 (MQ=255)
                                                                           ccTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGaa                                                                              >  6:263884/1‑133 (MQ=255)
                                                                                 cTGTCAATGTAGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTAGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTc                                                        <  7:27934/149‑1 (MQ=255)
                                                                                    tCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGa                                                               <  2:156849/139‑1 (MQ=255)
                                                                                                  aaTGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGa                                       <  4:320379/149‑1 (MQ=255)
                                                                                                            aCCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGgcggcgcg                             >  6:306519/1‑149 (MQ=255)
                                                                                                               aTGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGcggcgcggcg                          <  5:228873/149‑1 (MQ=255)
                                                                                                                tGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCggcgcggcgc                         >  3:27094/1‑149 (MQ=255)
                                                                                                                 ggggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCt                        >  6:197195/1‑149 (MQ=255)
                                                                                                                  gggTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCt                        <  6:34714/148‑1 (MQ=255)
                                                                                                                    gTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGtt                     <  4:27094/149‑1 (MQ=255)
                                                                                                                         tGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTcg                >  2:66563/1‑149 (MQ=255)
                                                                                                                           aTGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTcgcg              <  3:278362/149‑1 (MQ=255)
                                                                                                                              gCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCgcgg             <  1:66563/147‑1 (MQ=255)
                                                                                                                               ccTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGa                                       <  3:316915/120‑1 (MQ=255)
                                                                                                                               ccTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGa                                       >  4:316915/1‑120 (MQ=255)
                                                                                                                                cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGAt         >  6:134205/1‑149 (MQ=255)
                                                                                                                                cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGAt         <  5:306519/149‑1 (MQ=255)
                                                                                                                                cTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGAt         >  5:221214/1‑149 (MQ=255)
                                                                                                                                    ttCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTg     >  7:293154/1‑149 (MQ=255)
                                                                                                                                      cAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTcg                <  3:231963/136‑1 (MQ=255)
                                                                                                                                      cAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTcg                >  4:231963/1‑136 (MQ=255)
                                                                                                                                       aGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCt   >  4:61630/1‑148 (MQ=255)
                                                                                                                                       aGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTc  <  5:197195/149‑1 (MQ=255)
                                                                                                                                        gggCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCt   >  4:89387/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                            
AGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGATGGCGGCGCGGCGCTGTTCTTCGCGGCGATCCTGCTC  >  NC_002947/4399423‑4399706

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: ‑