Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,062 T→A 11.6% intergenic (+29/+199) glxK → / ← allE glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,0620TA11.6% 81.5 / 3.3 36intergenic (+29/+199)glxK/allEglycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (11/19);  new base A (2/2);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

CTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAA  >  NC_000913/542924‑543196
                                                                                                                                          |                                                                                                                                       
cTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCg                                                                                                                                    >  1:216239/1‑143 (MQ=255)
     cTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGTAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATaa                                                                                                                               <  2:26351/143‑1 (MQ=255)
         cACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATAAATCGCAtccgctt                                                                                                                           <  2:376267/143‑7 (MQ=255)
          aCCTGTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTGTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTg                                                                                                                          <  2:62393/143‑1 (MQ=255)
                                  cGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCAt                                                                                                  <  1:197192/143‑1 (MQ=255)
                                       aaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGGATGGGTa                                                                                             >  1:48271/1‑143 (MQ=255)
                                       aaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTa                                                                                             >  1:618156/1‑143 (MQ=255)
                                                 ctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcctt                                                                                   >  1:125127/1‑143 (MQ=255)
                                                    gaaCAGCGCGGGAAATATTGCCTGCGCCATTAAGATAGGTGAGGGAATTAAAAACTAACCCTTACCTTTGAAGCGGAGGCGATTTATATCGCAGAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcct                                                                                <  2:251945/142‑1 (MQ=255)
                                                           tcgtgAAATATCGCCTGCTCCATGAAAATAGGTCCGGGAATTATAAACTAACCCTTACCTTTACAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGCGTACCATATCCTTCCTGAAAATc                                                                         <  2:159905/139‑1 (MQ=255)
                                                           gcgcgAAATATTGCCTGCGCCATGAAAATAGGTCAGGCAACTAAAAACTAACCCTTACCTTTAAAGCGGATCCGATTTATATCGCATAAGAGTGCACTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                         <  2:285296/143‑1 (MQ=255)
                                                           gcgcgAAATATTGCCTGCGCAATTAAATTAGGTCAGTGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                         <  2:71842/143‑1 (MQ=255)
                                                              cgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgcg                                                                      >  1:170856/1‑143 (MQ=255)
                                                                    cttGCCTGCTCCATTAAAACAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAAtt                                                                <  1:602339/142‑1 (MQ=255)
                                                                     ttGCCTGCGACATTAATATAGGACAGGGAATTAAAAACTAACCCTTAGCTTTAAAGCGGATGCGATTTATAGCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTc                                                               <  1:102970/143‑1 (MQ=255)
                                                                           gcgcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTGAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGTTGCGGCATGAGTACCATATCCTTCCTGAAAAACAAGTAAATTCtatata                                                         >  1:201250/1‑143 (MQ=255)
                                                                              ctaTTTAAATATGTCAGTGAATTAAAAACGAACCCTTACCATAGAATTGGATTCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTg                                                      <  2:390768/141‑1 (MQ=255)
                                                                                aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTGTATATTGCa                                                    >  1:495093/1‑143 (MQ=255)
                                                                                  gataagcggTCAGGGGATTAAAAACTGACCCTTAGCTTGACAGCGGATGCGATTTATATCGCATGAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAGCCTTCCTGGATATCGCGCAAATTCTATATATTGCaga                                                  <  2:12844/136‑1 (MQ=255)
                                                                                       tAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCCGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                             <  1:247762/143‑1 (MQ=255)
                                                                                          gtcgtgcaaTTAAGAGCTACCCCTTACCTTTAAAGCGGAGTCTATTTATCTCGCATAAGAGTGCAGTACGCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                          <  2:353447/136‑1 (MQ=38)
                                                                                           tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGGTGCGAACAAAATACACAAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa                                         <  1:149235/143‑1 (MQ=255)
                                                                                           tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGCAAAAAGAGCATATTCTTTATCTTGCAGCGCTCATGTa                                         >  1:444597/1‑143 (MQ=255)
                                                                                               ggAATTAAATACTAACCCTTAGCGTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGcc                                     <  2:127441/143‑1 (MQ=255)
                                                                                               ggAATTAAAAGCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGcc                                     <  2:127444/143‑1 (MQ=255)
                                                                                                gagttAAAAACTAACCCTTACCGTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt                                    <  2:313092/140‑1 (MQ=255)
                                                                                                gacttAACACCTAACCATTACCTTTAAAGCGGATGCGATGTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTGCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt                                    <  2:190248/140‑1 (MQ=38)
                                                                                                      aaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaag                              >  1:345774/1‑143 (MQ=255)
                                                                                                               ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTCGGGCTGATgagcggagcgga                     >  2:166708/1‑141 (MQ=255)
                                                                                                                      tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGCGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGcc              >  2:261408/1‑142 (MQ=255)
                                                                                                                       ttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa             <  1:265668/143‑1 (MQ=255)
                                                                                                                          actcttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGt          >  1:20816/7‑143 (MQ=255)
                                                                                                                            gcgattGCGAGTTACAGCGCATAAGAGTGCAGTACTCCTGCCGGATGCGGCATGAGTACCATATCCTCCCTGAAACTCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa        <  2:221978/138‑1 (MQ=37)
                                                                                                                               taTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaa     >  1:390331/2‑143 (MQ=255)
                                                                                                                               gATGCGATTTACATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaa     <  1:39768/143‑1 (MQ=255)
                                                                                                                                  gCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaa  <  1:349902/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
CTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAA  >  NC_000913/542924‑543196

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