Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,831 A→T 12.5% 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.5% 106.4 / 8.4 40K102M (AAG→ATG) ‡PP_4720RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (15/20);  new base T (3/2);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

CTCGGAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATA  >  NC_002947/5367682‑5367979
                                                                                                                                                     |                                                                                                                                                    
cTCGGAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACaa                                                                                                                                                      <  2:95684/150‑1 (MQ=255)
    gAAGACCGCGAAGAACGCGCCGAGACCATCCCCGAGCTGTGCAAGGCCGGCCGGGCCGAGCTCGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACTTCCACCGTTACAAGTaa                                                                                                                                                  <  2:127733/150‑1 (MQ=255)
        aCCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGt                                                                                                                                              <  2:974251/150‑1 (MQ=255)
          cgcgAAGAACACGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAg                                                                                                                                            <  1:535217/150‑1 (MQ=255)
                   cgcgCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCa                                                                                                                                   >  1:645965/1‑150 (MQ=255)
                     cgcCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGt                                                                                                                                 <  2:645965/150‑1 (MQ=255)
                                tCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACTTGTAACGGtg                                                                                                                                             <  2:893032/127‑2 (MQ=255)
                                tCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACTTGTAACGGtg                                                                                                                                             >  1:893032/1‑126 (MQ=255)
                                    cGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAAcgcg                                                                                                                  >  2:990390/1‑150 (MQ=255)
                                     gAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACgcgc                                                                                                                 <  1:285926/150‑1 (MQ=255)
                                           tgCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTCCCAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGAc                                                                                                           >  1:613647/1‑150 (MQ=255)
                                            gCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGAcc                                                                                                          <  1:964470/150‑1 (MQ=255)
                                                   cggccgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCg                                                                                                   <  1:990390/150‑1 (MQ=255)
                                                      ccgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACTTGTAACGGtggat                                                                                                                                          >  1:1117286/1‑104 (MQ=255)
                                                      ccgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACTTGTAACGGtggat                                                                                                                                          <  2:1117286/108‑5 (MQ=255)
                                                       cgggCCGAGCTGGTGCAGCCCCTCCGAACGAAATCCCTGCTCCACTGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTg                                                                                               <  2:851200/147‑1 (MQ=255)
                                                       cgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGAAAGAAAAAAATCACTAGAACGCCACCCTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTg                                                                                               <  1:860942/150‑1 (MQ=255)
                                                        gTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTgg                                                                                              <  2:619152/150‑1 (MQ=255)
                                                        gTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTgg                                                                                              <  2:550339/150‑1 (MQ=255)
                                                                   gTGCAGACCAGCAGCAAGAAAGCCCTGATCTACCGCAAGAACCTACAGCCGAACAAGCAGCTGGCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGGCACCAGCAc                                                                                   >  2:415245/1‑150 (MQ=255)
                                                                                  aaGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGcc                                                                    >  1:884815/1‑150 (MQ=255)
                                                                                   aGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGccc                                                                   >  1:746072/1‑150 (MQ=255)
                                                                                           cTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGAt                                                                                                                                     >  1:1021996/1‑76 (MQ=255)
                                                                                           cTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGAt                                                                                                                                     <  2:1021996/76‑1 (MQ=255)
                                                                                                cTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAg                                                                      >  1:711739/1‑134 (MQ=255)
                                                                                                cTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAg                                                                      <  2:711739/134‑1 (MQ=255)
                                                                                                    cGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGaa                                                  <  2:1106920/150‑1 (MQ=255)
                                                                                                         gAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCt                                             >  1:119706/1‑150 (MQ=255)
                                                                                                             ccACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCa                                         >  1:98184/1‑150 (MQ=255)
                                                                                                                   ccGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGccacca                                                                                       >  1:874349/1‑98 (MQ=255)
                                                                                                                   ccGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGccacca                                                                                       <  2:874349/98‑1 (MQ=255)
                                                                                                                              gcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGaaa                        <  1:741696/150‑1 (MQ=255)
                                                                                                                                tGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATa                      >  1:1015038/1‑150 (MQ=255)
                                                                                                                                 gTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATAc                     >  2:855619/1‑150 (MQ=255)
                                                                                                                                    cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCAc                                                                                   <  2:1089102/85‑1 (MQ=255)
                                                                                                                                    cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCAc                                                                                   >  1:1089102/1‑85 (MQ=255)
                                                                                                                                    cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGa                                                   <  1:214317/117‑1 (MQ=255)
                                                                                                                                    cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGa                                                   >  2:214317/1‑117 (MQ=255)
                                                                                                                                          ctaCCGTTACTTGTAACGGTGGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCCCCATTACGATGCCGGAAAAGCACACCACCAGTACGGAATACTGTTGCCATCGCTCGCCCACGGAAATACCGTAGCGCa            >  1:1097435/3‑150 (MQ=255)
                                                                                                                                            aCCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAg          <  2:746072/150‑1 (MQ=255)
                                                                                                                                            aCCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAg          <  1:855619/150‑1 (MQ=255)
                                                                                                                                                    aaGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATa  <  2:119706/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCGGAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATA  >  NC_002947/5367682‑5367979

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: