Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,077 G→C 8.5% intergenic (+59/‑117) 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,0770GC 8.5% 141.8 / 5.4 47intergenic (+59/‑117)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (2/2);  total (23/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 = 6.19e-01

ATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGACAGC  >  NC_002947/6004933‑6005226
                                                                                                                                                |                                                                                                                                                     
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                                                                                                                                                <  1:224015/139‑1 (MQ=255)
             tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                                                >  2:224015/1‑139 (MQ=255)
                       gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaaaa                                                                                                                           >  1:89756/1‑150 (MQ=255)
                          gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGCTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaaga                                                                                                                                                  >  2:334545/1‑122 (MQ=14)
                          gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGCTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaaga                                                                                                                                                  <  1:334545/124‑3 (MQ=14)
                           cAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                                                                       <  2:481187/150‑1 (MQ=255)
                                   ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCATACGCGCGGCGTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGacacac                                                                                                               >  1:776729/1‑150 (MQ=255)
                                         gTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGCCACACTCACtt                                                                                                         <  2:715106/150‑1 (MQ=255)
                                            ggggCGCTGGAGGAATGGATTGGCGAGCAGTTGCAGGTTTGCCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttt                                                                                                      <  2:770343/149‑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)
                                                                   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                                                                                                                                       >  2:640051/1‑72 (MQ=255)
                                                                                         ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                                                                       <  1:640051/72‑1 (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)
                                                                                                          ccGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCTAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGc                                        >  2:311682/1‑150 (MQ=38)
                                                                                                            gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                                      >  1:422578/4‑150 (MQ=39)
                                                                                                            gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                                           <  1:630743/94‑1 (MQ=17)
                                                                                                            gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                                           >  2:630743/4‑97 (MQ=17)
                                                                                                                cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGa                                  >  1:426637/1‑150 (MQ=255)
                                                                                                                  aCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATg                                >  1:456562/1‑150 (MQ=255)
                                                                                                                     gggCCCGCGTGCGCGCTGCTTCAGGCGGGAAAGGGACGTAGAGGTAGATCAAAAAAGTAATGACACACTCACTTTTTTTCGCTTGTGATCAAAGTGTATGCCCAGGTTTCATGTAAGCGTCCGACCCCAACCAGGAAGCTAGCGATGccc                             <  2:292762/150‑1 (MQ=255)
                                                                                                                                 gagcgGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGACCAATAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTAAAGTCCGCGTGGGAGACCGACCATGGGTCTAGCGATGCCCGCCAACGATTTc                 <  1:698857/148‑1 (MQ=255)
                                                                                                                                   gcgGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGt               >  2:831/1‑150 (MQ=255)
                                                                                                                                        tttAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGcc          <  2:377516/150‑1 (MQ=255)
                                                                                                                                        tttAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGcc          <  2:731310/150‑1 (MQ=255)
                                                                                                                                        ttgagacGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGcc          <  2:270134/144‑1 (MQ=255)
                                                                                                                                            gCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGa      <  2:611489/150‑1 (MQ=255)
                                                                                                                                                ggAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGACAGc  <  1:746993/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGACAGC  >  NC_002947/6004933‑6005226

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: