Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 531,251 C→T intergenic (‑25/+44) selU ← / ← allS tRNA 2‑selenouridine synthase/DNA‑binding transcriptional activator AllS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,2510CT96.0% 72.7 / NA 26intergenic (‑25/+44)selU/allStRNA 2‑selenouridine synthase/DNA‑binding transcriptional activator AllS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/14);  minor base G (0/1);  total (11/15)
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

ACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGAT  >  NC_000913/531015‑531474
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                               
accaGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTAtttt                                                                                                                                                                                                                   <  2:434550/251‑1 (MQ=255)
        aTGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTTTTTTACGCGCCa                                                                                                                                                                                                           <  2:254621/251‑1 (MQ=255)
            cccAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGGATTGCGGCGGGCTCAAAATCGCAAAGGGG‑GCGGACAACAATAATGGGCGGATCAGCAAACATCCGGGGACGGTAGGCCGGGTCCGGGTGTCTTTCTTGCATAGGTAAAAATGAAACTCAAATCAGATTGTGCGCCATTTTATGCGCCACcgc                                                                                                                                                                                                       >  2:128703/1‑251 (MQ=255)
               agcgccagctcTGCGTCTGAGCGTTGCTTTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTAGTCATTAACGGCAGATTGATCGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACCTCAATAATGGGCGTATCAGCAAACAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCTCCAACGCGCa                                                                                                                                                                                                    <  2:368095/251‑1 (MQ=255)
                                       tGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTATGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGc                                                                                                                                                                              <  1:252811/249‑1 (MQ=255)
                                              tataGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCCAACTCG‑ATTGGGGCGCGGACATCAAAAAAGGGCGTATCAGCAAACAGCAGGGCCCGGTAGTCCTGTTCCGGGGGTCTCCCTTGCATTGGTTAAAGAGAACCCCCAATCAGTTTGCTCCCTATTTTACCCGCCACAGCGCCAGGGAACTTGATTTTTTACTTCGTggggttc                                                                                                                                                                     >  2:79088/1‑247 (MQ=255)
                                                             cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑CTAAGGGCGCGAACAACAATAATGGGGGGATCAGCAAACAGCAGGGCACGATAGTCCTGTTCCCTGTGTCTCTCTTGTATAGTTAAAAGTGAACCTCAAATCCGTTTGCGCGCTCTTTTTCGCGCCAACCCGCAAGGAAAATTGATTTTTAACTGCGTGGGTTGCCGCAAATTTTTaaa                                                                                                                                                       >  2:229085/1‑250 (MQ=255)
                                                             cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGcc                                                                                                                                                                                                            <  1:144606/197‑1 (MQ=255)
                                                                    gcgTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCCAACTCG‑ATAGGGGCGCGAACCTCAATAATGGGCGGATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTTCATAGTTAAAAGTGAACCTCAAATCAATTTGCGCGCTATTTTTAGCCCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGGTGCCGAAAATTTCTAAAAAACCGCt                                                                                                                                               >  1:319560/1‑251 (MQ=255)
                                                                                   cATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGAGCGAACATCAATAATGGGCGTATCAGCAATCAGTAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCCACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTgc                                                                                                                                 <  1:177491/250‑1 (MQ=255)
                                                                                                         gcggGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAat                                                                                                           >  1:104003/1‑250 (MQ=255)
                                                                                                                             aaaCTCG‑ATAGGGGCGCGAACATCAATAATGGGGGGATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGAGTGTCTCTCTTGCATAGTTAAAAGTGAACATCAAATCAGTTTGCGCGCTATTTTACGCGACAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAACTTTCTAAAAATCCGCTGATTACCGGCCTGCGCTGGGAAAACAAGGTCACAATATCTTCGACCGCTTTGcgg                                                                                        >  2:372146/1‑247 (MQ=255)
                                                                                                                                                aaCATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCGAAATTTGCGCCATGc                                                                    <  1:111782/251‑1 (MQ=255)
                                                                                                                                                              gCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTTATTTTTAACTGCGTGGGTTGCCGAAAATTGCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCg                                                      <  1:39047/251‑1 (MQ=255)
                                                                                                                                                                 tATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGa                                                     >  1:399503/1‑249 (MQ=255)
                                                                                                                                                                     agcaatcagcaGGGCACGATAGTACTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGCGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAaggagg                                               >  2:212648/1‑251 (MQ=255)
                                                                                                                                                                        aatcagcaGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGgcg                                            >  2:246928/1‑251 (MQ=255)
                                                                                                                                                                        aatcagcaGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGgcg                                            >  1:130841/1‑251 (MQ=255)
                                                                                                                                                                                gggCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTg                                    >  1:363065/1‑251 (MQ=255)
                                                                                                                                                                                gggCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTg                                    <  1:200163/251‑1 (MQ=255)
                                                                                                                                                                                   cACGATAGTCCTGGTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTGTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGa                                 <  2:306123/251‑1 (MQ=255)
                                                                                                                                                                                          gTCCCGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAGCAAGGTCACAATATCTTCGACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCg                          <  2:336349/251‑1 (MQ=255)
                                                                                                                                                                                                  ccGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGAcc                   <  1:138879/250‑1 (MQ=255)
                                                                                                                                                                                                   cGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAg                 <  1:246928/251‑1 (MQ=255)
                                                                                                                                                                                                                  cATAGTTAAAAGTGAACCTCAAATCAGTTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTTCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGAt  <  1:138135/251‑1 (MQ=255)
                                                                                                                                                                                                                  cATAGTTAAAAGTGAACCTCAAATCAGGTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTGTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCGTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGAt  <  2:29828/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                               
ACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCG‑ATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGAT  >  NC_000913/531015‑531474

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