Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,210 A→T 10.2% 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,2100AT10.2% 152.7 / 8.4 49intergenic (+177/+51)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base T (2/3);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

ATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGC  >  NC_000913/543068‑543353
                                                                                                                                              |                                                                                                                                               
aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGt                                                                                                                                            <  4:54909/148‑1 (MQ=255)
      gTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTGGcttt                                                                                                                                      <  3:215419/148‑4 (MQ=255)
                cATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa                                                                                                                            >  3:207986/1‑148 (MQ=255)
                cATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCCGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAg                                                                                                                           <  1:23831/149‑1 (MQ=255)
                  tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATTGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACCACAAAAGCCACCGTATAAACCTTGGCAAAATCCAACATAGc                                                                                                                          >  4:146635/1‑148 (MQ=37)
                  tGCCGCATCCGGCATGAGTACTATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGc                                                                                                                          >  3:199135/1‑148 (MQ=255)
                    ccGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTa                                                                                                                        >  3:312215/1‑148 (MQ=255)
                         tGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAaata                                                                                                                  <  2:47084/149‑1 (MQ=255)
                              cATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAg                                                                                                             >  2:344831/1‑149 (MQ=255)
                                tGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGc                                                                                                           <  1:240833/149‑1 (MQ=255)
                                  aGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGa                                                                                                         <  2:236234/149‑1 (MQ=255)
                                   gTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAg                                                                                                        >  1:198651/1‑149 (MQ=255)
                                       cATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGt                                                                                                    <  2:151453/149‑1 (MQ=255)
                                               tcctGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACtt                                                                                             >  3:57675/1‑148 (MQ=255)
                                                  tGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAAt                                                                                         <  2:198651/149‑1 (MQ=255)
                                                        tCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTAc                                                                                   <  1:344831/149‑1 (MQ=255)
                                                                       atatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATAc                                                                                                                                               <  2:136280/74‑1 (MQ=255)
                                                                       atatTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATAc                                                                                                                                               >  1:136280/1‑74 (MQ=255)
                                                                           tGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                 <  4:199135/148‑1 (MQ=255)
                                                                              agagATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAg                                                             >  2:18782/1‑149 (MQ=255)
                                                                              agagATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAg                                                             <  2:304921/149‑1 (MQ=255)
                                                                              agagATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAg                                                             <  1:309394/149‑1 (MQ=255)
                                                                                agaTCATTTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTa                                                           >  2:218908/1‑149 (MQ=255)
                                                                                  aTCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTa                                                                                    >  2:322602/1‑122 (MQ=255)
                                                                                  aTCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTa                                                                                    <  1:322602/122‑1 (MQ=255)
                                                                                        tAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCa                                                                                                                                 >  4:366325/1‑71 (MQ=255)
                                                                                        tAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCa                                                                                                                                 <  3:366325/71‑1 (MQ=255)
                                                                                         aGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGc                                                   <  3:231762/148‑1 (MQ=255)
                                                                                             cTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAAc                                              >  2:126718/1‑149 (MQ=255)
                                                                                                 tAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCt                                           <  4:207986/148‑1 (MQ=255)
                                                                                                    gcgaagcgCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTcac                                       >  2:67729/1‑149 (MQ=255)
                                                                                                       aagcgCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTcaccac                                    <  2:386740/149‑1 (MQ=255)
                                                                                                         gcgcATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGc                                  <  1:126718/149‑1 (MQ=255)
                                                                                                                aGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGGACGCTTCACCACGCCCTacac                           >  2:127328/1‑149 (MQ=255)
                                                                                                                aGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTaca                            <  4:57675/148‑1 (MQ=255)
                                                                                                                  gCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGa                                                <  3:318765/126‑1 (MQ=255)
                                                                                                                  gCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGa                                                >  4:318765/1‑126 (MQ=255)
                                                                                                                   cAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCa                         >  3:21307/1‑148 (MQ=255)
                                                                                                                   cAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCa                         >  4:18344/1‑148 (MQ=255)
                                                                                                                        ttACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                 <  2:132190/103‑1 (MQ=255)
                                                                                                                        ttACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCttt                                                                 >  1:132190/1‑103 (MQ=255)
                                                                                                                             aaaaaaGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTaaa                                                         >  2:308618/1‑106 (MQ=255)
                                                                                                                             aaaaaaGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTaaa                                                         <  1:308618/106‑1 (MQ=255)
                                                                                                                               ttttgCCACGGTTTATACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTa            >  2:317718/5‑149 (MQ=255)
                                                                                                                                      aCGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAAt      <  4:21307/148‑1 (MQ=255)
                                                                                                                                      aCGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAAt      <  3:382803/148‑1 (MQ=255)
                                                                                                                                         gTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGGACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAgc  <  1:127328/149‑1 (MQ=255)
                                                                                                                                         gTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATaa    <  2:210557/147‑1 (MQ=255)
                                                                                                                                         gTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATaa    >  1:210557/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                               
ATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGC  >  NC_000913/543068‑543353

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