Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 427,920 C→G 16.5% S225C (TCT→TGT)  yclK → two‑component sensor histidine kinase [YclJ]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome427,9200CG16.5% 128.6 / 4.4 51S225C (TCT→TGT) yclKtwo‑component sensor histidine kinase [YclJ]
Reads supporting (aligned to +/- strand):  ref base C (17/24);  new base G (7/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 6.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTGAGCTGGCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACATGAATTAAAAACACCGCTGACGACGATCAAC  >  BS168_reference_genome/427784‑428047
                                                                                                                                        |                                                                                                                               
cTTGTGAGCTGGCTTCTGTCCTAATTCCATGTCACGCGGATTCAAAAGCTGAGCGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                <  5:507/138‑1 (MQ=255)
cTTGTGAGCTGGCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                <  5:318165/138‑1 (MQ=255)
cTTGTGAGCTGGCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                <  3:329119/138‑1 (MQ=255)
    tgAGCTGGCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACtt                                                                                                                           <  2:236194/139‑1 (MQ=255)
      aGCTGGCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACttt                                                                                                                          >  3:427722/1‑138 (MQ=255)
           gCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAatat                                                                                                                     >  8:295105/1‑138 (MQ=255)
           gCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACATTAatat                                                                                                                     >  8:218229/1‑138 (MQ=255)
                tGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTGTGCCTTTAATATCATGGc                                                                                                               >  8:265102/1‑139 (MQ=255)
                  tCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCg                                                                                                              >  3:172186/1‑138 (MQ=255)
                   ccAAATTCCATGTCAAGCGGATTCAAACGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGaa                                                                                                            <  1:339402/139‑1 (MQ=255)
                   ccAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGaa                                                                                                            <  2:257791/139‑1 (MQ=255)
                    cAAATTCCATGTCAAGCGGATTCAAAAGATGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGaa                                                                                                            <  5:357027/138‑1 (MQ=255)
                     aaaTTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTGGGGCTTTAATATAATGGCGaaa                                                                                                           >  2:339402/1‑138 (MQ=255)
                       aTTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTGGGGCGTTTATATAATTGGGaaaaa                                                                                                         >  1:326716/1‑138 (MQ=255)
                       aTTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATTGCGaaaaaa                                                                                                        >  5:155032/1‑139 (MQ=255)
                           taTGTCAAGCGTATTAAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTaa                                                                                                    <  1:274953/138‑1 (MQ=255)
                            aTGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                <  4:109146/110‑1 (MQ=255)
                            aTGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                >  3:109146/1‑110 (MQ=255)
                            aTGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTa                                                                                                                         <  5:335835/117‑1 (MQ=255)
                            aTGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTa                                                                                                                         >  6:335835/1‑117 (MQ=255)
                             tGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAg                                                                                                  <  6:187670/139‑1 (MQ=255)
                                 aaGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTGGGAAATTAAAATATTTGCGAAAAAACAAAATCAGt                                                                                              >  8:156929/1‑139 (MQ=255)
                                  aGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                <  7:305051/104‑1 (MQ=255)
                                  aGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                >  8:305051/1‑104 (MQ=255)
                                       aTTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTGGGACTTTAATATAATGGCGAAAAAAATAAATCAGTCAAggc                                                                                        >  2:274953/1‑138 (MQ=255)
                                          aaacaGCTTAGAGAAGCGACATATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACTGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATg                                                                                     <  2:326716/135‑1 (MQ=255)
                                                 tGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                >  4:267395/1‑89 (MQ=255)
                                                 tGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCt                                                                                                                                <  3:267395/89‑1 (MQ=255)
                                                   agagaAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTGGGACTTTAATATCATGGCGATAAAAATAAAGCCGTCAAGGGGTGGAATCGAgc                                                                            >  7:94185/1‑139 (MQ=255)
                                                   agagaAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGa                                                                              <  4:427722/137‑1 (MQ=255)
                                                   agagaAGCGACAGATAAAGTGGCTTCCGGAGATTATGACATCCAATTGGAAAACAGCTACTGGGACGAAATCTGAGTACTGGCGTCTTTACTTAATATAATGGCTAAAAAAATTAAACAGCAAAGGCATTaaaactacc                                                                            >  8:382491/1‑133 (MQ=255)
                                                     agaAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGCCTGTTATTTAATGGCGAAAAAACTAAAGCAGTCAAGGGATGGAATCGAgcgc                                                                          >  5:286761/1‑139 (MQ=255)
                                                      taaGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGcc                                                                         <  1:87004/138‑1 (MQ=255)
                                                               gATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGc                                                                <  8:132097/139‑1 (MQ=255)
                                                                             cGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGtt                                                   <  2:85572/138‑1 (MQ=255)
                                                                              ggCGATTATGACATCCATTTTGAAAACAGCTACGGGTACGAAATCGGCGTACTGGCGGCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTa                                                 <  8:71249/139‑1 (MQ=255)
                                                                                 gtttATGACATCCATTTCGAAAACAGCTACCGGGACGAACTCGGCGTACTGGCGTGTTAGTTTAATTTGATGGCGAAAAATATAAACCAGTCAAGCGATGTAATGGAGCGGCTTGAGATGCGGAGAAGGCAGTTTATCg                                              <  4:258320/137‑1 (MQ=255)
                                                                                   gtATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCt                                            <  5:421924/138‑1 (MQ=255)
                                                                                           tCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATAAAGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACgtgt                                     >  6:161607/1‑138 (MQ=255)
                                                                                            ccATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACATTAAAATCATTGCGAAAAAACTAAAGCACTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTca                                   >  4:106728/1‑139 (MQ=255)
                                                                                               tttGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACAt                                <  7:295105/139‑1 (MQ=255)
                                                                                               tttGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACAt                                <  4:172186/139‑1 (MQ=255)
                                                                                               tttGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACAt                                <  5:382593/139‑1 (MQ=255)
                                                                                               tttGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCGGAGATTAATATAATGGCGAAAAAAATAAAGAAGTCAAGGGATTAAATACAGCGCCAAGATAAGCGGAGAAGGCAGGTTATCCCTGACgtgtgaaa                                 >  6:421924/1‑134 (MQ=255)
                                                                                               tttGGAAAACAGCTAAGGGGACGAAATCGGCGTTCTGGCGTATGACTTTAATATCATTGAGAAAAAACTAAAGCAGTCAAGTGATGAAATCGAGCGCCTTGAGAAGAGGAGAAGGCATTTTATCGCTGACGTGTCACAt                                <  7:382491/139‑1 (MQ=255)
                                                                                                       aCAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATAGAGCGCCGTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACAttaattaa                         >  2:359296/1‑138 (MQ=255)
                                                                                                                   aCGAAATCGGCGTACTGGCGTCTGACTTTAATATCATTGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACATGAATTAAAAACACCGCTgac            >  3:148357/1‑139 (MQ=255)
                                                                                                                     gAAATGGGCGTCCTGGCGTCTGACTTTAATTTCATGGCTAAAAAACTAAAGCAGTCAAGGGATGAAATGGAGCGCCTTGAGAACCGGAGAAGGCAGTTTATCGCTGACGTGTAACATGAATTAAAAACACCGCTgacga          <  5:23175/139‑1 (MQ=255)
                                                                                                                     gAAATCGGCGTACTGGCGTGTGACTTTAATATCTTGGCGAAAAAACTAATGCAGTCAAGGGATGAAATCGTGGGCCTTGAGAAGGGGAGATGGCAGTTTATCGCTGACGTGTCACATGAATTAAAAACACCGCTgacga          <  7:156929/139‑1 (MQ=255)
                                                                                                                             cGTACTGGCGTGTGTCTTTTATATCATGGGGCAAAAAATTAAGGAGGCAAGGGATGTAATCGAGCGCCGGGTGTAGCGGAGAAGGGAGTTTTTCGCTGACGTGGCACTTGAATTAAAAACACCGCGGACGACGATCAAc  >  6:233934/1‑139 (MQ=255)
                                                                                                                             cGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACATGAATTAAAAACACCGCTGACGACGATCAAc  >  8:352431/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                               
CTTGTGAGCTGGCTTCTGTCCAAATTCCATGTCAAGCGGATTCAAAAGCTGAGAGAAGCGACAGATAAAGTGGCTTCCGGCGATTATGACATCCATTTGGAAAACAGCTACGGGGACGAAATCGGCGTACTGGCGTCTGACTTTAATATCATGGCGAAAAAACTAAAGCAGTCAAGGGATGAAATCGAGCGCCTTGAGAAGCGGAGAAGGCAGTTTATCGCTGACGTGTCACATGAATTAAAAACACCGCTGACGACGATCAAC  >  BS168_reference_genome/427784‑428047

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