Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,006,587 A→C 16.8% intergenic (+55/+44) ydcJ → / ← PP_5261 metalloprotein/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,006,5870AC16.8% 85.6 / 3.0 39intergenic (+55/+44)ydcJ/PP_5261metalloprotein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (19/13);  new base C (3/4);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGC  >  NC_002947/6006439‑6006734
                                                                                                                                                    |                                                                                                                                                   
gTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGa                                                                                                                                                     <  1:46227/149‑1 (MQ=255)
       gTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTgg                                                                                                                                              <  1:131038/149‑1 (MQ=255)
       gTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTgg                                                                                                                                              <  5:314792/149‑1 (MQ=255)
            ggATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCg                                                                                                                                         <  2:219191/149‑1 (MQ=255)
              aTGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGGAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCggg                                                                                                                                       <  5:321362/149‑1 (MQ=255)
                            tATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCGGGGCACAGGCGTTGACCCGGGGTTCGATGTAAATCTGTACCGGCCTCTGCGCGGGCGCGCCCGCTCCCACAGGTACAGGGCATACGCCCTGTGGGAGCGGGCGAGCCCGCGaa                                                                                                                           <  7:269294/147‑1 (MQ=38)
                                 gcAATGCGAGGGCCGGTCGCTCCTGGCCTGTGCCCCGGCGGTTAACGTGGGTGCCATTTAAATCGGTCCCGGCCTCTTCGGGGGCGCGCCCGCTCCCCCCGGTACTGGGCATCCGACCTGTGGGAGGGGGCGAGCCCGCGAAGAGGCCg                                                                                                                    <  8:109213/149‑1 (MQ=19)
                                                     tGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCGCCCCCAGGGACAGGGGATACGCCCTGTGGGAGCGGGCGGGCCCGCGGAGAGGCCGGTGCAGGTGATGAAAGAcgc                                                                                                >  4:46049/1‑149 (MQ=38)
                                                      gCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCCCTACCCCAGGTAAAGGGGATACGACCTAGGGGAGCGGGCGCCCCCGCGAAGAGGCCGGGGGAGGGGCTGAAAGAcccg                                                                                               >  7:315792/1‑149 (MQ=255)
                                                        aGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTGCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGc                                                                                                                                      >  3:240451/1‑108 (MQ=255)
                                                        aGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGc                                                                                                                                      <  4:240451/108‑1 (MQ=255)
                                                                 cacaGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCCCAGGTACAGGGCCTACGCCCTGTGGGAGCGGGCGAGGCCGCGAAGAGGCCGGTGCAGGGGATCAAAGCCGCGGCAGGATTTcc                                                                                    >  6:367980/1‑149 (MQ=38)
                                                                              ccTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCACACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAAt                                                                       <  2:62938/149‑1 (MQ=255)
                                                                                           tAAATCTGAACCGCCCTCTTCGCGGGCCCGCCCGCTCCCACCGGAACCGGGCAACCCACCTGGGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCAc                                                          <  3:46049/149‑1 (MQ=37)
                                                                                                             ttCGCGGGCTCGCCCGCTCCCCCAGGTAGAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGGAAATCCATTGCCGCCACCGCGCCGTCCCGGGACGc                                        >  2:230772/1‑149 (MQ=38)
                                                                                                             ttCGCGGGCTCGCCCGCTCCCACCGGTACAGGGCATACCACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAAc                                          <  7:119980/147‑1 (MQ=255)
                                                                                                             ttCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAAc                                          <  7:240355/147‑1 (MQ=255)
                                                                                                             tccgcgGGCCCCCCCGCACCCACAGGAACAGGCCATCCCCCCTGGGGCTGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAAc                                          <  3:39127/145‑1 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAAccgccg                            >  3:50096/1‑148 (MQ=255)
                                                                                                                            ggacccACTGGAACCGGGCATACCCCCTGGGGGAGGGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCaa                         <  7:228262/146‑1 (MQ=255)
                                                                                                                            gCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCaa                         >  5:224561/1‑149 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTc                                                 <  4:196170/121‑1 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTc                                                 >  3:196170/1‑121 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAAc                                 <  5:178208/137‑1 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAAc                                 >  6:178208/1‑137 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCg                      >  4:298604/1‑148 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCg                      >  2:333869/1‑148 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGc                     >  1:263442/1‑149 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGc                     >  6:213832/1‑149 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGc                     >  3:41879/1‑149 (MQ=255)
                                                                                                                                ccACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGc                     >  3:219028/1‑149 (MQ=255)
                                                                                                                                  acaGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTc                   >  3:85865/1‑149 (MQ=255)
                                                                                                                                   caGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTcg                  >  8:280705/1‑149 (MQ=255)
                                                                                                                                   caGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGAAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTcg                  >  3:139100/1‑149 (MQ=255)
                                                                                                                                    tcccaCAGGTCGTATGCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAAccgc                              >  7:221144/5‑136 (MQ=39)
                                                                                                                                    tcccaCAGGTCGTATGCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAAccgc                              <  8:221144/132‑1 (MQ=39)
                                                                                                                                          aGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTc           <  1:333869/149‑1 (MQ=255)
                                                                                                                                              cATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCg       >  3:101401/1‑149 (MQ=255)
                                                                                                                                               cctgtaCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCgg      <  6:200497/144‑1 (MQ=255)
                                                                                                                                                tACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCggg     >  6:130836/1‑149 (MQ=255)
                                                                                                                                                   gACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTgc  >  8:292275/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGCCGCAACCGCTCGCCGCCTCGCCGGGTGC  >  NC_002947/6006439‑6006734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑