Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,571,643 C→T 9.5% C224C (TGC→TGT)  yneJ → LysR family DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,571,6430CT 9.5% 114.7 / 3.5 42C224C (TGC→TGT) yneJLysR family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (2/2);  total (24/18)
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.56e-01

GGACGTGAACGGTGAAAGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCA  >  NC_002947/3571505‑3571791
                                                                                                                                           |                                                                                                                                                    
ggACGTGAACGGTGAAAGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCgg                                                                                                                                            <  6:139676/149‑1 (MQ=255)
    gTGAACGGTGAAAGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCACGCGGTACAGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGc                                                                                                                                        >  5:220/1‑149 (MQ=255)
        aCGGTGAAAGCATCTACTCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGc                                                                                                                                    <  3:195076/149‑1 (MQ=255)
            tGAAAGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCCGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCGGCGTCAGTGCGGGTGCCGGACTgg                                                                                                                                >  4:266587/1‑149 (MQ=255)
              aaaGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTGAGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGc                                                                                                                               >  1:290330/1‑149 (MQ=255)
               aaGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCa                                                                                                                             <  2:71068/149‑1 (MQ=255)
                 gCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCAtg                                                                                                                           <  2:244063/149‑1 (MQ=255)
                     cTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCAtgctgc                                                                                                                       >  2:224819/1‑149 (MQ=255)
                        cGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCAtgctgc                                                                                                                       <  4:71609/146‑1 (MQ=255)
                           cTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGTGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgcg                                                                                                                   >  7:12900/1‑147 (MQ=255)
                           cTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGTGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCgca                                                                                                                 >  7:345929/1‑149 (MQ=255)
                           cTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgcg                                                                                                                   >  6:218102/1‑147 (MQ=255)
                           cTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgcg                                                                                                                   >  5:161406/1‑147 (MQ=255)
                           cTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCcgcg                                                                                                                   >  2:343275/1‑147 (MQ=255)
                                   ccAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATggcc                                                                                                         <  1:224819/149‑4 (MQ=255)
                                     aaCTGTTCATACCGCCATCACTTCGAAAAGTGTTT‑CGTGCAGGACCAGGCGGTACCGGTCTAGATTCACTATATGGATTCGTACCACGTCATGTTGGCGTGCGTCAGTGCGTGTGCCGGc                                                                                                                                    <  7:58315/120‑1 (MQ=255)
                                          ttCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACaa                                                                                                    <  2:290330/147‑1 (MQ=255)
                                            cATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGTGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAatat                                                                                                 <  8:345929/148‑1 (MQ=255)
                                            cATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGTGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAatat                                                                                                 <  8:12900/148‑1 (MQ=255)
                                                  gCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGTCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCggg                                                                                          >  6:215405/1‑149 (MQ=255)
                                                          ttCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTc                                                                                                         <  4:185687/126‑1 (MQ=255)
                                                          ttCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTc                                                                                                         >  3:185687/1‑126 (MQ=255)
                                                          ttCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAc                                                                                  >  4:65984/1‑149 (MQ=255)
                                                             gaaacgtGGTT‑CGTGCTGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGt                                                                               <  6:220/143‑1 (MQ=255)
                                                             gAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGt                                                                               >  4:212600/1‑149 (MQ=255)
                                                                              aGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATg                                                               <  6:161406/149‑1 (MQ=255)
                                                                                   cAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                          <  3:212600/149‑1 (MQ=255)
                                                                                   cAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCgga                                                          <  1:343275/149‑1 (MQ=255)
                                                                                   cAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCg                                                            >  8:224893/1‑147 (MQ=255)
                                                                                       cGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGAc                                                      <  5:166489/149‑1 (MQ=255)
                                                                                         gTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACtt                                                    >  1:142446/1‑149 (MQ=255)
                                                                                                  aaGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGTGCCTGGCCAATGTCGGAGGACTTCCGCTATCt                                           >  1:370138/1‑149 (MQ=255)
                                                                                                       tCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGaaaa                                      >  3:143370/1‑149 (MQ=255)
                                                                                                         aCGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAAcc                                    >  3:360992/1‑149 (MQ=255)
                                                                                                           gAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctg                                  >  6:87215/1‑149 (MQ=255)
                                                                                                                 gAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCCTGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctggctgg                             >  2:15160/1‑148 (MQ=255)
                                                                                                                 gAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctggctgg                             <  2:142446/148‑1 (MQ=255)
                                                                                                                                aTGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTAc             <  3:65984/149‑1 (MQ=255)
                                                                                                                                  gTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACgg           >  7:269833/1‑149 (MQ=255)
                                                                                                                                       ccTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTcg      >  2:291117/1‑149 (MQ=255)
                                                                                                                                        cTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGAGTTCCGCTAGCTGAAAACCTGGCGGGTATGGCGACGAGGTACGGTCgcgc     >  5:258/1‑149 (MQ=255)
                                                                                                                                        cTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTcgc     >  1:253276/1‑149 (MQ=255)
                                                                                                                                           cGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCgcgc   <  5:215405/148‑1 (MQ=255)
                                                                                                                                           cGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCa  <  8:269833/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GGACGTGAACGGTGAAAGCATCTACGCCTTCCGCGCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTT‑CGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCA  >  NC_002947/3571505‑3571791

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