Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,285 A→G 21.5% intergenic (+31/+60) 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,2850AG21.5% 87.1 / 19.9 47intergenic (+31/+60)ibpA/PP_1983small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base G (5/5);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

AAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246141‑2246431
                                                                                                                                                |                                                                                                                                                    
aaGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTTAACTGATCGAGCTGCAACGAAAAAAGGGCACCTCCAGAGGTg                                                                                                                                                >  1:570325/1‑150 (MQ=255)
   gccgccGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                              >  2:1227505/1‑150 (MQ=255)
        cgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttt                                                                                                                                       <  2:264757/150‑1 (MQ=255)
        cgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑tttt                                                                                                                                       >  2:1274954/1‑150 (MQ=255)
          gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                     <  2:526644/150‑1 (MQ=255)
          gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                     <  1:560851/150‑1 (MQ=255)
          gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                     <  1:297362/150‑1 (MQ=255)
            cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTTAACTGATCGAGCTGCAACGAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTg                                                                                                                                  <  2:570325/150‑1 (MQ=255)
            cTGGCCAATGGCCTGCTCAGCACCGACCGTGTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCAGTCCGATCAATGGCGAAAAACCAGAGCTGAACTGATCGAGCTGCAACGAACCGAAAGTTACCTCCAGCGGTGCCC‑‑TTTTTCGt                                                                                                                                   <  1:37862/150‑1 (MQ=255)
             tGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGTCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACG‑‑AAAAAGGGCACCTCTGGAGGTGCCC‑‑tttttttcgtt                                                                                                                                <  2:1346242/150‑7 (MQ=38)
             tGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACG‑‑AAAAAGGGCACCTCTGGAGGTGCCC‑‑tttttttcgtt                                                                                                                                >  1:1114771/1‑144 (MQ=38)
             tGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACG‑‑AAAAAGGGCACCTCTGGAGGTGCCC‑‑tttttttcgtt                                                                                                                                <  2:1114771/150‑7 (MQ=38)
                       ccTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCCTGCTTATAACga                                                                                                                        >  2:944569/1‑149 (MQ=255)
                         tGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAg                                                                                                                      >  1:1312778/1‑150 (MQ=255)
                              aGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCcg                                                                                                                 >  1:146546/1‑150 (MQ=255)
                                 aTCGACCTGCTGCGCATCGTGCCTGAAAAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGt                                                                                                              >  1:1076646/1‑150 (MQ=255)
                                  tCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                             >  2:562883/1‑150 (MQ=255)
                                    gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAgg                                                                                                           >  1:326092/1‑150 (MQ=255)
                                        tgctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTgcg                                                                                                       >  2:729052/1‑150 (MQ=255)
                                             cgcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGcc                                                                                                  <  1:127920/150‑1 (MQ=255)
                                                        tGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                       >  1:20823/1‑150 (MQ=255)
                                                         gaagaaGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                      <  1:986324/150‑1 (MQ=255)
                                                         gaagaaGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                      <  2:281913/150‑1 (MQ=255)
                                                               gccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATc                                                                                <  2:1312778/150‑1 (MQ=255)
                                                                  aaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGc                                                                             <  2:20823/150‑1 (MQ=255)
                                                                      ccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tt                                                                                                                                         <  1:1150357/86‑1 (MQ=255)
                                                                      ccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tt                                                                                                                                         >  2:1150357/1‑86 (MQ=255)
                                                                        aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                       <  1:562883/150‑1 (MQ=255)
                                                                             cATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCa                                                                  <  1:1227505/150‑1 (MQ=255)
                                                                              aTTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAg                                                                 >  2:615008/1‑150 (MQ=255)
                                                                                  cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAAc                                                             <  1:1274954/150‑1 (MQ=255)
                                                                                        aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCATCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATTCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       >  2:675467/1‑150 (MQ=255)
                                                                                        aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       <  1:944569/150‑1 (MQ=255)
                                                                                                        gcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                       <  1:615008/150‑1 (MQ=255)
                                                                                                               tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGGAACGGTAATgg                                >  2:517097/1‑150 (MQ=255)
                                                                                                                      gCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtg                         >  1:1304781/1‑150 (MQ=255)
                                                                                                                        tGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                       >  2:1065360/1‑150 (MQ=255)
                                                                                                                         gCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTAAAACGCCGTCCCGTGAAGGTGCGCCACGTTGATGCGAGCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGc                      >  1:1242825/1‑150 (MQ=255)
                                                                                                                              gAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATa                 >  1:1083674/1‑150 (MQ=255)
                                                                                                                               aaaaaaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATaa                <  1:729052/150‑1 (MQ=255)
                                                                                                                                aaaaaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAc               <  1:261985/150‑1 (MQ=255)
                                                                                                                                  aaaaGGGCACCTCCCGAGGTGCCC‑‑TTTTTCTTGCTTCTCACTCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAccc             <  2:1083674/150‑1 (MQ=255)
                                                                                                                                   aaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATAGCCGGCTATGCAGATAACCCt            >  2:983063/1‑150 (MQ=255)
                                                                                                                                      gggCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTGCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg           >  1:828592/1‑150 (MQ=255)
                                                                                                                                      gggCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGt         >  1:595960/1‑150 (MQ=255)
                                                                                                                                            cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt   <  2:326092/150‑1 (MQ=255)
                                                                                                                                             tCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg  >  2:120371/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246141‑2246431

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