Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,474,890 C→T 100% L45L (TTG→TTA kre ← regulator of transcription factor ComK function via modulation of mRNA stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,474,8900CT100.0% 90.4 / NA 28L45L (TTG→TTAkreregulator of transcription factor ComK function via modulation of mRNA stability
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/16);  total (12/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGCAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGT  >  BS168_reference_genome/1474769‑1475006
                                                                                                                         |                                                                                                                    
gCATGTGAAAATAGTAGTCTCCGTTTGAGATAAGCACGTCCGATTGCTGTTGGCTGAATCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcag                                                                                                     <  5:265446/139‑1 (MQ=255)
   tgtgAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACTAGCTTTAGTAAAGCCCGTTTTTTctg                                                                                                     >  4:176749/1‑134 (MQ=255)
   tgtgAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTATTGTAAAGCCAGATTTTTcag                                                                                                     >  8:160097/1‑136 (MQ=255)
   tgtgAAAATAGTAGTCTCCGATTGAGATAAGAACGTCAGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAATCGCGttatttacg                                                                                                     >  5:131511/1‑133 (MQ=255)
      gAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGata                                                                                               <  4:173749/139‑1 (MQ=255)
              tagtCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTAtt                                                                                       >  6:69280/1‑139 (MQ=255)
                   tcCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAggg                                                                                  >  5:147080/1‑139 (MQ=255)
                           agaTAAGCACGTGCGATTGATGTTGTCTGAACCTTGGATTTTTTGAAAAATTAAGCCCTATTTTCTTTCATTTACCTTCTTTGACAATCTTTTGTATAGCACGTTTTTTCAGCAGATATAGGTATTTAGGGTTggggg                                                                           <  6:131511/138‑1 (MQ=255)
                                     gTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTAACTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGAAGTTAGGGTTGGGGGCGGTTTTTGCa                                                                >  1:76570/1‑139 (MQ=255)
                                      tCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcag                                                                                                     <  8:25544/101‑1 (MQ=255)
                                      tCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcag                                                                                                     >  7:25544/1‑101 (MQ=255)
                                               tGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGttt                                                                     >  8:144484/1‑124 (MQ=255)
                                               tGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGttt                                                                     <  7:144484/124‑1 (MQ=255)
                                               tGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTca                                                      <  4:23943/139‑1 (MQ=255)
                                                gTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                     <  5:69280/139‑1 (MQ=255)
                                                 ttcgCTGAACCTTGGATTATTAGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTGca                                                      <  4:38815/134‑1 (MQ=255)
                                                 ttGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                     <  3:176749/138‑1 (MQ=255)
                                                 ttGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                     <  2:252704/138‑1 (MQ=255)
                                                              ggATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATg                                                              >  6:119941/1‑116 (MQ=255)
                                                              ggATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATg                                                              <  5:119941/116‑1 (MQ=255)
                                                                                       tgttctttcCTTTACCTTCTTTGACAAGCTTTTGTAAAGCGCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTATACTGCTTTTGAAAGATTTTCTACGGTTGGCTg              <  2:76570/137‑1 (MQ=255)
                                                                                            tttcCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGa          >  7:169208/1‑138 (MQ=255)
                                                                                                 tttACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGt                                                                       >  2:110190/1‑72 (MQ=255)
                                                                                                 tttACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGt                                                                       <  1:110190/72‑1 (MQ=255)
                                                                                                 tttACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTc    >  6:30473/1‑139 (MQ=255)
                                                                                                    aCCTTCTTTTACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGGTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGt  <  7:160097/138‑1 (MQ=255)
                                                                                                    aCCTTCTTTGACTAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTctt  <  6:147080/138‑3 (MQ=255)
                                                                                                    aCCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGt  <  8:169208/138‑1 (MQ=255)
                                                                                                                         |                                                                                                                    
GCATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGCAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGT  >  BS168_reference_genome/1474769‑1475006

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: