Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 583,396 T→C 10.2% V421A (GTG→GCG) ‡ selB → selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947583,3960TC10.2% 123.5 / 5.3 49V421A (GTG→GCG) ‡selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (3/2);  total (25/24)
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.46e-01

GCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTT  >  NC_002947/583256‑583543
                                                                                                                                            |                                                                                                                                                    
gctgcGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCAccc                                                                                                                                              <  2:153529/149‑1 (MQ=255)
    cgcgATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCgg                                                                                                                                            <  3:259080/147‑1 (MQ=255)
           cAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCaa                                                                                                                                   <  5:132603/149‑1 (MQ=255)
           cAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCa                                                                                                                                    <  7:271915/148‑1 (MQ=255)
               ggATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAActgc                                                                                                                               >  1:173475/1‑149 (MQ=255)
                          aGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAAc                                                                                                                   >  2:3322/1‑149 (MQ=255)
                              aCTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGCCACCACCCgcac                                                                                                                                          <  2:256883/123‑4 (MQ=255)
                              aCTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGCCACCACCCgcac                                                                                                                                          >  1:256883/1‑120 (MQ=255)
                                    aGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGc                                                                                                                 >  5:113517/1‑141 (MQ=255)
                                    aGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGc                                                                                                                 <  6:113517/141‑1 (MQ=255)
                                      cACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAg                                                                                                       <  7:122138/149‑1 (MQ=255)
                                        cTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGgc                                                                                                     >  3:11767/1‑149 (MQ=255)
                                        cTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCAGGCAGTTGCCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGgc                                                                                                     >  6:226480/1‑149 (MQ=255)
                                              gCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTc                                                                                                 >  7:166253/1‑147 (MQ=255)
                                                cAGTGCACCGGGCGGCATTAGCCCGTAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAag                                                                                              >  1:332837/1‑148 (MQ=255)
                                                 aGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATgg                                                                                                          <  6:201362/135‑1 (MQ=255)
                                                 aGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATgg                                                                                                          >  5:201362/1‑135 (MQ=255)
                                                 aGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAag                                                                                              >  1:184065/1‑147 (MQ=255)
                                                  gTGGACCGGGCGGGATTATCCCGCATCGACTCGAGCGGGATTTCAACCGCCTGGGGGAAACCTTGCATTTACCGGCTGATGTGCGGGTGGTGGCCACCCGGGATGGCCAACTGCTGTTTCGCCAACGCTCAATGGCAGGCGCTCAagca                                                                                            <  1:230363/149‑1 (MQ=255)
                                                      cgttggcggcATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAgg                                                                                       <  3:346234/145‑1 (MQ=255)
                                                      aCCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATgg                                                                                                          >  2:219404/1‑130 (MQ=255)
                                                      aCCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATgg                                                                                                          <  1:219404/130‑1 (MQ=255)
                                                           gcggcATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTg                                                                                  <  2:173475/149‑1 (MQ=255)
                                                                aTTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACg                                                                                                                  >  8:352161/1‑112 (MQ=255)
                                                                aTTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACg                                                                                                                  <  7:352161/112‑1 (MQ=255)
                                                                     cccGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTgg                                                                        >  4:195554/1‑149 (MQ=255)
                                                                      ccGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgc                                                                       >  8:210652/1‑149 (MQ=255)
                                                                       cGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgcg                                                                      >  5:242959/1‑149 (MQ=255)
                                                                        gCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgcgc                                                                     <  1:3322/149‑1 (MQ=255)
                                                                          aGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgcgcgc                                                                   <  1:180428/149‑1 (MQ=255)
                                                                                  gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCAc                                                             <  3:195554/147‑1 (MQ=255)
                                                                                  gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                           <  2:184065/149‑1 (MQ=255)
                                                                                       gCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGc                                                         >  6:296211/1‑146 (MQ=255)
                                                                                             cAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGccc                                                 <  6:45482/148‑1 (MQ=255)
                                                                                              aaCCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGa                                               >  3:167867/1‑149 (MQ=255)
                                                                                                    cTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCt                                         <  5:296211/149‑1 (MQ=255)
                                                                                                         ggAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGcc                                    <  5:173878/149‑1 (MQ=255)
                                                                                                                    aGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCACTCCAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCGCGCGCAAGAGCCCGACCAGCTTGGGCCTGACCGTGAt                         >  8:136173/1‑149 (MQ=255)
                                                                                                                       tACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGc                      >  8:227808/1‑149 (MQ=255)
                                                                                                                          cGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTg                   <  8:166253/149‑1 (MQ=255)
                                                                                                                           ggATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTgc                  >  5:162670/1‑149 (MQ=255)
                                                                                                                                cctgccGGGTGGCCACCACCC‑GCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCt            >  8:283271/6‑149 (MQ=255)
                                                                                                                                 tgtgCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGgcgc             <  7:227808/148‑1 (MQ=255)
                                                                                                                                 tgtgCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCt            >  3:353814/1‑149 (MQ=255)
                                                                                                                                  tgccGGGTGGCCACCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCtt           >  7:298059/4‑149 (MQ=255)
                                                                                                                                   tgCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCACTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCttt          >  8:178725/1‑149 (MQ=255)
                                                                                                                                   gccGGGTGGCCACCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCttt          <  8:298059/147‑1 (MQ=255)
                                                                                                                                    gCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTg         >  5:110296/1‑149 (MQ=255)
                                                                                                                                     cGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGt        <  4:167867/149‑1 (MQ=255)
                                                                                                                                        gtggtgGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGc     <  7:210652/149‑1 (MQ=255)
                                                                                                                                          ggtggCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGcc                                    <  4:363505/116‑1 (MQ=255)
                                                                                                                                          ggtggCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGcc                                    >  3:363505/1‑116 (MQ=255)
                                                                                                                                           gtggCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCtt  <  5:226480/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCGGGTGGTGGCCACCCGGCAGGGCCAACTGCTG‑TTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTT  >  NC_002947/583256‑583543

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: ‑