Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,758,201 G→A A554V (GCC→GTC)  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,2010GA81.2% 19.3 / 3.2 16A554V (GCC→GTC) selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (5/8);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

AATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCACGGTTCGTCACCAAACAGTGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTCTTCGCTGAAGCCCGCTTTGTGATCTGGCAGATG  >  NC_000913/3757965‑3758427
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
aaTTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTACACGACGACGAGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATc                                                                                                                                                                                                                      >  1:308212/1‑251 (MQ=255)
                                                               gAAATAAGAATGCGGCGCGCAATAAATGATCATTTCCACTACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCTAGCTTTCGGCCTACGCCTAAGCGATCGCGAAAACCCGCCGCGCAGGTTGAACCACACTCCTGATCGCGATCGCGGATCATATGGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTTCGCCGCCTgg                                                                                                                                                       <  1:77023/251‑1 (MQ=255)
                                                                                cgcAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAggn                                                                                                                                          <  1:204471/247‑2 (MQ=255)
                                                                                                        cgacgCGTAAAGCCAATGCGGTCAAAATATGCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACGCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGtc                                                                                                                <  1:246555/249‑1 (MQ=255)
                                                                                                                cccgCCACTGCGGGCAACATATTCCAGAATCTGAATTGCCAGCTGTCGGCCTACGCCTAAGCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGNCAANNNNNTc                                                                                                                                                                                                                 >  2:192844/4‑144 (MQ=255)
                                                                                                                aaaGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATAGTGACAAACTCGACAATCCTATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTGGCTGCGCCGTCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGcc                                                                                                       <  1:159835/250‑1 (MQ=255)
                                                                                                                                                gAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTTCCAGGTCACGCACCCACCACGGTTCGTCACCAAAc                                                                      <  1:127281/251‑1 (MQ=255)
                                                                                                                                                                            gCGATCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAg                                                                                                     >  1:155346/1‑192 (MQ=255)
                                                                                                                                                                               aTCGCGAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGAGCGTTACGGGAATAACGATCTTTAACGATCGCGGTAATTATTCCGTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTg                                                                                                                                <  1:143657/162‑1 (MQ=255)
                                                                                                                                                                                 cgcgAAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGCCAATCCGATCGTTACGGTCATAACGATATTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCACGGTTCGTCACCAAACAGTGGCTCTGCTTTTTGCCAGATGGCCtgctgc                                     >  1:101488/1‑251 (MQ=255)
                                                                                                                                                                                     gAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGCTCGCGGATCATATTGGCAAAc                                                                                                                                                                                                                               >  1:184995/1‑61 (MQ=255)
                                                                                                                                                                                              cgcgCAGGTTGAACCACACTCCTTATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTGAACGATCGCGGTAATTATCCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTGTTGCCAGGTCACGCACCCACCACGGTTCGTCACCAAACAGTGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTCTTCGCTGAAGc                        <  1:181842/251‑1 (MQ=255)
                                                                                                                                                                                               gcgcAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCg                                                                                                                                                                                                                            >  1:4195/1‑54 (MQ=255)
                                                                                                                                                                                               gcgcAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCg                                                                                                                                                                                                                            <  2:4195/54‑1 (MQ=255)
                                                                                                                                                                                                           ccACACTCCTGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCACGGTTCGTCACCAAACAGTGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTCTTCGCTGAAGCCCGCTTTGTGAtn           <  1:65818/251‑2 (MQ=255)
                                                                                                                                                                                                                    tGATCGAGATCGCGGATCATATTGACAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCACGGTTCGTCACCAAACAGTGGATCTGCTTTTTGCCAGATGGCCTGCTGCTCTTCGCTGAAGCCCGCTTTGTGATCTGGCAGATg  >  1:18064/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
AATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGGCAAACTCGACAATCCGATCGTTACGGTAATAACGATCTTTAACGATCGCGGTAATTATTCCTTGCTGCGCCGCCTGGCGTAGAGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCACGGTTCGTCACCAAACAGTGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTCTTCGCTGAAGCCCGCTTTGTGATCTGGCAGATG  >  NC_000913/3757965‑3758427

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