Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,470,656 Δ1 bp 12.1% coding (178/1182 nt) PP_2165 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,470,6560C.12.1% 87.1 / 8.9 33coding (178/1182 nt)PP_2165two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (16/12);  new base . (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

ACGACGACGTGGTCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCA  >  NC_002947/2470516‑2470794
                                                                                                                                             |                                                                                                                                          
acgacgacgTGGTCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCGTTGCGCAGATgg                                                                                                                                   <  2:463722/150‑1 (MQ=255)
      acgTGGTCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGcc                                                                                                                             >  2:211860/1‑150 (MQ=255)
        gTGGTCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTc                                                                                                                           >  2:501039/1‑150 (MQ=255)
           gTCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGaa                                                                                                                        <  1:501039/150‑1 (MQ=255)
            tCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAAc                                                                                                                       <  2:196468/150‑1 (MQ=255)
            tCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAAc                                                                                                                       <  2:577276/150‑1 (MQ=255)
                          gacgccTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGTCTCGAACTTATCCgtcaggtc                                                                                                         >  2:354092/3‑150 (MQ=255)
                                    ttGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTg                                                                                               >  1:439719/1‑150 (MQ=255)
                                            aGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGTGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAg                                                                                       >  1:710593/1‑150 (MQ=255)
                                                 tttCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGGCATGCGCAGAtc                                                                                                                                    <  1:123193/100‑2 (MQ=255)
                                                 tttCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGGCATGCGCAGAtc                                                                                                                                    >  2:123193/1‑99 (MQ=255)
                                                 tttCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGcc                                                                                  >  2:272195/1‑150 (MQ=255)
                                                     aGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTg                                                                              >  1:263107/1‑150 (MQ=255)
                                                        gTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATc                                                                           <  1:272195/150‑1 (MQ=255)
                                                           ctccAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGTGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCgtg                                                                        <  2:710593/150‑1 (MQ=255)
                                                              cAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTg                                                                              >  2:141867/1‑141 (MQ=255)
                                                              cAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTg                                                                              <  1:141867/141‑1 (MQ=255)
                                                                 gCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTcc                                                                  >  2:773253/1‑150 (MQ=255)
                                                                             cagcagGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTc                                                      <  1:211860/150‑1 (MQ=255)
                                                                              agcagGGACGTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCa                                                     <  1:511816/150‑1 (MQ=255)
                                                                                cagGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATg                                                   >  2:468095/1‑150 (MQ=255)
                                                                                  ggggCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAg                                                 <  2:439719/150‑1 (MQ=255)
                                                                                    ggCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTg                                               <  2:263107/150‑1 (MQ=255)
                                                                                                     gaaCACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTg                              >  1:753832/1‑150 (MQ=255)
                                                                                                      aaCACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTgc                             <  2:660676/150‑1 (MQ=255)
                                                                                                           cAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCTTg                        >  1:321650/1‑150 (MQ=255)
                                                                                                           cAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTg                        >  1:98612/1‑150 (MQ=255)
                                                                                                                             tgcgatctgcgCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCgtggt                                                                      <  1:579759/86‑1 (MQ=255)
                                                                                                                             tgcgatctgcgCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCgtggt                                                                      >  2:579759/1‑86 (MQ=255)
                                                                                                                             tgcgatctgcgCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCGGACGACTACCtc      >  2:178357/1‑150 (MQ=255)
                                                                                                                               ccatctgcgGCATGCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGt                                        <  2:336688/112‑1 (MQ=255)
                                                                                                                               ccatctgcgGCATGCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGt                                        >  1:336688/3‑114 (MQ=255)
                                                                                                                                 atctgcgCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTTCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGTCTGGGCGCCGCCGCCTCCCtcatcc  >  2:345914/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                          
ACGACGACGTGGTCCGTGCCAGCCTCGCCGCCTATCTTGAAGACAGTGGTTTCAGCGTCCTCCAGGCCGGTAATGGCCAGCAGGGGCTTCAGGTCTTCGAAGAACACCAGCCCGACCTCGTGATCTGCGATCTGCGCATGCCGCAGATGGGCGGCCTCGAACTGATCCGTCAGGTCAGCGAGCGTGCGCCGCAGTTGCCGGTGATCGTGGTGTCCGGTGCCGGCGTCATGAGTGATGCGGTGGAAGCGTTGCGCCTGGGCGCCGCCGACTACCTCATCA  >  NC_002947/2470516‑2470794

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: