Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,399,270 Δ1 bp 8.2% coding (620/972 nt) 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,2700G. 8.2% 153.6 / 7.3 49coding (620/972 nt)PP_3878minor capsid protein C
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base . (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGT  >  NC_002947/4399124‑4399417
                                                                                                                                                  |                                                                                                                                                   
cAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGGCTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTtgtg                                                                                                                                                   <  6:222277/149‑1 (MQ=255)
cAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTtgtg                                                                                                                                                   <  8:149892/149‑1 (MQ=255)
      cACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCtt                                                                                                                                              >  7:303257/1‑148 (MQ=255)
           ccGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACAcc                                                                                                                                        >  2:14505/1‑149 (MQ=255)
             gcgGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACAcccg                                                                                                                                      >  3:382710/1‑148 (MQ=255)
                gTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGCCCGTTGTGACCTTCACACCCATGc                                                                                                                                   <  5:88507/149‑1 (MQ=255)
                 tATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGcc                                                                                                                                  >  7:261723/1‑149 (MQ=255)
                     cGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTAACACCCATGCCTTc                                                                                                                               >  3:32766/1‑148 (MQ=255)
                      ggCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCg                                                                                                                              <  4:32766/148‑1 (MQ=255)
                       gCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGat                                                                                                                            >  6:383764/1‑149 (MQ=255)
                       gCAACGTGGGTCGTGAACAGCTGATAGCTTTCTAGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGa                                                                                                                             >  1:54737/1‑148 (MQ=255)
                              gggTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGCTCATCTTc                                                                                                                     >  6:219729/1‑149 (MQ=255)
                               ggTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCGGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTc                                                                                                                     <  2:240558/148‑1 (MQ=255)
                                    tGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGAt                                                                                                              >  8:18135/1‑149 (MQ=255)
                                    tGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGAt                                                                                                              >  8:12215/1‑149 (MQ=255)
                                      aaCAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGAt                                                                                                              <  8:261723/148‑1 (MQ=255)
                                                 gCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGGCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTgg                                                                                                  <  6:100614/149‑1 (MQ=255)
                                                   tttCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGcc                                                                                                >  4:284804/1‑149 (MQ=255)
                                                   tttCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGc                                                                                                 >  4:151234/1‑148 (MQ=255)
                                                         gCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCAt                                                                                                          <  1:164293/133‑1 (MQ=255)
                                                         gCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCAt                                                                                                          >  2:164293/1‑133 (MQ=255)
                                                           aCGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCacac                                                                                         >  4:291024/1‑148 (MQ=255)
                                                                   gTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGAcc                                                                                <  3:315977/149‑1 (MQ=255)
                                                                      gAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCt                                                                              <  4:14932/148‑1 (MQ=255)
                                                                            ttGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGCCCTTCACACACATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCg                                                                        >  4:87852/1‑148 (MQ=255)
                                                                              gCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGa                                                                      <  3:151234/148‑1 (MQ=255)
                                                                              gCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGa                                                                      <  3:291024/148‑1 (MQ=255)
                                                                              gCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAg                                                                     >  8:204741/1‑149 (MQ=255)
                                                                               cAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAgg                                                                    <  3:284804/149‑1 (MQ=255)
                                                                               cAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAg                                                                     <  5:383764/148‑1 (MQ=255)
                                                                                 gagaTTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTc                                                                  <  7:204741/149‑1 (MQ=255)
                                                                                     ttGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGc                                                              <  8:303257/149‑1 (MQ=255)
                                                                                        cTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGa                                                                                  >  3:385682/1‑126 (MQ=255)
                                                                                               aTCGGCTCGTTCGGGGTCAGGTCATTTTTGTGTGCGCCGGCAATGACCGTTGTGACCTTCACAGCCATGCCTTCGATCATCTTCGTGCGATCCATTTGACTTGCCACCACACCTATGGACCCTACACCGGATGTCAGGCTAACCACGa                                                     <  2:75916/148‑1 (MQ=255)
                                                                                                   gCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCa                                                >  8:205627/1‑149 (MQ=255)
                                                                                                   gCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCa                                                <  6:369973/149‑1 (MQ=255)
                                                                                                    cTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAg                                               <  3:177432/149‑1 (MQ=255)
                                                                                                    cTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGt                                              <  7:12215/149‑1 (MQ=255)
                                                                                                    cTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGt                                              <  7:18135/149‑1 (MQ=255)
                                                                                                       gtgcGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCGTGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTg                                             <  4:382710/145‑1 (MQ=255)
                                                                                                       gTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACGCCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGc                                            <  2:54737/149‑1 (MQ=255)
                                                                                                           gggCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCTATTTATTTTATTGTTAAGCAGACCTATGTTATCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAg                                         <  8:62093/148‑1 (MQ=255)
                                                                                                               tCAGGTCATTTTTGTGTGCGCCGGAAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACAACGACTTCAGTGCAAGCAGAc                                    >  6:30771/1‑149 (MQ=255)
                                                                                                                    tCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAAt                               >  2:70666/1‑149 (MQ=255)
                                                                                                                       tttttGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAg                            <  1:14505/149‑1 (MQ=255)
                                                                                                                        ttttGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCCACCCGATGGCCCGTACCCCGCAGGTCGGGGTACCCACGCCTGCAGTGAAAGCAGCCGCAATCggg                           >  4:84341/1‑149 (MQ=255)
                                                                                                                              gtgCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAAccc                     >  5:2888/1‑149 (MQ=255)
                                                                                                                                gcgcCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAAccccc                   >  4:2734/1‑149 (MQ=255)
                                                                                                                                     ggCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGt              <  1:70666/149‑1 (MQ=255)
                                                                                                                                       cAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGTAc            >  8:122823/1‑149 (MQ=255)
                                                                                                                                         aaGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGc          <  6:2888/149‑1 (MQ=255)
                                                                                                                                                  gtgACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGt  <  7:205627/148‑1 (MQ=255)
                                                                                                                                                  gtgACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGt  >  1:222547/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAGTGCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGT  >  NC_002947/4399124‑4399417

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