Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,222,944 A→G 80.4% intergenic (‑251/+309) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,222,9440AG80.4% 103.3 / 18.4 51intergenic (‑251/+309)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base G (20/21);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GCGCCCCTCCGTGCCCCGTTGCGAATTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATG  >  NC_002947/3222808‑3223083
                                                                                                                                        |                                                                                                                                              
gcgcCCCTCCGTGCCCCGTTGCGAATTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGc                                                                                                                                    <  1:162721/149‑1 (MQ=255)
    cccTCCGTGCCCCGTTGCGAATTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGc                                                                                                                                >  5:124496/1‑149 (MQ=255)
                     cgtaatcgccgccCGATTGCGTGGAGGGTCCGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTCATGTCGTCAATTTCAGGCAACCAGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAgcgc                                                                                                            <  6:83216/142‑1 (MQ=255)
                      gAATTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGcc                                                                                                           >  1:71639/1‑149 (MQ=255)
                        aTTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGcc                                                                                                                                             >  6:7358/1‑116 (MQ=255)
                        aTTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGcc                                                                                                                                             <  5:7358/116‑1 (MQ=255)
                          tGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGAGCCTTgg                                                                                                       >  7:83503/1‑149 (MQ=255)
                            gccgccCGATTGCGTGGAGTGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGt                                                                                                      <  4:170016/148‑1 (MQ=255)
                             ccgccCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTTGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGc                                                                                                    <  4:39821/149‑1 (MQ=255)
                              cgccCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCt                                                                                                   <  1:116975/149‑1 (MQ=255)
                              cgccCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGGGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGc                                                                                                    >  1:69523/1‑148 (MQ=255)
                                     ttGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGcc                                                                                             >  2:267544/1‑148 (MQ=255)
                                          tGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCAtttt                                                                                       >  6:86175/1‑149 (MQ=255)
                                            gAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGAGCCTTGGTGCTAAAGCCCATTTTCg                                                                                     <  8:83503/149‑1 (MQ=255)
                                            gAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGATCGGAAATTTCCGGCAACCTGCGTCGTTGGGGAATAAGGATCCCGTGGCCCTTGCATTCCTGCTAGGCCCGCG‑‑ACAGCGGCTTGGTGGTAAAGCCCAttttag                                                                                     >  1:235259/1‑147 (MQ=255)
                                                   aGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGc                                                                              <  2:71639/149‑1 (MQ=255)
                                                             cAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTtg                                                                    >  7:169960/1‑149 (MQ=255)
                                                             cAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTtg                                                                    >  3:164301/1‑149 (MQ=255)
                                                                     gagaGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTAt                                                            <  2:69523/149‑1 (MQ=255)
                                                                      agagGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACt                                                                          >  2:277637/1‑134 (MQ=255)
                                                                      agagGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACt                                                                          <  1:277637/134‑1 (MQ=255)
                                                                      agagGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTAtt                                                           >  4:134868/1‑149 (MQ=255)
                                                                         ggTGGATGAAGGTGAGGTCGTCAATTTCAGGCTACCTGCGTCTTTGTGTAATAAGGATCCGGTGGACCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTAttttt                                                        <  8:169960/149‑3 (MQ=255)
                                                                          gTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGtt                                                       <  5:222415/149‑1 (MQ=255)
                                                                          gTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGt                                                        >  7:30009/1‑148 (MQ=255)
                                                                          gTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGt                                                        >  7:237787/1‑148 (MQ=255)
                                                                          gTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGt                                                        >  2:99488/1‑148 (MQ=255)
                                                                           tGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTta                                                      >  3:177053/1‑149 (MQ=255)
                                                                             gaggaAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTtat                                                     <  5:86175/148‑1 (MQ=255)
                                                                                      gtggtggtCAATTTCAGGGAACCTGCGTCTTTTTTGAATAATGTTCCGGTGGGCCTTGCATGCCTGCTCGGAACGCGCTCCGAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAgt                                              <  1:267544/147‑1 (MQ=255)
                                                                                         ggcgtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGGGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGtttt                                         <  2:235259/146‑1 (MQ=255)
                                                                                          tcgtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTAttt                                                          >  8:38091/1‑130 (MQ=255)
                                                                                          tcgtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTAttt                                                          <  7:38091/130‑1 (MQ=255)
                                                                                          tcgtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTg                                        >  8:254851/1‑148 (MQ=255)
                                                                                          tcgtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGt                                       <  3:134868/149‑1 (MQ=255)
                                                                                            gtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTc                                                                                  <  6:134276/104‑1 (MQ=255)
                                                                                            gtcAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTc                                                                                  >  5:134276/1‑104 (MQ=255)
                                                                                               aaTTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTaaaa                                   >  6:127622/1‑148 (MQ=255)
                                                                                                      ggCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAg                              >  1:323653/1‑146 (MQ=255)
                                                                                                      ggCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGa                            <  2:323653/148‑1 (MQ=255)
                                                                                                      ggCAACCTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAg                           <  6:124496/149‑1 (MQ=255)
                                                                                                        cAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGtt                         <  8:237787/149‑1 (MQ=255)
                                                                                                            cTGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGc                     <  8:30009/149‑1 (MQ=255)
                                                                                                             tGCGTCTTTGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCt                    >  6:57281/1‑149 (MQ=255)
                                                                                                               cGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCtt                   <  4:164301/148‑1 (MQ=255)
                                                                                                                   tttGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGtttt                                         >  7:255884/1‑122 (MQ=255)
                                                                                                                   tttGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGtttt                                         <  8:255884/122‑1 (MQ=255)
                                                                                                                   gttGTGGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACAt               <  4:177053/147‑1 (MQ=255)
                                                                                                                     tgtgGAATAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGc             <  7:254851/148‑1 (MQ=255)
                                                                                                                            tAAGGATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTGTATTTGTTATATCAAGTGTGTTGTAAAATACAGGAGTTTTGCTTACATGCAAGAAGc      >  7:288792/1‑148 (MQ=255)
                                                                                                                               ggATCCGGTGGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATAAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATg  >  7:280108/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                              
GCGCCCCTCCGTGCCCCGTTGCGAATTGGCCGCCCGATTGCGTGGAGGGTCAGGCCGAGGCCAGGTTTGGAGAGGTGGAGGAAGGTGAGGTCGTCAATTTCAGGCAACCTGCGTCTTTGTGGAATAAGGATCCGGTAGCCCTGGCATGCCTGCTCGGACCGCG‑‑ACAGCGCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATG  >  NC_002947/3222808‑3223083

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