Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,604,205 C→G 8.5% P12R (CCG→CGG) ‡ PP_4076 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,604,2050CG 8.5% 127.7 / 2.1 44P12R (CCG→CGG) ‡PP_4076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base G (2/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAC  >  NC_002947/4604059‑4604353
                                                                                                                                                   |                                                                                                                                                    
cGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGt                                                                                                                                                    >  2:296734/1‑149 (MQ=255)
   aCGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCg                                                                                                                                                 <  1:296734/149‑1 (MQ=255)
     ggggCGGGTGATCCGGGCGGGGACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTAGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGc                                                                                                                                                <  4:174272/146‑1 (MQ=18)
          gggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                          <  4:301580/149‑3 (MQ=255)
               aTCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAg                                                                                                                                      <  7:28819/148‑1 (MQ=255)
                tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGc                                                                                                                                    >  5:286395/1‑149 (MQ=255)
                       tGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGc                                                                                                                             >  5:264779/1‑149 (MQ=255)
                        ggTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGc                                                                                                                             >  7:5161/1‑148 (MQ=255)
                            ccAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGa                                                                                                                        <  8:58016/149‑1 (MQ=255)
                             cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGACGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAgg                                                                                                                         >  1:59308/1‑148 (MQ=255)
                             cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGa                                                                                                                        >  2:48956/1‑148 (MQ=255)
                             cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGa                                                                                                                        <  2:59308/148‑1 (MQ=255)
                               ggTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGGAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCgg                                                                                                                                         <  1:72425/129‑1 (MQ=255)
                               ggTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCgg                                                                                                                                         >  2:72425/1‑129 (MQ=255)
                                     gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                               >  7:41708/1‑149 (MQ=255)
                                     gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGCTAGGCGGGCGCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCtt                                                                                                                <  1:99548/148‑1 (MQ=255)
                                                  ccGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                                  <  4:225202/149‑1 (MQ=255)
                                                           aaGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCaa                                                                                          <  6:286395/148‑1 (MQ=255)
                                                                cATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCa                                                                                    <  1:48956/149‑1 (MQ=255)
                                                                cATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCa                                                                                    <  4:272058/149‑1 (MQ=255)
                                                                     gcg‑cGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTc                                                                                >  4:131192/1‑148 (MQ=255)
                                                                              ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                       <  6:50143/149‑1 (MQ=255)
                                                                               tCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  3:7126/1‑149 (MQ=255)
                                                                                cGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGATAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                       <  1:72127/147‑1 (MQ=255)
                                                                                cGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      <  7:86947/148‑1 (MQ=255)
                                                                                              tAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCt                                                                                 >  3:38746/1‑123 (MQ=255)
                                                                                              tAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCt                                                                                 <  4:38746/123‑1 (MQ=255)
                                                                                                ggCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAgg                                                      <  1:160365/148‑1 (MQ=255)
                                                                                                ggCCGGCTCGACACGCAATGAGGGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAgg                                                      <  2:279387/148‑1 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCagag                                          >  1:202603/1‑148 (MQ=255)
                                                                                                             cGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                        >  3:293802/1‑149 (MQ=255)
                                                                                                                       ggTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCt                              >  3:122742/1‑149 (MQ=255)
                                                                                                                         tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                            <  4:122742/149‑1 (MQ=255)
                                                                                                                            cttcGATCGCTTGATGGAGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGCCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGcc                         >  8:275491/1‑149 (MQ=255)
                                                                                                                                  tCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCTCCAGGAAATCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCATCCCCGAGATGGGCCTgcgc                    >  3:124477/1‑148 (MQ=255)
                                                                                                                                   tcccggATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGa                  >  2:188248/6‑149 (MQ=255)
                                                                                                                                   tcccggATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCATCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTCGGCCTGCGCGa                  >  8:68057/6‑149 (MQ=255)
                                                                                                                                    gCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGAc                 <  4:293802/149‑1 (MQ=255)
                                                                                                                                    gCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGAc                 <  8:5161/149‑1 (MQ=255)
                                                                                                                                     cTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAc                                                                                   <  5:187772/82‑1 (MQ=255)
                                                                                                                                     cTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAc                                                                                   >  6:187772/1‑82 (MQ=255)
                                                                                                                                        gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc              >  1:157029/1‑148 (MQ=255)
                                                                                                                                            gCGGACCCGTCCGCCATCCGGGAGGCGTCTCACCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGc                          >  2:109541/1‑132 (MQ=255)
                                                                                                                                                aCCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTGCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAt     <  8:41708/149‑1 (MQ=255)
                                                                                                                                                  ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGa   >  5:232911/1‑149 (MQ=255)
                                                                                                                                                   cGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAc  >  6:92454/1‑149 (MQ=255)
                                                                                                                                                   cGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAc  <  5:92454/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCG‑CGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAC  >  NC_002947/4604059‑4604353

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