Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,094 G→C 7.4% intergenic (+61/+167) 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,0940GC 7.4% 119.6 / 3.2 49intergenic (+61/+167)glxK/allEglycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base C (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAG  >  NC_000913/542957‑543232
                                                                                                                                         |                                                                                                                                          
gctgTGAAACCATTCTCTTCAACATCGCGCGAATTCTTGCCTGCGCCATTAAAATAGGTCAGGCAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCa                                                                                                                                       <  1:122886/140‑1 (MQ=255)
 cggttgaaCTAATCGCTTCTACCGCGCGTGAAATCTTGCCTGCGTCATTAGAATAGGTCAGTGAATTATAAACTTATTCTTACCTTTAAAGCGGATGCGATTTATCTCGCTTTAGAGTGCAGTCCTCATGCCGGAGGCGGCAt                                                                                                                                      <  2:771621/137‑1 (MQ=255)
   gttaaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTCACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGCGCCGTACCCATGCCGGAggcggccaga                                                                                                                                    >  2:661400/4‑139 (MQ=255)
   gTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                    >  2:356901/1‑143 (MQ=255)
   gTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                    >  2:356902/1‑143 (MQ=255)
   gTGAAACCAATCTCTTCAACAGCGCGCGAAATAGTGGCTGCGCCATTAAAATAGGTCAGGGAAGTAAAAACTAACCCTTACCTTTAAAGCGGATGAGATTTATATCGCATAAGACCGCCGTACTCATGCCGGCTGCGGCGTGa                                                                                                                                    >  2:887750/1‑143 (MQ=255)
         ccAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCTTGGAATTAAAATCTTACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTActg                                                                                                                              <  1:572558/143‑3 (MQ=255)
          cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCCTGAGTACCat                                                                                                                             >  1:769867/1‑143 (MQ=255)
             tctcTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCAtctc                                                                                                                          >  1:120098/1‑143 (MQ=255)
                 ttCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                                                      >  1:475412/1‑143 (MQ=255)
                        gcgcgcGAAAAACTGCCTGCGCCATTAAAATAGGTCAGCGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGaaaa                                                                                                               <  2:1012480/143‑1 (MQ=255)
                         cgcgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAACCTTCCTGCAAAt                                                                                                              >  2:489096/1‑143 (MQ=255)
                            gcgAAATATTGCCTGCGGCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGTAAATcgc                                                                                                           <  2:907158/143‑1 (MQ=255)
                            gcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGAGGCATGAGGACCATATCCTTCCTGAAAATcgc                                                                                                           >  1:571667/1‑143 (MQ=255)
                                  tatTGCCTGCGCCATTAAAATAGGTCAGGGAATTACAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGAGGCGGCATGAGTACCGTATCCTTCCTGAAAATcggggcagt                                                                                                     >  1:963996/1‑136 (MQ=255)
                                  tatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCCCAAAt                                                                                                     >  1:981136/1‑143 (MQ=255)
                                      gCCTGCGCCATCAAAATAGGTCAGGGAATTAAACACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                 <  1:213978/143‑1 (MQ=255)
                                       ccTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACACTTACACGAAGTACCGAGACTATTTTAGTTGTGGCACATCGCTGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtat                                                                                                >  1:744410/1‑143 (MQ=255)
                                         tGCGCCATTAAAATAGGTCAGGGAAGTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCCAAATCtatca                                                                                              >  2:997540/1‑141 (MQ=255)
                                               aTTAAAATATGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                        <  1:182568/143‑1 (MQ=255)
                                               aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATATATAGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGAACCATATCCTTCCTGAAAATCGCGCAAAATCTATATATTGCa                                                                                        >  2:669375/1‑143 (MQ=255)
                                                 tAAAATAGGTCAGGGCATTAAAAACTAACCCTTACCTCTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCaga                                                                                      <  2:919396/143‑1 (MQ=255)
                                                 tAAAATAGGTCAGGGAGTTAAAAACTAACCCTTACCTTTAAAGCGGGTGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCaga                                                                                      <  2:420899/143‑1 (MQ=255)
                                                 gaagagaGGTGAGGGAATTAAAGACTAACCCTTACCGTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTGCCTGAAAATCGCGCAAATTCTATATATTGCaga                                                                                      <  2:232296/137‑1 (MQ=255)
                                                  aaaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCagag                                                                                     <  1:390297/143‑1 (MQ=255)
                                                               gAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATCTTGCAGAGATCATGTAGGCCt                                                                        >  1:889488/1‑143 (MQ=255)
                                                                       aaCGAACCCTTCCCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcg                                                                <  2:652097/143‑1 (MQ=255)
                                                                         cTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaa                                                              <  1:601255/143‑1 (MQ=255)
                                                                         cTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaa                                                              <  1:48935/143‑1 (MQ=255)
                                                                          tcaCCCGTACCTTCGCAGCGCGTGGGATTGATATCGCATAAGACTGCAGTACTTATGCTGGCTCCGGCATGAGTACCATATCCTTCCTGAAAATTGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaag                                                             <  2:762061/141‑1 (MQ=25)
                                                                             gccTTACCTTTAACGCGGATGCGATTTATGTCGCATAAGCGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAgcgc                                                          <  2:268039/142‑1 (MQ=255)
                                                                             cccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAgcgc                                                          <  2:590056/143‑1 (MQ=255)
                                                                                        aaaGCGGACGGGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATg                                               <  2:591158/143‑1 (MQ=255)
                                                                                            cGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGATAc                                           >  2:943939/1‑143 (MQ=255)
                                                                                              gATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGAGAAGCGCATCAGGCAATGTAAcga                                         >  2:471337/1‑141 (MQ=255)
                                                                                                   gATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAAAGTTACaaaaaaa                                    >  2:802820/1‑143 (MQ=255)
                                                                                                    aTTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCGATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAg                                   >  1:471132/1‑143 (MQ=255)
                                                                                                            cGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGta                           >  2:87632/1‑143 (MQ=255)
                                                                                                             gCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                          >  1:575311/1‑143 (MQ=255)
                                                                                                             gCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                          <  1:489096/143‑1 (MQ=255)
                                                                                                              cATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAACCGAAGCGCATCAGGGAATGTGACCAAAAAAGCTACCGtata                         >  2:57378/1‑143 (MQ=38)
                                                                                                                gaaGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATaaa                       <  1:879397/142‑1 (MQ=255)
                                                                                                                        tggTACTCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATGAACCggggaa               >  1:414817/3‑141 (MQ=37)
                                                                                                                             cTCATGTCGGATGTGGCATGATTATCATATCCTTCCGGAAAATCGGGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCTCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATc          <  2:109577/143‑1 (MQ=37)
                                                                                                                              tCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAACAAAGCCACGGTATAAACCCTGGCAGAATcc         >  2:422251/1‑143 (MQ=37)
                                                                                                                                gtGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa       <  2:194078/142‑1 (MQ=255)
                                                                                                                                aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa       <  2:301301/143‑1 (MQ=255)
                                                                                                                                 tggtggATTCGGCAGGAGTACCAGATCCTCCCTGAAAATGGCGTAAATTCTATATATTGCAGAGATCATGTAGGTCTGATAGGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAATCCGTGGCAAAATCCAAc      <  2:140642/139‑1 (MQ=37)
                                                                                                                                 tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAAc      >  2:380569/1‑143 (MQ=255)
                                                                                                                                     ggATGCGGCATGAGTACCATATCCTTCCTGAACATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACCAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAg  <  2:556892/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
GCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAG  >  NC_000913/542957‑543232

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: