Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,029 G→C 8.3% intergenic (+11/‑165) PP_5259 → / → ydcJ LysR family transcriptional regulator/metalloprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,0290GC 8.3% 136.4 / 5.0 49intergenic (+11/‑165)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base C (2/2);  total (25/24)
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

CGACGGTAGCCTGGTGGTGCCGTGGAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAC  >  NC_002947/6004880‑6005175
                                                                                                                                                     |                                                                                                                                                  
cGACGGTAGCCTGGTGGTGCCGTGGAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTcg                                                                                                                                                    >  2:751734/1‑150 (MQ=255)
 gACGGTAGCCTGGTGGTGCCGTGGAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTcgc                                                                                                                                                   <  1:611756/150‑1 (MQ=255)
         ccTGGTGGTGCCGTGGAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGCTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCTGCTAAAg                                                                                                                                           <  2:710482/150‑1 (MQ=255)
                        gAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGa                                                                                                                            >  1:481187/1‑150 (MQ=255)
                         aaCCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAc                                                                                                                           >  1:553215/1‑150 (MQ=255)
                            cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                        <  1:125481/150‑1 (MQ=255)
                             aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgcg                                                                                                                       >  1:241155/1‑150 (MQ=255)
                             aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgcg                                                                                                                       <  2:126630/150‑1 (MQ=255)
                                  aTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGTAgga                                                                                                                  <  1:751734/150‑2 (MQ=255)
                                       cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTGCGGGGGCCGGGTAGGCGCGGc                                                                                                             >  1:12937/1‑150 (MQ=255)
                                              ggAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGcc                                                                                                      >  1:270134/1‑150 (MQ=255)
                                                     aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGTTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGGAGGCGCGGCTTTAGCAGGGAAAg                                                                                               >  1:731310/1‑150 (MQ=255)
                                                     aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                               >  1:770343/1‑150 (MQ=255)
                                                     aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTCAGCAAGGCATTCGGTTTCTACGGAGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                               >  2:20559/1‑150 (MQ=255)
                                                               gTATGCCGAGCATGCGGCAGAGGTGCCAAAGGGGCGGGCGCTGGTGGAATGGATTGGCGCGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGCAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                     <  2:600644/150‑1 (MQ=255)
                                                                  tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                             >  2:224015/1‑139 (MQ=255)
                                                                  tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                             <  1:224015/139‑1 (MQ=255)
                                                                            gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaaaa                                                                        >  1:89756/1‑150 (MQ=255)
                                                                               gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGCTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaaga                                                                                               <  1:334545/124‑3 (MQ=14)
                                                                               gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGCTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaaga                                                                                               >  2:334545/1‑122 (MQ=14)
                                                                                cAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                    <  2:481187/150‑1 (MQ=255)
                                                                                       gCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTa                                                                                                                     >  2:551723/1‑94 (MQ=25)
                                                                                       gCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTa                                                                                                                     <  1:551723/94‑1 (MQ=25)
                                                                                        ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCATACGCGCGGCGTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGacacac                                                            >  1:776729/1‑150 (MQ=255)
                                                                                              gTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGCCACACTCACtt                                                      <  2:715106/150‑1 (MQ=255)
                                                                                                 ggggCGCTGGAGGAATGGATTGGCGAGCAGTTGCAGGTTTGCCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttt                                                   <  2:770343/149‑1 (MQ=255)
                                                                                                  gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCtt                                                                                                           >  2:274497/1‑93 (MQ=255)
                                                                                                  gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCtt                                                                                                           <  1:274497/93‑1 (MQ=255)
                                                                                                  gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtttttt                                                  >  1:333085/1‑150 (MQ=255)
                                                                                                            ggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                        <  2:553215/150‑1 (MQ=255)
                                                                                                             gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAg                                       <  2:241155/150‑1 (MQ=255)
                                                                                                               aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCTCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCa                                     >  2:495099/1‑150 (MQ=255)
                                                                                                                 ggATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCaaa                                   <  1:299367/150‑1 (MQ=255)
                                                                                                                   aTTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                             >  1:374116/1‑74 (MQ=255)
                                                                                                                   aTTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                             <  2:374116/74‑1 (MQ=255)
                                                                                                                        tgAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGATAAAGCGGCTCCTACGCAGACCGCGTAGGCGCAACTTTAGCCGGGAAAGGGCCGTCGAGGTATATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAtt                            <  2:12937/149‑1 (MQ=255)
                                                                                                                        cGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAGCGCGTAGGCGCGGCGTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTATCTTGTGAGGAAAGTGCagg                            >  1:292762/1‑148 (MQ=255)
                                                                                                                            cAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTGGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCa                        >  1:377516/1‑150 (MQ=255)
                                                                                                                                   aGGTGTGAAGGCCTCTTGGCGGATAAAGCCGCTCTTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCGGTCGAGGTTTATCAAAAAAGGAATGACACACTCACTTTTTTTCGTTTGTGAGCAAAGTGCATTCCCAGCTTTCa                 <  2:331303/150‑1 (MQ=255)
                                                                                                                                     gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATg               >  1:611489/1‑150 (MQ=255)
                                                                                                                                      tttGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTTCGGTTAAAGTGCATTCCCAGCTTTCATGt              <  1:20559/150‑1 (MQ=255)
                                                                                                                                          ccgatctcctcGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGTCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGc          <  2:420016/145‑1 (MQ=255)
                                                                                                                                          aCGGCCTCTTCGCGGGTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTTGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAgg          >  2:595017/1‑149 (MQ=255)
                                                                                                                                             gCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTc       >  1:54881/1‑150 (MQ=255)
                                                                                                                                             gCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTc       <  2:776729/150‑1 (MQ=255)
                                                                                                                                              ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                    <  1:640051/72‑1 (MQ=255)
                                                                                                                                              ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                    >  2:640051/1‑72 (MQ=255)
                                                                                                                                                 ctTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGaa   <  2:333085/150‑1 (MQ=255)
                                                                                                                                                  ttCGCGGCTAAAGCCGCTCCTACGGGGTCCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCgagatcgg  >  1:420016/1‑143 (MQ=255)
                                                                                                                                                  ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                       >  2:478011/1‑81 (MQ=255)
                                                                                                                                                  ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                       <  1:478011/81‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGACGGTAGCCTGGTGGTGCCGTGGAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAC  >  NC_002947/6004880‑6005175

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