Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,284 C→T 16.7% intergenic (+30/+61) ibpA → / ← PP_1983 small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,246,2840CT16.7% 74.7 / 9.4 36intergenic (+30/+61)ibpA/PP_1983small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (14/16);  new base T (3/3);  total (17/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 = 7.57e-01

AAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGG  >  NC_002947/2246136‑2246432
                                                                                                                                                    |                                                                                                                                                      
aaGTGAGGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGAAACGAAAAAAAGGGCACCTCCa                                                                                                                                                       >  1:40878/1‑150 (MQ=255)
   tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                    <  1:227659/150‑1 (MQ=255)
           gccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCACATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tt                                                                                                                                          <  2:98732/150‑1 (MQ=255)
                  tGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTg                                                                                                                                   <  1:300448/150‑1 (MQ=255)
                      cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTtc                                                                                                                               >  1:791966/1‑150 (MQ=255)
                           gCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACg                                                                                                                          <  2:215424/150‑1 (MQ=255)
                                    gCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCcgc                                                                                                                 <  1:307512/150‑1 (MQ=255)
                                        cGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                             <  2:40878/150‑1 (MQ=255)
                                                        gTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                             >  2:148091/1‑150 (MQ=255)
                                                            cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                         >  1:483201/1‑150 (MQ=255)
                                                            cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                         >  1:119944/1‑150 (MQ=255)
                                                             tGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                        <  1:794042/150‑1 (MQ=255)
                                                                  aagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTaa                                                                                   <  1:637830/150‑1 (MQ=255)
                                                                        aggcgaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCAC‑‑TTTTTGCGGCTTCTAACGCTTATCCGCGTAGGTGCGCGGCCGTCATGCGATGAGTAGACTGcc                                                                             <  2:626089/150‑1 (MQ=255)
                                                                            caagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                         >  2:791896/1‑150 (MQ=255)
                                                                             aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                        >  1:411099/1‑150 (MQ=255)
                                                                                cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                                     <  2:411099/150‑1 (MQ=255)
                                                                                  cATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCa                                                                   <  1:343696/150‑1 (MQ=255)
                                                                                        gATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACa                                                             >  1:28383/1‑150 (MQ=255)
                                                                                          tCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATc                                                           <  2:483201/150‑1 (MQ=255)
                                                                                                gCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGa                                                     >  1:228933/1‑150 (MQ=255)
                                                                                                gCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGa                                                     <  2:228933/150‑1 (MQ=255)
                                                                                                      tctcAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTcc                                               <  1:167642/147‑1 (MQ=255)
                                                                                                       aCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcg                                              <  2:28383/150‑1 (MQ=255)
                                                                                                          aGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCGTCTGGCGCTTCCAGCAACATCTGCCGTCACTcccaca                                           >  2:167642/1‑146 (MQ=255)
                                                                                                                  aCTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAAt                                   <  2:791966/150‑1 (MQ=255)
                                                                                                                      aTCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCa                                                                   <  2:32219/114‑1 (MQ=255)
                                                                                                                      aTCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCa                                                                   >  1:32219/1‑114 (MQ=255)
                                                                                                                                  cacgaaaaaGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                    <  1:588419/145‑1 (MQ=255)
                                                                                                                                  cacgaaaaaGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                    >  2:588419/5‑149 (MQ=255)
                                                                                                                                     aaaaaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTAACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCCGATAAc                >  1:759368/1‑150 (MQ=255)
                                                                                                                                         aaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTTCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGTCGGCGCCTCCAGCAACATTTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg            >  2:267101/1‑150 (MQ=255)
                                                                                                                                          aGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                      <  2:388847/139‑1 (MQ=255)
                                                                                                                                          aGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                      >  1:388847/1‑139 (MQ=255)
                                                                                                                                            ggCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg           >  2:96335/1‑150 (MQ=255)
                                                                                                                                             gCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTAc        >  1:617816/1‑150 (MQ=255)
                                                                                                                                                   ccAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGGGCCGGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATgg  <  1:247383/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGG  >  NC_002947/2246136‑2246432

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