Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 530,540 G→A C229C (TGC→TGT)  selU ← tRNA 2‑selenouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913530,5400GA91.7% 24.3 / ‑5.3 12C229C (TGC→TGT) selUtRNA 2‑selenouridine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (2/9);  minor base C (0/1);  total (2/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGGCATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTCCAGCACCCACAGGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCAGGTTTTCAAAACTCGCCTGGCTAAGTTGTGGTTGTAACGTGCGACCAAACGCCGAACCGCGATGACGCGCCAACCCTTCCAGATCAACACCGTTCGGCTGTTGCTGCACTAACAGCGTTTTACCGCTGCCGG  >  NC_000913/530329‑530775
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
cGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTGCGTTCAGGCGCTCAAGACCGATCTCAAACGGATCTTCGACCACCGCAATCGCCGCCTTAGTCATGCGCTCGCGCAGCCATTCCGGCAGGTGATTCGAACCGATCATCCGGCtttag                                                                                                                                                                                                      <  1:154410/251‑3 (MQ=255)
      aGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTcc                                                                                                                                                                                                <  1:155595/251‑1 (MQ=255)
       gcgcTGTAGCCCCAGCCGACGATTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCCAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTcc                                                                                                                                                                                                >  1:54273/1‑250 (MQ=255)
                                  cGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGTGTAAAATCAGGATGCATACGCAAGAAATACTCTGCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGGCAGTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTACGTCTTCCAGCACCCACAGGCGCAAATTCTGACgg                                                                                                                                                                     <  1:137658/250‑1 (MQ=255)
                                                                  aCTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATTATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCGTTCGTCTTCCAGCACCCACATGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCAgg                                                                                                                                    <  1:68095/251‑1 (MQ=255)
                                                                  aCTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTGTCGTCTTCCAGCACCCACAGGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCAgg                                                                                                                                    <  1:77966/251‑1 (MQ=255)
                                                                     ccTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACGCTTCGTGCATGCGCTCAAGACGGATCGCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCg                                                                                                                                                                                                      <  1:93774/182‑1 (MQ=255)
                                                                                 gTTCGTCTCTGTACGCGTGTGTAAAATCATTATGCATACGCTATAAATACTCTTCGTTCATGCGCTCAAGACGGATCTCAAACTGTCCTTCTACCACCGCACTCGCCTCCTGAGTCATTCGCTCACTCAGACTTTCCGGCAGGTGATTCGAACCTATCATCCGGCTTTCGTCTTCCATCACCCCCAGTCTCAAATTCTGCCGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCATGTTTTCAAATCTCg                                                                                                                       <  1:127455/249‑1 (MQ=255)
                                                                                                cgTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGAGCGTCTTCCACCGCAATCTCCGCCTGATTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTCCAGCAcc                                                                                                                                                                                          <  1:112197/167‑1 (MQ=255)
                                                                                                                                                        cGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTCCAGCACCCACAGGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCAGGTTTTCAAAACTCGCCTGGCTAAGTTGTGGTTGTACCGTGCGACCAAACGCCGAACCGCGATGACGCGCCAACCCTTCCAGATCAac                                              >  1:133275/1‑251 (MQ=255)
                                                                                                                                                                                           cgccTGCGTCATTCGCTCGCGCAGACATTCCTGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTCCAGCACCCACAGGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTGCGGCAGCCAGCAGGTTTTCAAAACTCGCCTGGCTAAGTTGTGGTTGTAACGTGCGACCAAACGCCGAACCGCGATGACGCGCCAACCCTTCCAGATCAACACCGTTCGGCTGTTGCTGCACTAACAGCGTTTTAc           <  1:118506/251‑1 (MQ=255)
                                                                                                                                                                                                    cATTCGCTCGCGCAGACATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTCCAGCACCCACAGGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCAGGTTTTCAAAACTCGCCTGGCTAAGTTGTGGTTGTAACGTGCGACCAAACGCCGAACCGCGATGACGCGCCAACCCTTCCAGATCAACACCGTTCGGCTGTTGCTGCACTAACAGCGTTTTACCGCTGCCgg  <  1:140142/251‑1 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
CGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTCATTCGCTCGCGCAGGCATTCCGGCAGGTGATTCGAACCGATCATCCGGCTTTCGTCTTCCAGCACCCACAGGCGCAAATTCTGACGGGCGTCGGTTTTTAGCATTTCGGCAGCCAGCAGGTTTTCAAAACTCGCCTGGCTAAGTTGTGGTTGTAACGTGCGACCAAACGCCGAACCGCGATGACGCGCCAACCCTTCCAGATCAACACCGTTCGGCTGTTGCTGCACTAACAGCGTTTTACCGCTGCCGG  >  NC_000913/530329‑530775

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