Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,758,202 C→A 9.4% A554S (GCC→TCC)  selB ← selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,758,2020CA 9.4% 101.1 / 4.4 43A554S (GCC→TCC) selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base A (2/2);  total (20/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.48e-01

ACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCAC  >  NC_000913/3758068‑3758341
                                                                                                                                      |                                                                                                                                           
tcggcgcgTATAGCGATTGTGGTCAAAATACGCCAGAATCGGAATGGCCAGCTGTCGGCCCACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGa                                                                                                                                     <  1:772844/139‑1 (MQ=255)
  gacgCGTAAAGCGGGTGCGGTCAAAATAGTCCAGAATCTGAAGGGCCAGGTTTCGGCGTACGCCGAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACa                                                                                                                                   <  1:771806/143‑1 (MQ=255)
  gacgCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTCAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGCTCATGTTCGCGCACTCcaca                                                                                                                                   >  1:331387/1‑143 (MQ=255)
  gacgCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGACCCCCGCTCCTGATCTAGATCGCGGATCATATTGGCAAACTCGACa                                                                                                                                   >  1:460147/1‑143 (MQ=255)
          aaGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATc                                                                                                                           <  1:914376/143‑1 (MQ=255)
          aaGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATc                                                                                                                           <  1:199651/143‑1 (MQ=255)
          aaGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCGTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATc                                                                                                                           <  2:111751/143‑1 (MQ=255)
                           atatTCCAGTATCTGAGTTGACAGCTTTCGGCCTACGCCTAAGCGATCGCGGACATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGAt                                                                                                          <  2:300786/143‑1 (MQ=255)
                               tCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCttt                                                                                                      <  2:133786/143‑1 (MQ=255)
                                 cAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTaa                                                                                                    >  2:40468/1‑143 (MQ=255)
                                         gAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGAGCATATTGGCAAACTCGACAATCCGATGGTTACGGTAAGAACGATCTTTACAGATcgcg                                                                                            >  1:146281/1‑143 (MQ=255)
                                              accAGCTTGCGGCCTAAGCCCAAGCGATCGCGGAAAGCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAat                                                                                       <  1:729320/142‑1 (MQ=255)
                                               ccAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAatt                                                                                      >  2:171245/1‑143 (MQ=255)
                                               ccAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAatt                                                                                      >  2:627893/1‑143 (MQ=255)
                                               ccAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGACCATATTGGCAAACTCGACAATCCGATCGTTACGGTACTAACGCTCTGTAACGATCGCAGTAatc                                                                                      >  2:469358/1‑142 (MQ=255)
                                                      gcGGCCTACGCCTACGCCATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGTCTCATTTTGGCATCCTCTACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCtt                                                                               <  2:403862/142‑1 (MQ=255)
                                                        ggCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGAAAACTCGACAATCCGATCGTTACGGTAATAACGATCtt                                                                                                       >  2:131130/1‑117 (MQ=255)
                                                        ggCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGAAAACTCGACAATCCGATCGTTACGGTAATAACGATCtt                                                                                                       >  2:575291/1‑117 (MQ=255)
                                                        ggCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGAAAACTCGACAATCCGATCGTTACGGTAATAACGATCtt                                                                                                       <  1:131130/117‑1 (MQ=255)
                                                        ggCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGAAAACTCGACAATCCGATCGTTACGGTAATAACGATCtt                                                                                                       <  1:575291/117‑1 (MQ=255)
                                                            tACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTgcg                                                                         >  1:498874/1‑143 (MQ=255)
                                                               gCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTACCGATCGCGGTAATTATTCCTTGCTGcgccg                                                                      >  1:858119/1‑143 (MQ=255)
                                                                 cTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCgccgcc                                                                    <  1:203869/143‑1 (MQ=255)
                                                                 cTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTAAGGTACCACCGATCTTTACCGCTCGCGGTATTTATTCCTTGCTGCgccgcc                                                                    >  1:122164/1‑143 (MQ=255)
                                                                   caGCGATCGCGGAAATCCGCCGCGCAGGGTGCACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTg                                                                  <  1:379087/142‑1 (MQ=255)
                                                                   aaGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGGTCGGGGTAATTATACCTGACTGCGCCGCCTg                                                                  >  2:183431/1‑143 (MQ=255)
                                                                     gCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGc                                                                <  1:410309/143‑1 (MQ=255)
                                                                     gCGATCGCGGAAATCCGCCGCGCAGGGTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCATTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGc                                                                <  1:169236/143‑1 (MQ=255)
                                                                       gATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGt                                                              >  2:286044/1‑143 (MQ=255)
                                                                        aTCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCCTATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGCTCTTTCACGATCGCGGTACTTATTCCTTGCTGCGCCCCCTGGCGTa                                                             >  1:605519/1‑143 (MQ=255)
                                                                          cgcgGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTaga                                                           <  1:162129/143‑1 (MQ=255)
                                                                             ggAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTc                                                        >  1:759830/1‑143 (MQ=255)
                                                                             ggAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCCTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTc                                                        >  1:754389/1‑143 (MQ=255)
                                                                                    cgccgcGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCa                                                 <  1:416330/143‑1 (MQ=255)
                                                                                                        acTCGTGATCGAGAGTGCGGATCATATTGGCAATCTCGACAATCCGATCGTTACCGTAATAGCGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTccc                             <  2:473090/143‑1 (MQ=255)
                                                                                                          tCCTGATCGTGATCGCGGATCATATTTGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGt                           <  1:528717/143‑1 (MQ=255)
                                                                                                               aTCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCtt                      >  1:544252/1‑143 (MQ=255)
                                                                                                                        gcgGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGt             <  1:627893/143‑1 (MQ=255)
                                                                                                                        gcgGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGt             <  1:171245/143‑1 (MQ=255)
                                                                                                                        gcgGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTGCCCGTCTCTTTTGCCAGGt             >  1:766114/1‑143 (MQ=255)
                                                                                                                             tCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGc        <  2:460147/143‑1 (MQ=255)
                                                                                                                               atatTGGCAAACTCGACAATCCGATCGTTACGTTAATAACGATATTTAACGATCGCGATAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTCAGCGTCTCTTTTGCCAGGTCACGcac      >  2:475622/1‑143 (MQ=255)
                                                                                                                                   tGGCCAACCCAACAAGCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGcacccac  <  1:563359/143‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                           
ACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCAC  >  NC_000913/3758068‑3758341

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