Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,831 A→T 12.8% K102M (AAG→ATG) ‡ PP_4720 → RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,8310AT12.8% 127.1 / 11.0 47K102M (AAG→ATG) ‡PP_4720RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base T (4/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAA  >  NC_002947/5367689‑5367970
                                                                                                                                              |                                                                                                                                           
gACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACgg                                                                                                                                      <  1:63856/150‑1 (MQ=255)
  ccGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACTTGTAACGGtg                                                                                                                                    <  1:384623/150‑2 (MQ=255)
   cgcgAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAg                                                                                                                                   >  2:258942/1‑150 (MQ=255)
    gcgAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGc                                                                                                                                  >  1:168761/1‑150 (MQ=255)
             gcgcCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTcc                                                                                                                               <  2:664170/144‑1 (MQ=255)
             gcgcCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTcc                                                                                                                               >  1:664170/1‑144 (MQ=255)
                       cATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTg                                                                                                               >  2:201456/1‑150 (MQ=255)
                        aTCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGa                                                                                                              >  1:207323/1‑150 (MQ=255)
                         tCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGaa                                                                                                             >  1:88741/1‑150 (MQ=255)
                               aGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGcc                                                                                                       <  2:227030/150‑1 (MQ=255)
                               aGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGcc                                                                                                       <  1:258942/150‑1 (MQ=255)
                                 cTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCt                                                                                                     <  1:307013/150‑1 (MQ=255)
                                       aaGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCgg                                                                                               <  1:201456/150‑1 (MQ=255)
                                                 gTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTgg                                                                                     <  1:537695/150‑1 (MQ=255)
                                                           ggTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGTAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCa                                                                           <  1:168148/150‑1 (MQ=255)
                                                           ggTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCa                                                                           <  1:214061/150‑1 (MQ=255)
                                                            gTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCAc                                                                          <  2:508634/150‑1 (MQ=255)
                                                                  aCCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGcc                                                                    >  2:711471/1‑150 (MQ=255)
                                                                    cATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCt                                                                                                     >  1:809349/1‑115 (MQ=255)
                                                                    cATCGACAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCt                                                                                                     <  2:809349/115‑1 (MQ=255)
                                                                             gAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAAcgcg                                                                                                         >  2:690020/1‑102 (MQ=255)
                                                                             gAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAAcgcg                                                                                                         <  1:690020/102‑1 (MQ=255)
                                                                               aaGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGc                                                       >  2:293090/1‑150 (MQ=255)
                                                                                aGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCa                                                      <  2:168761/150‑1 (MQ=255)
                                                                                    cTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACCATGCCGGAAAAGCCCAGCACCAg                                                  >  2:457310/1‑150 (MQ=255)
                                                                                     tGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGccac                                                                                >  2:700989/1‑119 (MQ=255)
                                                                                     tGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGccac                                                                                <  1:700989/119‑1 (MQ=255)
                                                                                           aCCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACAGCCACAGTGACAAGTAACGGTACCTCCGATCAGTGGCGCGGAGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAg                                           <  2:159465/150‑1 (MQ=255)
                                                                                                aaGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACAt                                      <  1:293090/150‑1 (MQ=255)
                                                                                                      ccACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCa                                >  2:808124/1‑150 (MQ=255)
                                                                                                       cacaGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAg                               >  1:723163/1‑150 (MQ=255)
                                                                                                              gAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGc                        >  2:99168/1‑150 (MQ=255)
                                                                                                                  aaGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCAc                    <  2:207323/150‑1 (MQ=255)
                                                                                                                    gcagcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTg                                                                                      <  2:540934/82‑1 (MQ=255)
                                                                                                                    gcagcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTg                                                                                      >  1:540934/1‑82 (MQ=255)
                                                                                                                      agcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCt                                                                                                     >  2:387206/1‑65 (MQ=255)
                                                                                                                      agcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCt                                                                                                     <  1:387206/65‑1 (MQ=255)
                                                                                                                             cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGt         >  1:834483/1‑150 (MQ=255)
                                                                                                                             cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGt         >  1:697400/1‑150 (MQ=255)
                                                                                                                               aCATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAg       >  1:642002/1‑150 (MQ=255)
                                                                                                                                cATCAAGCGTTACAAGTAACGGTAGCTCCGTTCTGTGGCGCGGTGAACGCGCCCTGACCGGTACCTGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAgc      <  1:639856/150‑1 (MQ=255)
                                                                                                                                  gctaCCGTTACTTGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAgcgc    >  1:725356/4‑150 (MQ=255)
                                                                                                                                   ctaCCGTTACTTGTAACGGTGGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCa   >  1:553062/3‑150 (MQ=255)
                                                                                                                                    taCCGTTACTTGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCaa  >  1:728080/2‑150 (MQ=255)
                                                                                                                                    cACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCaa  <  2:723163/150‑1 (MQ=255)
                                                                                                                                      ccGTTACTTGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCa                                            >  2:145969/1‑106 (MQ=255)
                                                                                                                                      ccGTTACTTGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCa                                            <  1:145969/106‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAA  >  NC_002947/5367689‑5367970

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