Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,210 A→T 50.0% intergenic (+177/+51) glxK → / ← allE glycerate kinase II/S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,2100AT50.0% ‑3.8 / 14.4 12intergenic (+177/+51)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base T (3/3);  total (6/6)
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

ATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAG  >  NC_000913/543137‑543346
                                                                         |                                                                                                                                        
atatatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTg                                                                                                                                    <  3:329080/80‑1 (MQ=255)
atatatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTg                                                                                                                                    >  4:329080/1‑80 (MQ=255)
 tatatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAAtctc                                                                               <  3:684787/132‑1 (MQ=255)
 tatatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAAtctc                                                                               >  4:684787/1‑132 (MQ=255)
                            tAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTgg                                                                                              <  3:121034/90‑1 (MQ=255)
                            tAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTgg                                                                                              >  4:121034/1‑90 (MQ=255)
                                                                cACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAAtata                                              <  1:464370/102‑1 (MQ=255)
                                                                cACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAAtata                                              >  2:464370/1‑102 (MQ=255)
                                                                 aCGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAg  <  1:279285/145‑1 (MQ=255)
                                                                 aCGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAg  >  2:279285/1‑145 (MQ=255)
                                                                 aCGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAg  <  3:312008/145‑1 (MQ=255)
                                                                 aCGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAg  >  4:312008/1‑145 (MQ=255)
                                                                         |                                                                                                                                        
ATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAG  >  NC_000913/543137‑543346

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