Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,153,115 G→T 100% L192I (CTT→ATT)  walK ← two‑component sensor histidine kinase [YycG]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,153,1150GT92.1% 64.2 / 0.4 26L192I (CTT→ATT) walKtwo‑component sensor histidine kinase [YycG]
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (14/9);  total (16/10)
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.11e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCATGTCCATACTTTTTAACCTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAAGTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATCGCGCCCACAACCTCTTGGTTTTCAGTCATGACGGGCTTCGCGG  >  BS168_reference_genome/4152985‑4153243
                                                                                                                                  |                                                                                                                                
tcatGTCCATACTTTTTAACCTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCCAAACAATTGCCAAAc                                                                                                                          >  8:61170/1‑139 (MQ=255)
           cTTTTTAACCTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGa                                                                                                                <  3:97123/138‑1 (MQ=255)
              tttAACCTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGcc                                                                                                            >  4:124745/1‑139 (MQ=255)
              tttAACCTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAAATGCCAAACCGGTGCCGGACGc                                                                                                             >  6:280/1‑138 (MQ=255)
                   ccTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGGCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTAt                                                                                                       <  2:184973/139‑1 (MQ=255)
                      tCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCgt                                                                                                    >  1:72149/1‑139 (MQ=255)
                           gagaAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGa                                                                                               <  2:72149/139‑1 (MQ=255)
                                                          gCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTTCCGGACGCCAGTATAGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCAt                                                                >  6:202933/1‑139 (MQ=255)
                                                          gCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGGCTTGATTTGATTAAAGACATCTTCCa                                                                 >  4:57049/1‑138 (MQ=255)
                                                                gCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGCTCGTCTTCATTTGATTAAATCCATCTTCCATACTcgcg                                                         >  8:215052/1‑139 (MQ=255)
                                                                  atatCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGa                                                        <  3:77743/139‑1 (MQ=255)
                                                                     tCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGa                                                                                               <  8:77680/97‑1 (MQ=255)
                                                                     tCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGa                                                                                               >  7:77680/1‑97 (MQ=255)
                                                                                gtgggtAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATcgcg                                          >  3:47762/1‑139 (MQ=255)
                                                                                gtgggtAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATcgcg                                          >  2:49079/1‑139 (MQ=255)
                                                                                   ggtAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAAGTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGAt                                                                                              <  2:154631/84‑1 (MQ=255)
                                                                                   ggtAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAAGTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGAt                                                                                              >  1:154631/1‑84 (MQ=255)
                                                                                             cTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGTCGCCAGTTTGGTGTTGAGCGTGGTGGTTTGATTTAAGCCCTCTTCCATACGCGCGCCACCGTAGAGCGCGCCCCCAcccgcgg                             >  2:125503/1‑135 (MQ=255)
                                                                                             cTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGCCAACGTAGCTCGCGCCCACAACctcgt                             >  7:86751/1‑137 (MQ=255)
                                                                                                 ccAGAAAAATACCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATCGCGCCCACAACCTCTTGGtt                         <  4:47762/139‑1 (MQ=255)
                                                                                                           aCCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCg                                                                                                                  <  7:11380/40‑1 (MQ=255)
                                                                                                           aCCGAGAAGAGCCGTCAAAACAATTGCCAAACCTGTGCCg                                                                                                                  >  8:11380/1‑40 (MQ=255)
                                                                                                               agaagaGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACg                                                   >  4:225409/1‑99 (MQ=255)
                                                                                                               agaagaGCCGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACg                                                   <  3:225409/99‑1 (MQ=255)
                                                                                                                       cGTCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATCGCGCCCACAACCTCTTGGTTTTCAGTCATGACGGGCTTcgcg   <  7:61170/139‑1 (MQ=255)
                                                                                                                         tCAAAACAATTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATCGCGCCCACAACCTCTTGGTTTTCAGTCATGACGGGCTTCGCgg  >  7:66331/1‑138 (MQ=255)
                                                                                                                                  |                                                                                                                                
TCATGTCCATACTTTTTAACCTTCCTCGAGAAATTCCCTTTCGCAAGCTCCATCGCCTGCTTTCGCATATCTGAAAGCGGGTGGGTAATGGTTCTTGCCAGAAAAATACCGAGAAGAGCCGTCAAAACAAGTGCCAAACCTGTGCCGGACGCCAGTATCGTGTTGATCGTCTTCATTTGATTAAAGACATCTTCCATACTCGCGACAACGTAGATCGCGCCCACAACCTCTTGGTTTTCAGTCATGACGGGCTTCGCGG  >  BS168_reference_genome/4152985‑4153243

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: