Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,563,842 A→C intergenic (‑239/+80) yjiK ← / ← yjiL SdiA‑regulated family putative membrane‑anchored protein; putative phytase‑like esterase/putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,563,8420AC97.9% 172.1 / NA 47intergenic (‑239/+80)yjiK/yjiLSdiA‑regulated family putative membrane‑anchored protein; putative phytase‑like esterase/putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (21/25);  minor base T (0/1);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAG  >  NC_000913/4563695‑4563989
                                                                                                                                                   |                                                                                                                                                   
aaTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGcaaa                                                                                                                                                  <  2:821032/151‑1 (MQ=255)
   gCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTc                                                                                                                                               <  2:42151/151‑1 (MQ=255)
     aGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCAt                                                                                                                                             >  2:279647/1‑151 (MQ=255)
          cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCg                                                                                                                                        >  2:467867/1‑151 (MQ=255)
           aGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGc                                                                                                                                       >  2:63182/1‑151 (MQ=255)
                       gctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAAc                                                                                                                             <  2:779770/149‑1 (MQ=255)
                          ggCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCt                                                                                                                          <  2:65043/149‑1 (MQ=255)
                            cAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGc                                                                                                                      <  1:275520/151‑1 (MQ=255)
                                        tcagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATg                                                                                                           >  1:427840/1‑150 (MQ=255)
                                         cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGAc                                                                                                         <  1:86710/151‑1 (MQ=255)
                                          agtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACg                                                                                                        >  2:22018/1‑151 (MQ=255)
                                            tcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGcc                                                                                                      >  1:360068/1‑151 (MQ=255)
                                            tcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGcc                                                                                                      >  1:107813/1‑151 (MQ=255)
                                             cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGcc                                                                                                      <  1:424564/150‑1 (MQ=255)
                                             cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGcc                                                                                                      <  2:257204/150‑1 (MQ=255)
                                                cccAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGAcc                                                                                                  <  2:756478/151‑1 (MQ=255)
                                                  cAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGt                                                                                                <  1:42827/151‑1 (MQ=255)
                                                  cAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGt                                                                                                >  2:901959/1‑151 (MQ=255)
                                                   aGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTg                                                                                               <  2:600631/151‑1 (MQ=255)
                                                   aGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTg                                                                                               >  1:76379/1‑151 (MQ=255)
                                                        tCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTg                                                                                          >  1:482414/1‑151 (MQ=255)
                                                         cACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTg                                                                                          >  1:652140/1‑150 (MQ=255)
                                                            ccATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaat                                                                                      >  2:19137/1‑151 (MQ=255)
                                                                tGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaataaa                                                                                   >  1:403153/1‑150 (MQ=255)
                                                                 gTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaataaaa                                                                                  >  2:427625/1‑150 (MQ=255)
                                                                   aTTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaataaaa                                                                                  >  2:256496/1‑148 (MQ=255)
                                                                       ccTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAataaaaataaaa                                                                            <  1:828791/150‑1 (MQ=255)
                                                                       ccTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAAtaaaaataaaaa                                                                           <  1:782702/151‑1 (MQ=255)
                                                                       ccTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAAtaaaaataaaaa                                                                           <  1:430081/151‑1 (MQ=255)
                                                                           gCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAActct                                                                       <  1:190622/151‑1 (MQ=255)
                                                                           gCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAActct                                                                       <  2:785466/151‑1 (MQ=255)
                                                                            ccGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCtt                                                                      <  2:380193/151‑1 (MQ=255)
                                                                             cGGGGAGGCAGCTCAGCCAAATGCATACCCCTCGTCTGCCACGTGATAAACCGAGAATATCGCCGCCTCGTAAAGTCATGCCAGCTTTGTGAAAAGCGTTATCGGTTTGTATGACGGCGACCGTGGAGTGAaataaaaataaatccccttt                                                                     <  1:622085/151‑9 (MQ=255)
                                                                               gggggACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCa                                                                   <  2:287352/151‑1 (MQ=255)
                                                                                     cAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCg                                                             >  1:454896/1‑151 (MQ=255)
                                                                                     cAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCg                                                             <  2:60473/151‑1 (MQ=255)
                                                                                      aGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCgt                                                            >  2:356949/1‑151 (MQ=255)
                                                                                            gCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTc                                                      <  1:578061/151‑1 (MQ=255)
                                                                                                  tGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGt                                                >  2:809922/1‑151 (MQ=255)
                                                                                                           cGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGTCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAAtt                                       <  1:803415/151‑1 (MQ=255)
                                                                                                               gccgccGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAcc                                    >  1:806363/1‑150 (MQ=255)
                                                                                                                  gccgcGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAccgccg                                <  2:614150/151‑1 (MQ=255)
                                                                                                                     gcgTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAAATAacgcccccc                             >  1:842824/1‑147 (MQ=255)
                                                                                                                          aaTCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAAtt                        >  1:88819/1‑151 (MQ=255)
                                                                                                                             cACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTgc                      <  1:365421/150‑1 (MQ=255)
                                                                                                                                         cGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTgc         <  2:614256/151‑1 (MQ=255)
                                                                                                                                         cGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTgc         <  2:44542/151‑1 (MQ=255)
                                                                                                                                          gCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTgcg        <  1:266277/151‑1 (MQ=255)
                                                                                                                                                gcgcAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAg  <  2:762265/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAG  >  NC_000913/4563695‑4563989

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