Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,538 T→C V140A (GTG→GCG)  crp → DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,5380TC80.0% 46.1 / 10.1 30V140A (GTG→GCG) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (13/11);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

GTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCT  >  NC_000913/3486322‑3486776
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
gTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTg                                                                                                                                                                                                              >  1:90969/1‑251 (MQ=255)
     tttATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCaa                                                                                                                                                                                                          >  2:50774/1‑250 (MQ=255)
               aaCTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCAAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGAc                                                                                                                                                                                               >  1:17749/1‑251 (MQ=255)
                     gCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGCCTCTGCTGAATCTGGCAAAACAACCAGACGCTATg                                                                                                                                                                                         >  1:11239/1‑251 (MQ=255)
                                 ggggCCAGGAACGTCGCGCGTGGTTACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGAc                                                                                                                                                                             <  1:45723/250‑1 (MQ=255)
                                        ggAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGCCGGGCCGCATTGCACAGACTCTGCTGCATCTGGCAAAACAACCAGACGCTATGACTCACCCGAACGGTATGc                                                                                                                                                                      >  2:231451/1‑251 (MQ=255)
                                             gTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc                                                                                                                                                                 >  2:86901/1‑251 (MQ=255)
                                                      gggTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCgg                                                                                                                                                                               >  1:243707/1‑228 (MQ=255)
                                                          aCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGGCGGTATGCAAATCAAAATTACCCGTc                                                                                                                                                    >  1:223535/1‑251 (MQ=255)
                                                            gtacGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAACATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCCCAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTCCTGAATCTGGCACAACAACAAGAAGATATGACTCACCCGGCCGGTATGCCAAACAAAATTACCCGTccg                                                                                                                                                  >  2:277681/4‑249 (MQ=255)
                                                               aaaaatccGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCTTCTGCAAGTCATTTCAGAGAAAGTGGTCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGaa                                                                                                                                               <  1:169047/245‑1 (MQ=255)
                                                                  aaaCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGTACATTCTGAGGCGTTTGTCTGCACAGATGGCGCGTCGTCTGAAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTACTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCACACCACCCAGACGCTATAACACACCCGGCCGGTATGCAAATCCAAATTCACCGACAGGaactt                                                                                                                                            >  2:165638/1‑248 (MQ=255)
                                                                                       aaaTTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTtct                                                                                                                       <  1:59383/251‑1 (MQ=255)
                                                                                                              cAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACACCCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCt                                                                                                >  1:288633/1‑251 (MQ=255)
                                                                                                                   gATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCACAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGCGGGCCGCATTCTgc                                                                                              >  2:266400/1‑247 (MQ=255)
                                                                                                                         ggTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTgg                                                                                     <  1:104084/251‑1 (MQ=255)
                                                                                                                          gTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTCCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGa                                                                                    <  1:86901/251‑1 (MQ=255)
                                                                                                                           tAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGAATGCAAATCAAAATTACCCGTCAGGAAATTGGACAGATTGTCGGCTGTTCTCGTGAAACCGTGGGAAGCATTCTGAAGATGCTGGa                                                                                    >  2:62833/1‑250 (MQ=255)
                                                                                                                                                    gtcggCACAGATGGTGAGTTGTCTGCAAGTCATTACAGAGAAATTGGGCACCCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCGCCCGGACTGTATGCAAATCAAAATTCCCCGTCAGGAAATTGGGCAGATTGTCGGCTGTTTTCGTGAAACCGTGGGACGCATTCTGAAGATGCTTGAAGATCAGAACCTGATCTCCGCACACg                                                          <  2:16649/247‑1 (MQ=255)
                                                                                                                                                        gCACAGATGGCGCGCCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTa                                                       <  1:278638/250‑1 (MQ=255)
                                                                                                                                                          acaGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGGAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTaaa                                                     >  2:53612/1‑250 (MQ=255)
                                                                                                                                                           caGATGGCGCGTCGTCGGAGAATCATTTCAGAGAAAGTGGGCAACCTGGCTTTCCCGGACGCGTCGGGCCGCATTGCTCAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCTTCAGGAAATTGGTAAGATTGTCGGCTGTTCTCGTGAATCCGTGGGACGCATTCTGAAGATGCTGGAAGAACAGAACCTGATCTCCGAACACGGGAAAAc                                                   <  2:71872/251‑1 (MQ=255)
                                                                                                                                                                       cgtcTGCAAGTCATTTCAGAGAAAGTGGACAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGCCGGTCTGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGttt                                        >  1:295853/1‑250 (MQ=255)
                                                                                                                                                                            gCAAGTCATTTCAGAGAAAGTGGGCAACCGGGCGTTCCGCGACGCGACGGGCCGCAGTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACgg                                    <  1:165638/249‑1 (MQ=255)
                                                                                                                                                                             cAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCAc                                 <  1:38587/251‑1 (MQ=255)
                                                                                                                                                                                      tCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTaa                          <  1:50774/249‑1 (MQ=255)
                                                                                                                                                                                       cAGAGAAAGTGGGCAACCAGGCGTTCCGCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATcc                       <  2:66824/251‑1 (MQ=255)
                                                                                                                                                                                                     aaCCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGCCGTTTACGGCACTCGTTAGTCCCGACGGAGTGgcgc         >  1:194374/1‑251 (MQ=255)
                                                                                                                                                                                                             gTTCCTCGCCGCGACCGTCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt  <  2:84735/250‑1 (MQ=255)
                                                                                                                                                                                                             gTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt  >  1:84735/1‑250 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
GTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCT  >  NC_000913/3486322‑3486776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: