Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 975,291 A→T 8.4% intergenic (+137/+32) ycfJ → / ← suhB flagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947975,2910AT 8.4% 138.9 / 4.7 49intergenic (+137/+32)ycfJ/suhBflagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase
Reads supporting (aligned to +/- strand):  ref base A (19/25);  new base T (2/2);  total (22/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 = 3.35e-01

AGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAAG  >  NC_002947/975150‑975440
                                                                                                                                             |                                                                                                                                                      
aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCt                                                                                                                                               >  2:535442/1‑150 (MQ=255)
 gTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCtt                                                                                                                                              >  2:75287/1‑150 (MQ=255)
     tttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTtttg                                                                                                                                          >  1:775045/1‑150 (MQ=255)
      ttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAGCACCTCTTGAGGTGCttttttgt                                                                                                                                        <  2:81131/150‑4 (MQ=38)
                tgtgCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTtttgtt                                                                                                                                        >  2:443970/1‑141 (MQ=255)
                tgtgCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTtttgtt                                                                                                                                        <  1:443970/141‑1 (MQ=255)
                     cGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgc                                                                                                                          <  1:78635/150‑1 (MQ=255)
                       gCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgcgc                                                                                                                        <  2:101056/150‑1 (MQ=255)
                            ttCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCCAAGTAGGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAgg                                                                                                                   >  1:8392/1‑149 (MQ=39)
                               gcgGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCt                                                                                                                <  2:122112/150‑1 (MQ=255)
                                   gTAAACCCGATCCTACAGGGGCAATGCGATCCGTTTTCGAACGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCAt                                                                                                            <  2:496164/150‑1 (MQ=255)
                                     aaaCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGc                                                                                                          <  2:433965/150‑1 (MQ=255)
                                       aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCa                                                                                                             <  1:718972/145‑1 (MQ=255)
                                       aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCa                                                                                                             >  2:718972/1‑145 (MQ=255)
                                        cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCttt                                                                                                       >  1:128001/1‑150 (MQ=255)
                                          cGATCCTATAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt                                                                                                     <  2:168030/150‑1 (MQ=255)
                                                         aaTGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTgg                                                                                      <  2:775045/150‑1 (MQ=255)
                                                                     gTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagtcagt                                                                          <  1:428837/150‑1 (MQ=255)
                                                                         gAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCg                                                                      >  2:271577/1‑150 (MQ=255)
                                                                           gCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgc                                                                                                                          <  1:788112/96‑1 (MQ=255)
                                                                           gCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgc                                                                                                                          >  2:788112/1‑96 (MQ=255)
                                                                               gTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGa                                                                >  1:852923/1‑150 (MQ=255)
                                                                                     ccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCttttt                                                                                                                                           <  1:790857/69‑1 (MQ=255)
                                                                                     ccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCttttt                                                                                                                                           >  2:790857/1‑69 (MQ=255)
                                                                                       gcgAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTg                                                        <  2:128001/150‑1 (MQ=255)
                                                                                                tCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc                                                <  2:786201/149‑1 (MQ=255)
                                                                                                tCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc                                                >  1:786201/1‑149 (MQ=255)
                                                                                                       aGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt                                                                                                     <  2:681519/89‑1 (MQ=255)
                                                                                                       aGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt                                                                                                     >  1:681519/1‑89 (MQ=255)
                                                                                                        ggCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt                                       <  2:653892/150‑1 (MQ=255)
                                                                                                        ggCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt                                       <  1:535442/150‑1 (MQ=255)
                                                                                                                  aaCCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTg                             >  1:105993/1‑150 (MQ=255)
                                                                                                                      accaGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGa                         <  2:804696/150‑1 (MQ=255)
                                                                                                                      accaGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGa                         <  2:852923/150‑1 (MQ=255)
                                                                                                                       ccaGGCACAAAAAAGAGCTTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAg                        <  2:326471/150‑1 (MQ=255)
                                                                                                                          ggCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGaa                     <  1:75287/150‑1 (MQ=255)
                                                                                                                           gCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTgg                                                                                      >  2:571403/1‑84 (MQ=255)
                                                                                                                           gCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTgg                                                                                      <  1:571403/84‑1 (MQ=255)
                                                                                                                           gCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAg                    >  2:371257/1‑150 (MQ=255)
                                                                                                                             acaAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTc                  >  2:433956/1‑150 (MQ=255)
                                                                                                                             aacaaaaaGCACCTCTTGAGGTGCTTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGt                   >  1:6697/4‑150 (MQ=255)
                                                                                                                               caaaaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                    <  1:111612/65‑1 (MQ=37)
                                                                                                                               caaaaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                    >  2:111612/2‑66 (MQ=37)
                                                                                                                               aaaaaaGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt                                       >  2:176542/1‑127 (MQ=255)
                                                                                                                               aaaaaaGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt                                       <  1:176542/127‑1 (MQ=255)
                                                                                                                                aaaaaGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTg               <  1:433956/150‑1 (MQ=255)
                                                                                                                                 aaaaGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTgg              <  2:105993/150‑1 (MQ=255)
                                                                                                                                         cTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGtt      >  1:453373/1‑150 (MQ=255)
                                                                                                                                           cAAGAGATGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTAATGATTTCTGGAACGAGAGGAGAAATGGCTGTCAGTACCGGCTTGGAGAAGGTGAGTTTGTTTAGATGGATGTGGCCTTTGTCGAGGAAGTCGTGGCCGCCGTTGa    >  2:735753/1‑150 (MQ=255)
                                                                                                                                            aaGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGaa   <  1:501307/150‑1 (MQ=255)
                                                                                                                                             agagGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAAg  <  2:189666/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
AGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAAG  >  NC_002947/975150‑975440

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