Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 381,856 T→A 10.7% D230V (GAC→GTC)  PP_0317 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947381,8560TA10.7% 144.9 / 8.9 56D230V (GAC→GTC) PP_0317methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (3/3);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

GCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTGG  >  NC_002947/381715‑382000
                                                                                                                                             |                                                                                                                                                
gccgccGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGAc                                                                                                                                           >  3:195058/1‑149 (MQ=255)
 ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCAAGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCTGCGACt                                                                                                                                          >  6:67653/1‑149 (MQ=255)
  cgccgAGCCCTGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGAcct                                                                                                                                         <  4:180039/149‑3 (MQ=255)
           cAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGgcg                                                                                                                                 >  1:138388/1‑148 (MQ=255)
            aGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCa                                                                                                                               >  6:130452/1‑149 (MQ=255)
              cTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAtt                                                                                                                             >  3:249073/1‑149 (MQ=255)
                 ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                            >  2:46208/1‑147 (MQ=255)
                     gCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGc                                                                                                                      >  6:65625/1‑149 (MQ=255)
                     gCGGAGCCGGCGATCAGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGc                                                                                                                      <  8:225879/149‑1 (MQ=255)
                       ggAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGctt                                                                                                                    <  7:138063/149‑1 (MQ=255)
                                gATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGcttct                                                                                                                  >  8:188714/1‑142 (MQ=255)
                                    tgctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCtt                                                                                                       >  1:82595/1‑149 (MQ=255)
                                     gctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCtt                                                                                                       >  5:121246/1‑148 (MQ=255)
                                      ctgcAAGGTATCGCGCAGGTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGTCTGCTTCTTCGGTGCCCTCTGCGGTGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGGTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                                     <  5:67653/149‑1 (MQ=255)
                                      ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                                     >  4:146546/1‑149 (MQ=255)
                                      ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                                     <  4:249073/149‑1 (MQ=255)
                                      ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                                     <  5:130452/149‑1 (MQ=255)
                                      ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                                     >  3:214739/1‑149 (MQ=255)
                                      ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                                     <  7:188714/149‑1 (MQ=255)
                                         cAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGAc                                                                                                    <  1:232267/147‑1 (MQ=255)
                                         cAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGa                                                                                                  >  7:151487/1‑149 (MQ=255)
                                         cAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCAGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGa                                                                                                  <  6:121246/149‑1 (MQ=255)
                                           aGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATa                                                                                                <  2:284318/149‑1 (MQ=255)
                                                tCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGc                                                                                            <  3:285331/148‑1 (MQ=255)
                                                  gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCggg                                                                                         <  2:82595/149‑1 (MQ=255)
                                                      aGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTcag                                                                                     <  1:46208/149‑1 (MQ=255)
                                                         ttGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTcagcag                                                                                  >  7:273385/1‑149 (MQ=255)
                                                             ccTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAt                                                                                                                              >  4:197740/1‑101 (MQ=255)
                                                             ccTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAt                                                                                                                              <  3:197740/101‑1 (MQ=255)
                                                              cTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTACTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAgg                                                                             >  1:103492/1‑149 (MQ=255)
                                                                tgttGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGa                                                                                                                          >  3:115586/4‑102 (MQ=255)
                                                                tgttGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGa                                                                                                                          <  4:115586/99‑1 (MQ=255)
                                                                tgttGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGca                                                                           >  3:144317/4‑149 (MQ=255)
                                                                        gcCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGGCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGaa                                                                   <  4:214739/149‑1 (MQ=255)
                                                                          cATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAgg                                                                 <  2:138388/149‑1 (MQ=255)
                                                                                    aaCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGAGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGt                                                       >  5:248052/1‑149 (MQ=255)
                                                                                     aCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTg                                                      >  7:66485/1‑149 (MQ=255)
                                                                                        gCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTggcg                                                   >  4:57590/1‑149 (MQ=255)
                                                                                                 cTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGt                                                                <  4:182474/127‑1 (MQ=255)
                                                                                                 cTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGt                                                                >  3:182474/1‑127 (MQ=255)
                                                                                                    cgtcgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACcagca                                         <  8:241575/147‑1 (MQ=255)
                                                                                                    cgtcgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACcagca                                         <  3:57590/147‑1 (MQ=255)
                                                                                                       cgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGcc                                     >  3:139024/1‑148 (MQ=255)
                                                                                                                   cgcgGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGc                         <  4:144317/148‑1 (MQ=255)
                                                                                                                            cGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCATGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCg               >  6:145949/1‑149 (MQ=255)
                                                                                                                                ggTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCGCAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCa            <  1:45058/148‑1 (MQ=255)
                                                                                                                                 gTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATa          <  8:151487/149‑1 (MQ=255)
                                                                                                                                  tCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGcac                                                                          >  2:187721/1‑84 (MQ=255)
                                                                                                                                  tCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGcac                                                                          <  1:187721/84‑1 (MQ=255)
                                                                                                                                  tCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATAc         <  6:248052/149‑1 (MQ=255)
                                                                                                                                   caagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACt        >  6:73522/4‑149 (MQ=255)
                                                                                                                                    aagTCGCCGACGGCGACTTGCTCGGCGCATTCGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAg                                                                  <  3:124159/88‑1 (MQ=255)
                                                                                                                                    aagTCGCCGACGGCGACTTGCTCGGCGCATTCGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAg                                                                  >  4:124159/3‑90 (MQ=255)
                                                                                                                                      gTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGt     >  8:206304/1‑149 (MQ=255)
                                                                                                                                      gTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGt     <  5:73522/149‑1 (MQ=255)
                                                                                                                                         gCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTgg  >  3:276258/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTGG  >  NC_002947/381715‑382000

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: ‑