Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,127,940 C→A intergenic (‑46/+138) yiiX ← / ← metJ putative lipid binding hydrolase/DNA‑binding transcriptional repressor MetJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,9400CA95.9% 168.3 / ‑2.0 49intergenic (‑46/+138)yiiX/metJputative lipid binding hydrolase/DNA‑binding transcriptional repressor MetJ
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (21/26);  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 = 5.87e-01

TGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCG  >  NC_000913/4127809‑4128062
                                                                                                                                   |                                                                                                                          
tGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCgttgtt                                                                                                                    >  1:630999/1‑140 (MQ=255)
       tgatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCg                                                                                                             >  2:604707/1‑140 (MQ=255)
       tgatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCg                                                                                                             >  2:403736/1‑140 (MQ=255)
         atgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGa                                                                                                           <  1:155995/140‑1 (MQ=255)
         atgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAACAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGa                                                                                                           >  2:24511/1‑140 (MQ=255)
            tCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATgc                                                                                                        <  1:397761/140‑1 (MQ=255)
            tCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATgc                                                                                                        >  2:54866/1‑140 (MQ=255)
             cGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATgcg                                                                                                       <  2:769011/140‑1 (MQ=255)
              gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATgcgg                                                                                                      >  2:164833/1‑140 (MQ=255)
                  gTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTg                                                                                                  <  2:363711/140‑1 (MQ=255)
                  gTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTg                                                                                                  <  1:391001/140‑1 (MQ=255)
                        ggCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                            <  2:674781/140‑1 (MQ=255)
                             ccAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATc                                                                                       <  1:4498/140‑1 (MQ=255)
                                          ggTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGtt                                                                          <  1:164833/140‑1 (MQ=255)
                                                 aaaCCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                     >  1:684055/1‑122 (MQ=255)
                                                 aaaCCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                     <  2:684055/122‑1 (MQ=255)
                                                      agcagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCa                                                              >  2:252559/1‑140 (MQ=255)
                                                      agcagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCa                                                              >  1:38765/1‑140 (MQ=255)
                                                        cagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAAt                                                            >  2:790662/1‑140 (MQ=255)
                                                        cagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAAt                                                            <  2:630999/140‑1 (MQ=255)
                                                        cagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAAt                                                            <  1:54866/140‑1 (MQ=255)
                                                         agGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATa                                                           <  1:125321/140‑1 (MQ=255)
                                                           gCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATaaa                                                         <  1:370912/140‑1 (MQ=255)
                                                           gCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATaaa                                                         <  1:403736/140‑1 (MQ=255)
                                                             tGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAAtt                                                       >  1:420757/1‑140 (MQ=255)
                                                               aGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGc                                                     >  1:144531/1‑140 (MQ=255)
                                                                     agcagcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATg                                               >  2:715010/1‑140 (MQ=255)
                                                                        agcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACg                                            <  1:604707/140‑1 (MQ=255)
                                                                        agcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACg                                            <  1:252559/140‑1 (MQ=255)
                                                                          ccTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTa                                          <  1:715010/140‑1 (MQ=255)
                                                                            tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAgg                                        <  1:24511/140‑1 (MQ=255)
                                                                            tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAgg                                        >  1:251475/1‑140 (MQ=255)
                                                                            tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAgg                                        <  2:38765/140‑1 (MQ=255)
                                                                             gTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGc                                       <  2:420757/140‑1 (MQ=255)
                                                                                        ttAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCt                                     <  1:3784/131‑1 (MQ=255)
                                                                                        ttAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCt                                     >  2:3784/1‑131 (MQ=255)
                                                                                        ttAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGt                            >  1:242479/1‑140 (MQ=255)
                                                                                        ttAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGt                            >  1:123303/1‑140 (MQ=255)
                                                                                               tCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCAt                     <  2:144531/140‑1 (MQ=255)
                                                                                               tCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCAt                     <  1:790662/140‑1 (MQ=255)
                                                                                                 cATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCa                   >  2:619511/1‑140 (MQ=255)
                                                                                                 cATGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCa                   >  1:421280/1‑140 (MQ=255)
                                                                                                  aTGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAg                  >  1:621865/1‑140 (MQ=255)
                                                                                                   tGTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAgg                 <  2:151400/140‑1 (MQ=255)
                                                                                                    gTATAAAAAGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGc                <  2:123303/140‑1 (MQ=255)
                                                                                                         aaaaGCGGTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTc           >  1:215078/1‑140 (MQ=255)
                                                                                                              cGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTc      >  1:329376/1‑140 (MQ=255)
                                                                                                                gTGGGTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcg    <  2:421280/140‑1 (MQ=255)
                                                                                                                  gggTCGCAGACAACGTGATCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgcg  <  2:370443/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                          
TGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCG  >  NC_000913/4127809‑4128062

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