Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,779 G→A 9.7% intergenic (+29/‑10) 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,7790GA 9.7% 128.7 / 6.4 52intergenic (+29/‑10)crp/yhfKcAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/19);  new base A (3/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-02

GAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCTTCTCCCTC  >  NC_000913/3486631‑3486926
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                >  2:17505/1‑149 (MQ=255)
                                                      aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATc                                                                                                <  1:716545/149‑1 (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)
                                                                                                                                                   aGCGCGCCATTTTGTTTCCGCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACCAACGCTGGTGCTATGTTTGACCCGGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCTTCTCCCTc  >  1:569132/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAC‑TATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCTTCTCCCTC  >  NC_000913/3486631‑3486926

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: ‑