Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,776,092 C→T 9.2% G42G (GGC→GGT)  PP_1583 → structural protein P5

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,776,0920CT 9.2% 108.0 / 2.0 43G42G (GGC→GGT) PP_1583structural protein P5
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (2/2);  total (24/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 = 8.30e-01

AGCCAGACCTACCCAGCGCCCGCCATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTCATAAACCGCT  >  NC_002947/1775943‑1776240
                                                                                                                                                     |                                                                                                                                                    
aGCCAGACCTACCCAGCGCCCGCCATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGgtgggc                                                                                                                                                      <  2:506389/150‑1 (MQ=255)
  ccAGACCTACCCAGCGCCCGCCATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGgtgggcgt                                                                                                                                                    <  1:876321/150‑1 (MQ=255)
    aGACCTACCCAGCGCCCGCCATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGCTTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGgtgggcgtgg                                                                                                                                                  >  2:870556/1‑150 (MQ=255)
         tACCCAGCGCCCGCCATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAg                                                                                                                                             >  1:101583/1‑150 (MQ=255)
                    cgccATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGc                                                                                                                                  >  1:1080601/1‑150 (MQ=255)
                        aTGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGa                                                                                                                              <  1:1111465/150‑1 (MQ=255)
                               ggCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGAcc                                                                                                                       <  1:870556/150‑1 (MQ=255)
                                  tttttttCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAg                                                                                                                    >  1:978565/1‑150 (MQ=255)
                                          gCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAAc                                                                                                            >  1:313825/1‑150 (MQ=255)
                                           ccTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTTGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACg                                                                                                           <  1:174671/150‑1 (MQ=255)
                                            cTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACgg                                                                                                          >  1:918499/1‑150 (MQ=255)
                                             tGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGt                                                                                                         >  1:505069/1‑150 (MQ=255)
                                              ggAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTa                                                                                                        <  2:588715/150‑1 (MQ=255)
                                               gAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTAt                                                                                                       <  1:104203/150‑1 (MQ=255)
                                                aGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATc                                                                                                      >  2:1055150/1‑150 (MQ=255)
                                                 gTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATcc                                                                                                     >  1:819612/1‑150 (MQ=255)
                                                  tAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGATCTGCCGGAGCTGTAGGTGGGCGTGGGCACGCCGCGTTTCGCCGGTTCAGTCCAGGCCGAGTACGGTATCCg                                                                                                    >  1:299522/1‑150 (MQ=255)
                                                    aaCCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGTGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAg                                                                                                  >  2:344444/1‑150 (MQ=255)
                                                    aaCCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGTGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAg                                                                                                  >  2:345023/1‑150 (MQ=255)
                                                      cccATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGTGTTGGCAAGCCGCGTTCCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGcc                                                                                                <  1:344444/150‑1 (MQ=255)
                                                      cccATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGTGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGcc                                                                                                <  1:345023/150‑1 (MQ=255)
                                                         aTGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTg                                                                                             >  1:832620/1‑150 (MQ=255)
                                                                 aCGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCt                                                                                     >  1:710942/1‑150 (MQ=255)
                                                                         tCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGc                                                                                                 <  1:1077375/130‑1 (MQ=255)
                                                                         tCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGc                                                                                                 >  2:1077375/1‑130 (MQ=255)
                                                                          ccGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGATCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTACACaa                                                                            <  2:835720/150‑1 (MQ=255)
                                                                           cGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAAc                                                                           <  2:489604/150‑1 (MQ=255)
                                                                           cGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAAc                                                                           >  2:721253/1‑150 (MQ=255)
                                                                                caacaacCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCGGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTAACCGAGCCCTGGGCAAGCTGTTGCTCAACTACCg                                                                      <  1:614409/150‑1 (MQ=255)
                                                                                 aacaacCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGACAAGCCGTTGCTCAACTACCgg                                                                     >  2:179538/1‑150 (MQ=255)
                                                                                         ggggAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAgga                                                             >  1:81772/1‑150 (MQ=255)
                                                                                           ggAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAggatg                                                           <  1:848304/150‑1 (MQ=255)
                                                                                                 tCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGc                                                     >  1:962268/1‑150 (MQ=255)
                                                                                                          aTCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTg                                            <  2:1091167/150‑1 (MQ=255)
                                                                                                              aCGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCAACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCACAATTCGGGGCAAGGATGGGATGCCTGGCGTTGGTCg                                        <  1:641406/150‑1 (MQ=255)
                                                                                                                    cGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTgg                                  <  2:832620/150‑1 (MQ=255)
                                                                                                                         ggCAGGGTCAGCTCGGGCTGGAGGTGGGCGCGGGCCAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCg                             <  2:299522/150‑1 (MQ=255)
                                                                                                                          gCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGa                            >  2:741690/1‑150 (MQ=255)
                                                                                                                             gggTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGacac                         >  1:223722/1‑150 (MQ=255)
                                                                                                                               gTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACAAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTc                       >  1:642915/1‑150 (MQ=255)
                                                                                                                                tCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAgga                                                             >  2:896668/1‑111 (MQ=255)
                                                                                                                                tCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAgga                                                             <  1:896668/111‑1 (MQ=255)
                                                                                                                                   gCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGACATCGACACTCCGCt                   >  2:520242/1‑150 (MQ=255)
                                                                                                                                   gCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCCTCCCGCCCGAGTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCt                   <  2:819612/150‑1 (MQ=255)
                                                                                                                                       gggCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAg               >  2:414645/1‑150 (MQ=255)
                                                                                                                                          cTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTc            <  1:149716/150‑1 (MQ=255)
                                                                                                                                                gtgggcgtgggcAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTCATAAAc      >  1:507062/1‑150 (MQ=255)
                                                                                                                                                  gggcgtgggcAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTCATAAACCg    <  2:918499/150‑1 (MQ=255)
                                                                                                                                                    gcgtgggcAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAggatgg                                                          >  2:368861/1‑94 (MQ=255)
                                                                                                                                                    gcgtgggcAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAggatgg                                                          <  1:368861/94‑1 (MQ=255)
                                                                                                                                                    gcgtgggcAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTCATAAACCGCt  >  1:137293/1‑150 (MQ=255)
                                                                                                                                                    gcgtgggcAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTCATAAACCGCt  <  2:962268/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCCAGACCTACCCAGCGCCCGCCATGAGCGGGCTTTTTTTCGCCTGGAGTAAACCCATGACTGCACGCGGTATCCGCAACAACAACCCGGGGAACATCGATTTCAATCCACGCAACGCCTGGCAGGGTCAGCTCGGGCTGGAGGTGGGCGTGGGCAAGCCGCGTTTCGCCCGATTCGACCAAGCCGAGAACGGTATCCGAGCCCTGGGCAAGCTGTTGCTCAACTACCGGGGCAAGGATGGGATGCCTGGCGTTGGTCGGCCTGGCATCGACACTCCGCTGGAGTTCATAAACCGCT  >  NC_002947/1775943‑1776240

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