Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,079,889 G→A 11.1% A61V (GCC→GTC)  PP_2689 ← endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,079,8890GA11.1% 91.8 / 3.0 36A61V (GCC→GTC) PP_2689endoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (17/15);  new base A (2/2);  total (19/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 = 7.12e-01

CCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGACGACCGCCCATTCCAATGGCGCCCTGGA  >  NC_002947/3079741‑3080037
                                                                                                                                                    |                                                                                                                                                     
ccGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCgg                                                                                                                                                       <  1:97476/149‑1 (MQ=255)
        ggCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAgg                                                                                                                                               <  5:17510/149‑1 (MQ=255)
        ggCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAgg                                                                                                                                               <  1:406001/149‑1 (MQ=255)
          ccagccaCCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGgt                                                                                                                                             >  8:103718/1‑149 (MQ=255)
                 ccacTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGgcg                                                                                                                                       >  8:82560/1‑148 (MQ=255)
                    cTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAgg                                                                                                                                   <  2:132854/149‑1 (MQ=255)
                        ggCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGt                                                                                                                               >  4:153930/1‑149 (MQ=255)
                              gcgGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAg                                                                                                                                    <  1:219324/138‑1 (MQ=255)
                              gcgGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAg                                                                                                                                    >  2:219324/1‑138 (MQ=255)
                                     gggTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGc                                                                                                                   <  7:103718/148‑1 (MQ=255)
                                           ggTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCgggcctg                                                                                                            >  8:145804/1‑149 (MQ=255)
                                            gTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGGCGGTCACGTAGATCATCACCTGGGTCACGGCATCCAGCGAGCTGCCAGAGGCTTCCAGGGTGTGGCTCAGGTTGTCCAGTGTCTGGCgggccttg                                                                                                           <  5:366926/149‑3 (MQ=255)
                                                 tGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTc                                                                                                      >  1:313277/1‑149 (MQ=255)
                                                   gAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGa                                                                                                    <  7:82560/149‑1 (MQ=255)
                                                        aCTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGACTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATgccg                                                                                               >  4:354073/1‑149 (MQ=255)
                                                          tCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATgccgcc                                                                                             >  7:265396/1‑149 (MQ=255)
                                                            ggCTTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCTATgccgc                                                                                              <  4:185800/146‑1 (MQ=255)
                                                            ggCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCgg                                                                                           >  8:171604/1‑149 (MQ=255)
                                                            ggCGTACAAGGCATTTACAGTTTCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGAGGCTTCAAGTGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGCGCCTCGATGacgccgc                                                                                           >  1:164294/1‑148 (MQ=255)
                                                                             cAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGc                                                                         <  7:171604/149‑1 (MQ=255)
                                                                                  gCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGACTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGccgccg                                                                                            >  5:71985/1‑126 (MQ=255)
                                                                                  gCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGACTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGccgccg                                                                                            <  6:71985/126‑1 (MQ=255)
                                                                                       aTAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTg                                                               >  4:174512/1‑149 (MQ=255)
                                                                                       aTAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTg                                                               <  1:193224/149‑1 (MQ=255)
                                                                                             gcggTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGACTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCt                                                         <  3:354073/149‑1 (MQ=255)
                                                                                                    gtcACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAgt                                                  >  2:15590/1‑149 (MQ=255)
                                                                                                     tcACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGTAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAgtg                                                 >  2:320235/1‑149 (MQ=255)
                                                                                                     tcACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCa                                                    <  5:253726/146‑1 (MQ=255)
                                                                                                         gTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTAc                                              >  6:31925/1‑148 (MQ=255)
                                                                                                          tAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACaa                                            <  5:31925/149‑1 (MQ=255)
                                                                                                                 gCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCg                                     >  2:117948/1‑149 (MQ=255)
                                                                                                                     cTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGc                                 <  3:153930/149‑1 (MQ=255)
                                                                                                                            aCGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGGGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGTAACCCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCcgacga                           <  1:320235/149‑1 (MQ=255)
                                                                                                                              ggCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGACGACCg                        >  2:65821/1‑149 (MQ=255)
                                                                                                                                 aTCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGAcgcccgcc                      >  4:65213/1‑148 (MQ=255)
                                                                                                                                          gCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGACGACCGCCCATTCCAATg            <  1:117948/149‑1 (MQ=255)
                                                                                                                                                aGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGACGACCGCCCATTCCAATGGCGccc      >  1:28097/1‑149 (MQ=255)
                                                                                                                                                    gCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGACGACCGCCCATTCCAATGGCGCCCTGGa  <  1:15590/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTGGAAATACTCGGCGTACACGGCATTTACAGTTGCCAGATAACGGCGGTCGGTCACGTAGATCAGCACCTGGGTCACGGCATCCAGCGAGCTGCCAGCGGCTTCCAGGGTGTGGCGCAGGTTGTCCAGTGTCTGGCGGGCCTGGGCCTCGATGCCGCCGGCAA‑CCACTTGGCCCTGCGTATCGATTGGGATCTGCGCAGTGTACAACGTACCGTTGCCGACGACCGCCCATTCCAATGGCGCCCTGGA  >  NC_002947/3079741‑3080037

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