Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,420,641 Δ1 bp 9.6% intergenic (+81/+26) ctpH → / ← PP_t62 methyl‑accepting chemotaxis protein CtpH/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,420,6410A. 9.6% 157.8 / 11.3 52intergenic (+81/+26)ctpH/PP_t62methyl‑accepting chemotaxis protein CtpH/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base . (2/3);  total (26/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 = 6.64e-01

AGGCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCTT  >  NC_002947/2420493‑2420787
                                                                                                                                                    |                                                                                                                                                   
aGGCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGcac                                                                                                                                                    >  2:946045/1‑150 (MQ=255)
     gggCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGc                                                                                                                                              >  2:246305/1‑150 (MQ=255)
         aGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGc                                                                                                                                          <  2:706749/150‑1 (MQ=255)
                        tGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGAtt                                                                                                                           <  1:1038469/150‑1 (MQ=255)
                          tCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACAt                                                                                                                                                  >  1:1057748/1‑126 (MQ=255)
                          tCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACAt                                                                                                                                                  <  2:1057748/126‑1 (MQ=255)
                              ggCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCg                                                                                                                     >  1:1102467/1‑150 (MQ=255)
                                   ttGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCATACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGgaga                                                                                                                 >  1:1149626/1‑149 (MQ=255)
                                   ttGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGgaga                                                                                                                 <  2:1149626/149‑1 (MQ=255)
                                    tGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGgagaga                                                                                                               <  1:904211/150‑1 (MQ=255)
                                       tcGAAGTCACCTCCCGCCTGAGTGCCTGCATCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAg                                                                                                            <  2:1102467/149‑1 (MQ=255)
                                           aGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGTAtt                                                                                                        <  1:544022/150‑1 (MQ=255)
                                            gTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCttt                                                                                                                                       <  1:105408/118‑1 (MQ=255)
                                            gTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCttt                                                                                                                                       >  2:105408/1‑118 (MQ=255)
                                                 cTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCCGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGGCGGAGAGACAGGATTCGAAcc                                                                                                  >  1:45051/1‑150 (MQ=255)
                                                    ccGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGCtt                                                                                                                                        <  2:228058/109‑1 (MQ=255)
                                                    ccGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGCtt                                                                                                                                        >  1:228058/1‑109 (MQ=255)
                                                       ccTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGc                                                                                            <  1:301819/150‑1 (MQ=255)
                                                           aGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGCTTtgt                                                                                                                                     <  2:503696/105‑1 (MQ=39)
                                                           aGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCCATGTGCGGGCTTtgt                                                                                                                                     >  1:503696/1‑105 (MQ=39)
                                                                cTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACCTATCATTCGCACAGAAAACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTc                                                                                   >  1:147606/1‑150 (MQ=255)
                                                                 tGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTcc                                                                                  <  1:946045/150‑1 (MQ=255)
                                                                   aCTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCg                                                                                >  1:657405/1‑150 (MQ=255)
                                                                       ggTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGt                                                                                                                       <  1:528570/107‑1 (MQ=255)
                                                                       ggTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGt                                                                                                                       >  2:528570/1‑107 (MQ=255)
                                                                              aGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGga                                                                                                                   <  1:872711/104‑1 (MQ=255)
                                                                              aGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGga                                                                                                                   >  2:872711/1‑104 (MQ=255)
                                                                               ggAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCGTCTCGTTCGGGACgagaaggcc                                                                    <  2:853552/150‑6 (MQ=255)
                                                                                   ttGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTAcc                                                                >  1:1137447/1‑150 (MQ=255)
                                                                                     gCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTg                                                              >  1:279140/1‑150 (MQ=255)
                                                                                       aGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCt                                                                                                 >  1:746212/1‑113 (MQ=255)
                                                                                       aGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCt                                                                                                 <  2:746212/113‑1 (MQ=255)
                                                                                       aGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGc                                                            <  2:657405/150‑1 (MQ=255)
                                                                                            ccATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTgcgc                                                       <  1:246305/150‑1 (MQ=255)
                                                                                              aTAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACACAGCACGCACATGGCGGGCTTTGACAAAAGCATTTGGTCGAAGAGACAGGAGGCTAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTa                                                     >  2:739218/1‑150 (MQ=255)
                                                                                                agagGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTaca                                                   <  2:279140/150‑1 (MQ=255)
                                                                                                                    gACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTa                               <  2:147606/150‑1 (MQ=255)
                                                                                                                     aCGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTAc                              >  2:942480/1‑150 (MQ=255)
                                                                                                                     aCGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTAc                              >  2:344758/1‑150 (MQ=255)
                                                                                                                       gTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTAccg                            >  1:763139/1‑148 (MQ=255)
                                                                                                                       gTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTg                            >  1:542108/1‑150 (MQ=255)
                                                                                                                           ttttCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACCGCATc                        <  2:763139/150‑1 (MQ=255)
                                                                                                                                  acacaacacaAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTgg                                                             >  2:505971/1‑106 (MQ=255)
                                                                                                                                  acacaacacaAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTgg                                                             <  1:505971/106‑1 (MQ=255)
                                                                                                                                   cacaacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGttttt                <  2:737941/150‑1 (MQ=255)
                                                                                                                                      aacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTc                                                                                   <  2:1086221/80‑1 (MQ=255)
                                                                                                                                      aacacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTc                                                                                   >  1:1086221/1‑80 (MQ=255)
                                                                                                                                       acacaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTc            <  1:942480/150‑1 (MQ=255)
                                                                                                                                         acaaAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGc          >  2:411240/1‑150 (MQ=255)
                                                                                                                                             aGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGc      >  1:611097/1‑150 (MQ=255)
                                                                                                                                              gCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACaaaa                                    >  2:602556/1‑119 (MQ=255)
                                                                                                                                              gCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACaaaa                                    <  1:602556/119‑1 (MQ=255)
                                                                                                                                                 cGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTACGGTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCtt  >  2:886367/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGCGGGGCAGCGCCTGCATGGTTTGTCAGGGCAGTTGGGCGAAGTCACCTCCCGCCTGAGTGCCTGACTCGGTCTGCAGGAGTTGCAGGTGCCATAGAGGCTGAAGCCGCGCACTGACGTATTTTTCGCACACAACACAAAGCCCGCACATGGCGGGCTTTGTGACAAGGATTTGGTCGGAGAGACAGGATTCGAACCTGCGGCCTTCTCGTCCCGAACGAGACGCGCTACCTGGCTGCGCTACACTCCGATGAACAAAATCCTACTGCATCCCGTTTTTGCTCGCAAGCCCTT  >  NC_002947/2420493‑2420787

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: