Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,774 T→A 9.3% intergenic (+24/‑15) crp → / → yhfK cAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,7740TA 9.3% 140.2 / 7.8 55intergenic (+24/‑15)crp/yhfKcAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base A (3/2);  total (31/24)
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.65e-01

AGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCT  >  NC_000913/3486629‑3486918
                                                                                                                                                 |                                                                                                                                                 
aGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGt                                                                                                                                                <  2:549552/149‑1 (MQ=255)
aGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGt                                                                                                                                                <  1:278265/149‑1 (MQ=255)
aGGAAATTGGTCAGATTGTCGGCTGTGCTCGTGAAACCGTGGGGCGCATGCCGAAGATGCTGGAAGATCAGAACCTGATCCCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGt                                                                                                                                                <  1:549651/149‑1 (MQ=255)
  gAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAg                                                                                                                                              >  1:716499/1‑149 (MQ=255)
         gTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACACGCCa                                                                                                                                       <  2:668840/149‑1 (MQ=255)
         gTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCCGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACACGCCa                                                                                                                                       <  2:55352/149‑1 (MQ=255)
               ttGTCGGCTGTCCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGt                                                                                                                                 <  1:149384/149‑1 (MQ=255)
                     gcatttttcGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGACCCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCGTTTTGTTTcccc                                                                                                                           >  2:97076/8‑149 (MQ=255)
                         ttCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATGCCGTCTGTGTGGCGCGTTTCCTGGTAGCGCGCCCTTTTGTTTCCCCCGAt                                                                                                                       >  1:300176/1‑149 (MQ=255)
                           ctcGTGAAACCGTGGGACGCATTCTGAAAATGCTCGAAGATCAGAACCTGATCGCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAtgt                                                                                                                     >  1:430998/1‑149 (MQ=255)
                            tcGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAtgtg                                                                                                                    >  1:269071/1‑149 (MQ=255)
                            tcGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAtgtg                                                                                                                    >  1:561478/1‑149 (MQ=255)
                             cGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAGCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTgg                                                                                                                   <  2:536907/149‑1 (MQ=255)
                               tGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGgcg                                                                                                                 <  1:553841/149‑1 (MQ=255)
                               tGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGACCAGAACCAGATCTCCGCACACGGTAACAAACTAGGCTGGTACGGCGATGGTTTGTCCCTTCGGAGTGGCGCGTTACCTTGTAGCGCGCCATTTTGTTTCCCCCGATGTGgcg                                                                                                                 <  1:573667/149‑1 (MQ=255)
                                      gTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAACCCGTCGTCGTTTACGGCACTCGTTAATCCCGTCGGAATGGGGCGTGACCGGGTAGGGCGCCCTTTTGTTTCCCCCGATGTGGCGCAGACTg                                                                                                          >  1:180046/1‑149 (MQ=255)
                                        gggACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                                                                        <  2:130450/149‑1 (MQ=255)
                                          gACGCATTCTGCAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAttt                                                                                                      <  2:114112/149‑1 (MQ=255)
                                          gACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAttt                                                                                                      <  1:170635/149‑1 (MQ=255)
                                            cGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTAt                                                                                                    <  2:499696/149‑1 (MQ=255)
                                             gCATTCTGGAGATGCTGGAAGATCAGAACCTTATCTGCGCACACGGTAAACCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATc                                                                                                   <  2:21289/149‑1 (MQ=255)
                                                       gATGCCGGAAGATCAGAACCTGATCTACGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGatat                                                                                         <  1:322529/149‑1 (MQ=255)
                                                        aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATc                                                                                        <  1:716545/149‑1 (MQ=255)
                                                        aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATc                                                                                        >  2:17505/1‑149 (MQ=255)
                                                          gCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCaa                                                                                      >  1:169268/1‑149 (MQ=255)
                                                             ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑Ta                                                                                  >  2:574526/1‑149 (MQ=255)
                                                               aaGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATg                                                                                <  2:32972/149‑1 (MQ=255)
                                                                ataTCAGACCCGGTTCCCCGCACAGGGTAAAACCCTCGTGGTTTACGGCCCTCGTTAATCCCGCCGGAGTGGCGCGTTACCGGGTAGCTCGCCATTTTTTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGc                                                                               <  1:199175/147‑1 (MQ=255)
                                                                aGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGc                                                                               >  2:554266/1‑149 (MQ=255)
                                                                  aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCAc                                                                             >  2:12686/1‑149 (MQ=255)
                                                                     aGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCCAC‑TATGCACTTc                                                                          >  2:132401/1‑149 (MQ=255)
                                                                      gAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGGCTGATTTATCACCCCGATATCAAC‑TATGCGCTTCg                                                                         >  2:580189/1‑149 (MQ=255)
                                                                      gAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGCCTGATTTATCACCCCGATATCAAC‑TATGCACTTCg                                                                         >  2:228019/1‑149 (MQ=255)
                                                                                tcCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTCCCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGAGATCAAC‑TATGCACTTCGACAAACGCTg                                                               >  1:626690/1‑149 (MQ=255)
                                                                                     acacGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGGCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCt                                                          >  2:34107/1‑149 (MQ=255)
                                                                                                 cATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTATCTGGTAGCGGGATATTTTGTTTCTACCGATGTGGCGCAGAGTGATTTATCACTCCGATATCAAC‑TATGCACTCCGACAAACCCTGGTGCTATGTTTGCCCGt                                              >  2:573852/1‑149 (MQ=255)
                                                                                                    cgtcgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGc                                           >  2:87147/1‑149 (MQ=255)
                                                                                                    cgtcgtTTACGGCACTCGTTAATCCCGTCGGAGTGACGCGGTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGCCTGATTTATCACCCCGCTATCAAC‑TATGCCCTTCGCCAAACGCTGGTGCTATGTTTGCCCGTGGc                                           >  2:367696/1‑149 (MQ=255)
                                                                                                      tcgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCg                                         >  2:769753/1‑149 (MQ=255)
                                                                                                            aCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGt                                   >  2:701596/1‑149 (MQ=255)
                                                                                                              ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGCTTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTa                                 >  2:274399/1‑149 (MQ=255)
                                                                                                               gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTaa                                <  2:691366/149‑1 (MQ=255)
                                                                                                                 gcTCGTTAATCCCGGCGGAGTGGCGCGTTACCTGGTGGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATg                              <  1:213790/148‑1 (MQ=255)
                                                                                                                      ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCCCTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTgg                         >  2:638922/1‑149 (MQ=255)
                                                                                                                          tCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACGTCGACAAACGCTGGTGCCATGTTTGCCCGTGGCCGGTGGGTTAATGCTTGGCGaa                     >  2:379478/1‑149 (MQ=255)
                                                                                                                              gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCCTTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTGTGCCCGTGGCCGCTGGGTTAATGCTTGGCGAATTAc                 >  2:295637/1‑149 (MQ=255)
                                                                                                                              gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTAc                 >  2:92134/1‑149 (MQ=255)
                                                                                                                                   aatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTc            >  1:410683/3‑149 (MQ=255)
                                                                                                                                   aatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCCATGTTTGCCCGGTGCCGGTGGGTTAATGCTTGCCGAATTACGATTc            >  1:262283/3‑149 (MQ=255)
                                                                                                                                    gTGGCGCGTTACCTGGTAGCGCGCCATTTGGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TCTGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCg           <  1:85391/149‑1 (MQ=255)
                                                                                                                                      ggCGCGTTACCTGGGAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGt         <  1:203276/149‑1 (MQ=255)
                                                                                                                                       gcgcgcTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAACGCGTGGCGAATTACGATTCGGTc        >  2:784786/1‑149 (MQ=255)
                                                                                                                                         gcgTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTg      <  1:12686/149‑1 (MQ=255)
                                                                                                                                            ttACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGctc   <  1:554266/149‑1 (MQ=255)
                                                                                                                                             tACCTGGTAGCGCGCAATCTGGATTCCCCCGATGTGGCGCAATNCTATGGNTCACAAGAATGCAAACAGGTGC‑CATCGACAAACGTTGGTGCTACGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGctct  >  2:166425/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCT  >  NC_000913/3486629‑3486918

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