Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,222,249 C→A 9.8% G321C (GGC→TGC)  PP_1067 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,222,2490CA 9.8% 107.7 / 4.5 41G321C (GGC→TGC) PP_1067sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (21/16);  new base A (2/2);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

TTCGGCAATCCGCCGCTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCGGGAAGGTCAG  >  NC_002947/1222106‑1222396
                                                                                                                                               |                                                                                                                                                   
ttCGGCAATCCGCCGCTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCa                                                                                                                                                >  3:295020/1‑149 (MQ=255)
ttCGGCAATCCGCCGCTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCa                                                                                                                                                >  4:458/1‑149 (MQ=255)
        tCCGCCGCTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGAt                                                                                                                                        <  2:228311/149‑1 (MQ=255)
              gcTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGAtt                                                                                                                                       <  1:300990/144‑1 (MQ=255)
              gcTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGAtt                                                                                                                                       >  2:300990/1‑144 (MQ=255)
                  cAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCa                                                                                                                              <  6:157976/149‑1 (MQ=255)
                        gctgTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCACTGCACCATGGATTCGCGACGCAGGTCCt                                                                                                                        >  2:303153/1‑149 (MQ=255)
                        gctgTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCACTGCACCATGGATTCGCGACGCAGGTCCt                                                                                                                        >  2:307941/1‑149 (MQ=255)
                          tgTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCTCGACGCAGGTCCTg                                                                                                                       >  8:147427/1‑148 (MQ=255)
                               gCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTgg                                                                                                                 <  7:153366/149‑1 (MQ=255)
                                    ttCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGt                                                                                                              <  1:302818/147‑1 (MQ=255)
                                             ttCTCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTg                                                                                                    >  2:133769/1‑148 (MQ=255)
                                               cgcgGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTggg                                                                                                                <  4:31981/134‑1 (MQ=255)
                                               cgcgGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTggg                                                                                                                >  3:31981/1‑134 (MQ=255)
                                                  gggCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGCATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGc                                                                                              >  1:131257/1‑149 (MQ=255)
                                                      cgccAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAgg                                                                                          >  4:194310/1‑149 (MQ=255)
                                                       gccAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAggg                                                                                         >  2:257329/1‑149 (MQ=255)
                                                          aGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTTCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGgt                                                                                      <  3:233827/149‑1 (MQ=255)
                                                          aGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGgt                                                                                      >  4:37557/1‑149 (MQ=255)
                                                                     gATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGgc                                                                           >  5:71426/1‑149 (MQ=255)
                                                                       tCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGgcgc                                                                         <  4:295020/149‑1 (MQ=255)
                                                                           tttGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTgg                                                                     >  4:209493/1‑149 (MQ=255)
                                                                             tGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCACTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGt                                                                    <  1:307941/148‑1 (MQ=255)
                                                                             tGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCACTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGt                                                                    <  1:303153/148‑1 (MQ=255)
                                                                                  cgcTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTgcgc                                                              >  7:190810/1‑149 (MQ=255)
                                                                                   gcTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTgcgc                                                              >  1:73818/1‑148 (MQ=255)
                                                                                       gTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGa                                                          <  2:73818/148‑1 (MQ=255)
                                                                                          ccAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAgg                                                      >  5:2091/1‑149 (MQ=255)
                                                                                                           tCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTg                                     <  6:71426/149‑1 (MQ=255)
                                                                                                                    aGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCg                            <  3:37557/149‑1 (MQ=255)
                                                                                                                        aaaGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGttt                        >  4:54419/1‑149 (MQ=255)
                                                                                                                         aaGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTc                       <  1:257329/149‑1 (MQ=255)
                                                                                                                          aGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGTGTGTCGACCAGGCGCCTGGTTTGCTCCAGATAGGCCACTTCGTGCGTTTTGCCATGAGCGGTTTCg                      <  1:133769/149‑1 (MQ=255)
                                                                                                                                 gCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCg                      <  2:74755/142‑1 (MQ=255)
                                                                                                                                 gCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCg                      >  1:74755/1‑142 (MQ=255)
                                                                                                                                 gCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTgg               >  8:206859/1‑149 (MQ=255)
                                                                                                                                    caccaGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCg            <  3:209493/149‑1 (MQ=255)
                                                                                                                                       caGAATGCCCTGCACAATGGATTCGCGACGCAGGTCCTGCAGTGGGGACAGCAGGGTGAAATGCCAGGGGGAGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTAGGTGTTGTCCggg          >  7:293034/1‑148 (MQ=255)
                                                                                                                                        aGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCggg          >  5:70429/1‑147 (MQ=255)
                                                                                                                                         gAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCGGGAAg       >  5:50802/1‑149 (MQ=255)
                                                                                                                                           aTGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCGGGAAGGt     <  3:194310/149‑1 (MQ=255)
                                                                                                                                              cccTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCGGGAAGGTCAg  <  7:147427/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTCGGCAATCCGCCGCTCCAGCTGGCTGTTGGCCTCTTCAAGGGCTTCGCGGGCCGCCAGGCGGGTGGCGATCACTTTGCGCCGCTCGTTCCAGGCAATGCCGAGGATCGCCAGCAAGGCAAAGGCCACGCCCACCAGAATGCCCTGCACCATGGATTCGCGACGCAGGTCCTGCAGTGGGGTCAGCAGGGTGAAATGCCAGGGGGTGTCGACCAGGCGCCTGGTTTGCGCCAGATAGGCCACTTCGTGCGGTTTGCCATGAGCGGTTTCGGTGTTGGCCGGGAAGGTCAG  >  NC_002947/1222106‑1222396

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