Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,571,720:1 +G 13.5% coding (706/867 nt) 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,7201.G13.5% 119.4 / 16.2 37coding (706/867 nt)yneJLysR family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (16/16);  new base G (2/3);  total (18/19)
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.82e-01

GCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGA  >  NC_002947/3571582‑3571862
                                                                                                                                           |                                                                                                                                              
gcCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTc                                                                                                                                     <  1:218666/150‑1 (MQ=255)
 ccAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCATGCTGCGCGGCAGCATggc                                                                                                                                      >  2:505487/1‑145 (MQ=255)
 ccAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCATGCTGCGCGGCAGCATggc                                                                                                                                      <  1:505487/147‑3 (MQ=255)
  cAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGa                                                                                                                                   <  1:510919/150‑1 (MQ=255)
   aaCTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGAc                                                                                                                                  >  1:314794/1‑150 (MQ=255)
        ttCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATg                                                                                                                                        >  1:1206432/1‑139 (MQ=255)
        ttCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATg                                                                                                                                        <  2:1206432/139‑1 (MQ=255)
           aTACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACAACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGcc                                                                                                                          <  1:884012/150‑1 (MQ=255)
                       cTTCGAAAACGGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAc                                                                                                              <  1:687922/150‑1 (MQ=255)
                        ttCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCAcc                                                                                                             >  2:732464/1‑150 (MQ=255)
                                   gTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                  <  1:732464/150‑1 (MQ=255)
                                   gTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAAGGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATAGGCCCGGGTGCAGCACCGTCAGCGCCTg                                                                                                  >  1:147590/1‑150 (MQ=255)
                                                cAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCgcagc                                                                                                                                           >  1:1003081/1‑96 (MQ=255)
                                                cAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCgcagc                                                                                                                                           <  2:1003081/96‑1 (MQ=255)
                                                      gTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTc                                                                               >  2:55099/1‑150 (MQ=255)
                                                        aCCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCg                                                                             <  1:55099/150‑1 (MQ=255)
                                                            ggCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTAt                                                                         >  2:335210/1‑150 (MQ=255)
                                                              cAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCt                                                                       <  1:1187829/150‑1 (MQ=255)
                                                                  aTTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGaaa                                                                   >  2:128472/1‑150 (MQ=255)
                                                                    tCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTcc                                                                              <  2:1034742/137‑1 (MQ=255)
                                                                    tCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTcc                                                                              >  1:1034742/1‑137 (MQ=255)
                                                                       cGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCgcagca                                                                                                                                          <  1:849447/74‑1 (MQ=255)
                                                                       cGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCgcagca                                                                                                                                          >  2:849447/1‑74 (MQ=255)
                                                                        gAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACctgg                                                             <  2:536283/150‑1 (MQ=255)
                                                                                gTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATggcc                                                                                                                                     >  2:919933/1‑67 (MQ=38)
                                                                                gTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATg                                                     >  2:1259115/1‑150 (MQ=255)
                                                                                gTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCGCGGCAGCATggcc                                                                                                                                     <  1:919933/70‑4 (MQ=37)
                                                                                       cACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGcgacga                                              <  1:128472/150‑1 (MQ=255)
                                                                                         cGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAgg                                            >  1:700681/1‑150 (MQ=255)
                                                                                                          tCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGtt                           <  2:829969/150‑1 (MQ=255)
                                                                                                                cGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTAt                                                                         <  1:949844/98‑1 (MQ=255)
                                                                                                                cGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTAt                                                                         >  2:949844/1‑98 (MQ=255)
                                                                                                                        ggCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTg             <  1:826978/150‑1 (MQ=255)
                                                                                                                          ccTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGaa           <  2:700681/150‑1 (MQ=255)
                                                                                                                           cTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGCACACGTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAg          >  1:728519/1‑150 (MQ=255)
                                                                                                                           cTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAg          >  2:689327/1‑150 (MQ=255)
                                                                                                                                cATGCTGCGCTGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt     >  2:355170/1‑150 (MQ=255)
                                                                                                                                cATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCGAGCTTAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCt     <  1:374676/150‑1 (MQ=255)
                                                                                                                                 aTGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCtt    >  1:1213473/1‑150 (MQ=255)
                                                                                                                                      gcCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTga  <  2:497913/147‑1 (MQ=255)
                                                                                                                                      gcCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTga  >  1:497913/1‑147 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GCCAACTGTTCATACCGCCATCACTTCGAAAACTGGTTCGTGCAGGACCAGGCGGTACCGGGCAAGATTCACGAAATGGAGTCGTACCACGGCATGTTGGCCTGCGTCAGTGCGGGTGCCGGCCTGGCCATGCTGCCGCGCAGCATGCTCGACAATATGCCCGGGTGCAGCACCGTCAGCGCCTGGCCAATGTCGGAGGACTTCCGCTATCTGAAAACCTGGCTGGTATGGCGACGAGGTACGGTCTCGCGCAGTTTGAGCATGTTCGTGAAGTTGCTTGA  >  NC_002947/3571582‑3571862

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