Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,761,969 A→T 20.7% N359Y (AAC→TAC) ‡ PP_1567 → capsid protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,761,9690AT20.7% 60.4 / 11.9 29N359Y (AAC→TAC) ‡PP_1567capsid protein
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base T (2/4);  total (12/17)
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.82e-01

TGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTT  >  NC_002947/1761825‑1762112
                                                                                                                                                |                                                                                                                                               
tGGGCGGACATCGAGCTGACCAAGGACGCCAAAAACGCTTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGCGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATc                                                                                                                                            >  1:452872/1‑150 (MQ=255)
    cGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCATcgcgc                                                                                                                                        <  2:171828/150‑2 (MQ=255)
     ttaCATCGAGCTGACCAAGGCCGTCACCGGCTCCTACATCTGGGTGCATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGc                                                                                                                                       <  2:283975/148‑1 (MQ=255)
       aCATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCATcgcgccta                                                                                                                                     <  2:196892/150‑5 (MQ=255)
         aTCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAg                                                                                                                                   <  2:11508/150‑1 (MQ=255)
                    cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAATGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACcg                                                                                                                        >  1:273242/1‑150 (MQ=255)
                    cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACcg                                                                                                                        >  2:104586/1‑150 (MQ=255)
                        gACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAg                                                                                                                    <  1:104586/150‑1 (MQ=255)
                             cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGaa                                                                                                                                                >  1:341012/1‑117 (MQ=255)
                             cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGaa                                                                                                                                                <  2:341012/117‑1 (MQ=255)
                              aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCg                                                                                                              >  2:533336/1‑150 (MQ=255)
                               aCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCgg                                                                                                             <  1:94916/150‑1 (MQ=255)
                                       tACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                     <  2:182619/150‑1 (MQ=255)
                                       tACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                     <  2:104550/150‑1 (MQ=255)
                                            cTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGtctc                                                                                                <  2:330787/150‑1 (MQ=255)
                                                         cAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGt                                                                                   >  2:111452/1‑150 (MQ=255)
                                                           ggAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                 >  1:408320/1‑150 (MQ=255)
                                                               ggcggcGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAAc                                                                             <  1:533336/150‑1 (MQ=255)
                                                                                                        cAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCct                                    <  1:399242/150‑1 (MQ=255)
                                                                                                        cAACGCGGTGCCAGAAGGCGAGTTCCTGCCAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCct                                    >  1:173775/1‑150 (MQ=255)
                                                                                                                    aGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGt                        >  1:539847/1‑150 (MQ=255)
                                                                                                                                 cTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag           <  2:408320/150‑1 (MQ=255)
                                                                                                                                   ggTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtc         <  1:111452/150‑1 (MQ=255)
                                                                                                                                     gggccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCTGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgaggag                                         >  2:28538/5‑116 (MQ=255)
                                                                                                                                     gggccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCTGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgaggag                                         <  1:28538/112‑1 (MQ=255)
                                                                                                                                       ggCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttc     <  2:381287/150‑1 (MQ=255)
                                                                                                                                        ccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctg                                   <  1:599880/118‑1 (MQ=255)
                                                                                                                                        ccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctg                                   >  2:599880/2‑119 (MQ=255)
                                                                                                                                          gcgATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTcgttcgtt  >  2:566785/1‑150 (MQ=255)
                                                                                                                                             aTGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                            >  1:533226/1‑105 (MQ=255)
                                                                                                                                             aTGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                            <  2:533226/105‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTT  >  NC_002947/1761825‑1762112

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