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,5980GA95.2% 65.1 / ‑3.5 21L117L (CTC→CTTudkuridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/12);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTC  >  NC_000913/2142352‑2142834
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
gATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCaaga                                                                                                                                                                                                                                           >  1:152352/1‑250 (MQ=255)
gATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGaa                                                                                                                                                                                                                                          <  1:212986/251‑1 (MQ=255)
     ggTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAatgatt                                                                                                                                                                                                                                     <  2:172744/251‑2 (MQ=255)
     ggTTTTTCCCGCCGCGCGGCACGATAATGGCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAatgatg                                                                                                                                                                                                                                     <  2:152352/251‑1 (MQ=255)
                        cACGATAATGTCCGCATATTGTTTAGAAAGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTACGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCg                                                                                                                                                                                                                  >  1:180984/1‑251 (MQ=255)
                            aTAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCAc                                                                                                                                                                                                              <  1:273127/251‑1 (MQ=255)
                                                  aaGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGt                                                                                                                                                                                        >  2:154954/1‑251 (MQ=255)
                                                              aTTGCAGGAACAGCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTg                                                                                                                                                                            <  2:265995/251‑1 (MQ=255)
                                                                                          ttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCcg                                                                                                                                                 <  1:259531/250‑1 (MQ=255)
                                                                                            gATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGt                                                                                                                                              <  2:17497/251‑1 (MQ=255)
                                                                                              tatTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGttt                                                                                                                                            >  1:310404/1‑251 (MQ=255)
                                                                                                                     tGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAagc                                                                                                                      >  1:293877/1‑250 (MQ=255)
                                                                                                                             cACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACtg                                                                                                             <  2:236972/251‑1 (MQ=255)
                                                                                                                                                   gATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGt                                                                                       >  1:235863/1‑251 (MQ=255)
                                                                                                                                                   gATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCGTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGt                                                                                       <  1:311864/251‑1 (MQ=255)
                                                                                                                                                                             gCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTcc                                                             <  2:235863/251‑1 (MQ=255)
                                                                                                                                                                               ggggTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCAt                                                           >  1:70067/1‑251 (MQ=255)
                                                                                                                                                                                            aaaTGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCTTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCtt                                              <  2:184414/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGCAACAAAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGAtg       <  2:293877/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aacaaAATGCCTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCAACGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGCGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTtc  >  1:205335/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  cTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGc                                                                                                                            <  1:6667/119‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTTCAAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGc                                                                                                                            >  2:6667/1‑119 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
GATGCGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTATGTTCAACATAGCTGTAAACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTC  >  NC_000913/2142352‑2142834

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