Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,142,598 G→A L117L (CTC→CTT)  udk ← uridine/cytidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,142,5980GA94.1% 43.0 / ‑5.4 17L117L (CTC→CTT) udkuridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (6/10);  minor base C (1/0);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCG  >  NC_000913/2142378‑2142814
                                                                                                                                                                                                                            |                                                                                                                                                                                                                        
cGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTc                                                                                                                                                                                            >  1:110866/1‑251 (MQ=255)
                               aaTGAATTGCAGGAACATCTGGCTCACGGTTTGTTGATAGTGCGCCATCACTGAATCCAATGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACTAAAATGGAGAAGTTAAGTGCGTCACTCAAACGCGCATTCTTCAGCAACAAAATGCCTTCAATAATCATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCCTACGCGTATGTTCAACAt                                                                                                                                                             <  1:304045/251‑1 (MQ=255)
                                        naGNAACNTNNGNNGCANGGTTTTTTGATATTGCGCCATNNCTGAATCCANTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATNCAGCGGGGTNNCANCGAAANTNNNGNNNTTNAGTTCGTCCCGCAAACGCGNATNNGTCAGCNTCAAAATGCCTTCAAGAATGATGAcc                                                                                                                                                                                                              <  2:74117/192‑1 (MQ=255)
                                           gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCg                                                                                                                                                 <  1:79936/251‑1 (MQ=255)
                                                  gggCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCTCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAcgcg                                                                                                                                                                         <  1:60379/220‑1 (MQ=255)
                                                            ttttttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCATCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCa                                                                                                                                              <  1:17572/237‑1 (MQ=255)
                                                                         cgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGTGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAAGGCCTTCAAGAATGATGACCTTCTTCGGCGCAACCGTCACCGGTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAAGTGCCGAGCCGCGTGTCAACg                                                                                                                   <  1:76238/251‑1 (MQ=255)
                                                                                                      gCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGACAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCCACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCAcc                                                                                                                                                                                         >  1:235188/1‑152 (MQ=255)
                                                                                                                      cTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCACCGAAAATGGAGCAGTTAAGTTCGTCACGCAAACGCGAATCCGTAAGCAACAAAATGCATTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGAAGGTCAATTGCCGAGACGCGTTTCa                                                                                                                      >  1:228612/1‑203 (MQ=255)
                                                                                                                         gATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGTCCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAAGGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATggn                                                                   <  1:33041/251‑2 (MQ=255)
                                                                                                                                  cATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTgg                                                          >  1:15947/1‑251 (MQ=255)
                                                                                                                                          aGATATCCAGCGGGGTATCAACGAAAATTGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTGCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGATCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTTTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAAcg                                                  <  1:253235/251‑1 (MQ=255)
                                                                                                                                                      nggTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATc                                      >  1:66803/2‑251 (MQ=255)
                                                                                                                                                                                     gTCACGCAAACGCGCATCCGTCAGCACCAACATGCCTTCCAGAATGATGCCCTTCTTCGGCTCAACCGCCACCGTTTCTTTCATACGCGTATGTTCAACATCGCTGCAAAACGGCAGGTCAATTGCCGAGCAGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGCTGGCTTTGATCTTTGTAATAGCAGt       >  1:139952/1‑251 (MQ=255)
                                                                                                                                                                                          gCAAACGCGCATCCGTCAGCAACACAATGACTTAAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGACTTCg  >  1:291098/1‑251 (MQ=255)
                                                                                                                                                                                                 cgcATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGc                                                                                                                               <  1:155364/119‑1 (MQ=255)
                                                                                                                                                                                                      ccGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtgtctttc                                                                                                                                                                                <  1:96223/65‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                        
CGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCG  >  NC_000913/2142378‑2142814

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