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% 204.5 / ‑6.6 59L45L (TTG→TTAkreregulator of transcription factor ComK function via modulation of mRNA stability
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/30);  minor base A (1/0);  total (29/30)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAGGCATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGCAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCATGGTCGTC  >  BS168_reference_genome/1474765‑1475021
                                                                                                                             |                                                                                                                                   
ggAGGCATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGttttt                                                                                                                        >  4:189883/1‑139 (MQ=255)
     cATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcagc                                                                                                                   <  6:309880/139‑1 (MQ=255)
     cATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAACCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcagc                                                                                                                   <  7:236238/139‑1 (MQ=255)
       tgtgAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcag                                                                                                                    >  4:228744/1‑136 (MQ=255)
       tgtgAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCACTTTTTTTcag                                                                                                                    >  3:144553/1‑136 (MQ=255)
            aaaTAGTAGTCTCCGATTGATATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGatat                                                                                                             <  3:32500/138‑1 (MQ=255)
             aaTAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAg                                                                                                           <  1:180720/139‑1 (MQ=255)
                  tagtCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTAtt                                                                                                      >  5:71995/1‑139 (MQ=255)
                   agtCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTAttt                                                                                                     >  4:111882/1‑139 (MQ=255)
                       tcCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcag                                                                                                                    <  6:306070/120‑1 (MQ=255)
                       tcCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTcag                                                                                                                    >  5:306070/1‑120 (MQ=255)
                          gATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTca                                                                                                                     <  7:124728/116‑1 (MQ=255)
                          gATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTca                                                                                                                     >  8:124728/1‑116 (MQ=255)
                          gATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTg                                                                                              >  6:145130/1‑139 (MQ=255)
                               agaTAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTggggg                                                                                          <  3:122986/138‑1 (MQ=255)
                               agaTAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTACCTTTACCTTCTTTGACAAGCTTTTGAAATGCCCGTTTGTTCAGCAGATATAAGTATTTAGGTTTgaggtc                                                                                         >  2:115555/1‑137 (MQ=255)
                                 aTAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCgg                                                                                       >  6:261930/1‑139 (MQ=255)
                                 aTAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCgg                                                                                       >  6:233010/1‑139 (MQ=255)
                                    aGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGttt                                                                                    <  8:284181/139‑1 (MQ=255)
                                     gCACGTCCGATTGCTGTTGGCTGACCCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTGGACAAGCTTTTGTAAAGCCCGTTGTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGtttt                                                                                   >  2:77503/1‑139 (MQ=255)
                                       aCGTCCGATTGCTGTTGGCTGAACCTTGTATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTg                                                                                 <  4:86471/139‑1 (MQ=255)
                                         gTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCa                                                                               >  8:250755/1‑139 (MQ=255)
                                         gTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCa                                                                               <  3:111882/139‑1 (MQ=255)
                                         gTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCa                                                                               >  8:131730/1‑139 (MQ=255)
                                          tCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCAt                                                                              >  7:199281/1‑139 (MQ=255)
                                            cGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATgg                                                                            <  3:228744/139‑1 (MQ=255)
                                              aTTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCTGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCg                                                                          >  7:236468/1‑139 (MQ=255)
                                              aTTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCg                                                                          <  3:189883/139‑1 (MQ=255)
                                                    gTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                                    <  6:71995/139‑1 (MQ=255)
                                                     ttGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                                    <  4:217998/138‑1 (MQ=255)
                                                     ttGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                                    <  7:131730/138‑1 (MQ=255)
                                                     ttGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGGTGGGGGCGGTTTTTGCATGGCGGTTcac                                                                    <  5:261930/138‑1 (MQ=255)
                                                       ggCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAgt                                                                 >  5:33154/1‑139 (MQ=255)
                                                       ggCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAgt                                                                 >  5:56931/1‑139 (MQ=255)
                                                            aaCCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTAttt                                                                                                     <  2:92594/98‑1 (MQ=255)
                                                            aaCCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTAttt                                                                                                     >  1:92594/1‑98 (MQ=255)
                                                                        tttGAAAAATGATGCCCTATTTTCTTTCCTTTACGGTCTTTGACAAGCTGTTGTAAAGAACGGGTTTGCAGGATATATAGGTATTTAGGGTGGGTGGCGGGTTTTGAATGTGGGTTCACAGTGTAAACTCCTTTTg                                                   >  2:284144/1‑136 (MQ=255)
                                                                        tttGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGaaa                                                >  6:155898/1‑139 (MQ=255)
                                                                        tttGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGaa                                                 <  1:237046/138‑1 (MQ=255)
                                                                        tttGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGaa                                                 <  5:145130/138‑1 (MQ=255)
                                                                        tttGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTATTGCATGGCGGTTCACAGTGTAAACTGCTTTTGa                                                  >  3:218902/1‑137 (MQ=255)
                                                                        tatGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTCACAAGCATTTGTAAGGCCCGATTTTTCAGCAGATATAGGTATTAAGGGTTGGGGGCGGTTTTTGCATGGCCGATCACAGTGTATACTGCTTTTGa                                                  >  6:242374/3‑137 (MQ=255)
                                                                                     gCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGGTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGt                                   <  2:108672/139‑1 (MQ=255)
                                                                                           ttttCTTTCCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTg                             <  7:107233/139‑1 (MQ=255)
                                                                                                  tcCTTTACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCt                      >  8:222989/1‑139 (MQ=255)
                                                                                                     tttACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                                    <  4:156227/90‑1 (MQ=255)
                                                                                                     tttACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTcac                                                                    >  3:156227/1‑90 (MQ=255)
                                                                                                       tACCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTATTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCg                  <  4:218902/138‑1 (MQ=255)
                                                                                                        aCCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGt                 <  5:155898/138‑1 (MQ=255)
                                                                                                        aCCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGt                 <  7:250755/138‑1 (MQ=255)
                                                                                                        aCCTTCTTTGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTGTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGt                 <  5:233010/138‑1 (MQ=255)
                                                                                                              tttGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGTGGCGGTTTTTTAATGGCGTTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGATGGCTGCAGATCTTACGTATATg            >  1:78446/1‑137 (MQ=255)
                                                                                                              tttGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCa          >  4:80482/1‑139 (MQ=255)
                                                                                                                tGACAAGCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCa          <  7:222989/137‑1 (MQ=255)
                                                                                                                      gCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCATGgtcgtc  <  6:33154/139‑1 (MQ=255)
                                                                                                                      gCTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCATGgtcgtc  <  6:56931/139‑1 (MQ=255)
                                                                                                                       cTTTTGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCATGgtcgtc  >  3:232360/1‑138 (MQ=255)
                                                                                                                          ttGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCa          <  2:141406/127‑1 (MQ=255)
                                                                                                                          ttGTAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCa          >  1:141406/1‑127 (MQ=255)
                                                                                                                             |                                                                                                                                   
GGAGGCATGTGAAAATAGTAGTCTCCGATTGAGATAAGCACGTCCGATTGCTGTTGGCTGAACCTTGGATTTTTTGAAAAATGAAGCCCTATTTTCTTTCCTTTACCTTCTTTGACAAGCTTTTGCAAAGCCCGTTTTTTCAGCAGATATAGGTATTTAGGGTTGGGGGCGGTTTTTGCATGGCGGTTCACAGTGTAAACTGCTTTTGAAAGATTTTCTACGGTTGGCTGCAGATCTTTCGTATATGCATGGTCGTC  >  BS168_reference_genome/1474765‑1475021

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