Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,207 A→T 12.1% intergenic (+174/+54) 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,2070AT12.1% 101.0 / 6.2 33intergenic (+174/+54)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (15/14);  new base T (2/2);  total (17/16)
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.75e-01

ATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCG  >  NC_000913/543110‑543354
                                                                                                 |                                                                                                                                                   
atCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                                                                                                                                            <  1:1268995/107‑1 (MQ=255)
atCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                                                                                                                                            >  2:1268995/1‑107 (MQ=255)
       ctGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                                                                                                                                                    >  1:338436/1‑92 (MQ=255)
       ctGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                                                                                                                                                    <  2:338436/92‑1 (MQ=255)
                            tatatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa                                                                                                                                 >  4:972190/1‑90 (MQ=255)
                            tatatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa                                                                                                                                 <  3:972190/90‑1 (MQ=255)
                               atTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTa                                                                          <  3:504324/142‑1 (MQ=255)
                               atTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTa                                                                          >  4:504324/1‑142 (MQ=255)
                                           aTGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGa                                                                                                 <  1:401390/107‑1 (MQ=255)
                                           aTGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGa                                                                                                 >  2:401390/1‑107 (MQ=255)
                                               aGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                  >  3:562803/1‑134 (MQ=255)
                                               aGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                  <  4:562803/134‑1 (MQ=255)
                                                    tGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGa                                                                                                 >  1:761132/1‑98 (MQ=255)
                                                    tGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGa                                                                                                 <  2:761132/98‑1 (MQ=255)
                                                    tGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTAc                                                                         >  2:934333/1‑122 (MQ=255)
                                                    tGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTAc                                                                         <  1:934333/122‑1 (MQ=255)
                                                          gcgaagcgCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGt                                                                                                     >  4:71748/1‑88 (MQ=255)
                                                          gcgaagcgCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGt                                                                                                     <  3:71748/88‑1 (MQ=255)
                                                                   aTCAGGCAATGTTACAAAAAAAGCCACGGTATaaa                                                                                                                                                 >  1:951370/1‑35 (MQ=255)
                                                                   aTCAGGCAATGTTACAAAAAAAGCCACGGTATaaa                                                                                                                                                 <  2:951370/35‑1 (MQ=255)
                                                                    tCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTg                                                                                                       <  1:726834/76‑1 (MQ=255)
                                                                    tCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTg                                                                                                       >  2:726834/1‑76 (MQ=255)
                                                                                   aaaaaaGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTg                                                                                                       >  4:1460522/1‑61 (MQ=255)
                                                                                   aaaaaaGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTg                                                                                                       <  3:1460522/61‑1 (MQ=255)
                                                                                          ccACGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCt                <  1:778812/141‑1 (MQ=255)
                                                                                          ccACGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCt                >  2:778812/1‑141 (MQ=255)
                                                                                           cACGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                  >  2:681635/1‑90 (MQ=255)
                                                                                           cACGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                  <  1:681635/90‑1 (MQ=255)
                                                                                           cACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAAt       >  2:461099/1‑149 (MQ=255)
                                                                                            aCGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTGCAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGcc                                  <  4:852338/121‑1 (MQ=255)
                                                                                            aCGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGcc                                  >  3:852338/1‑121 (MQ=255)
                                                                                            aCGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTaaa           >  3:843313/1‑144 (MQ=255)
                                                                                            aCGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTaaa           <  4:843313/144‑1 (MQ=255)
                                                                                               gTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTaaa                                                                                           >  2:1001280/1‑61 (MQ=255)
                                                                                               gTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTaaa                                                                                           <  1:1001280/61‑1 (MQ=255)
                                                                                                tataAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAgcg  >  3:158866/1‑149 (MQ=255)
                                                                                                 |                                                                                                                                                   
ATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCG  >  NC_000913/543110‑543354

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