Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,205 C→A T29K (ACG→AAG)  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,2050CA95.0% 48.6 / ‑4.7 20T29K (ACG→AAG) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (8/11);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTA  >  NC_000913/3485958‑3486446
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
tACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTAAAGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCCCATTCATAAGTACCCATCCAAGAGCAAGCt                                                                                                                                                                                                                                                 >  1:146356/1‑251 (MQ=255)
               atgaaCTGCATGTATGCAAAGGATGGCACATTCCCGGGCAGTACAGGTGATAGCCCCTGCCCAGGTAGCGGGAAGCATATGTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGa                                                                                                                                                                                                                                  <  2:140134/247‑1 (MQ=255)
                                                           agttgtcaGCCCTTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGATACAGCGGCGTTATCTGGCTCTGGAGAAAGCTGATAACAGAGTATAACCGCGCAGGTTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                      <  2:432404/244‑1 (MQ=255)
                                                             ttGATAGCACCTTCCCAGGAAGAGGGAAGCATATTTCAGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACAGCGAATGGTGCTTGGCAAACCGAAAACAGACACGACTCTCGAATGGTTCTAGGCTCATTGCCACATTAATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGGAAAAGAGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGAAGTg                                                                                                                                                                                    >  2:317740/1‑251 (MQ=255)
                                                                agtgcctccctctaggtgtggggAAGCATATTTGGGAAATGCAGAGACAGGGGAGTTACCTGACACGGGAGAAAGTTTATAACATATGATCACCGCTCATGGTTCTTTGCAAACCTAAAACAGACCTGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAATAGCAAGCTTTTTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACTTCGTTAAAGGCTCTGTGGCAGTGCTg                                                                                                                                                                                 <  2:146356/231‑1 (MQ=255)
                                                                      ccTTCCCAGGTAGAGGGAAGCATATTTCGTCAATACAGCGACCGCAGCGTTATCTGGCTATGGAGAAAGCTTATAACAGAGGATAACCGCACATGGTGCTTGGCAAACAGCAAACAGAACCGACTATCGAATGGTTCTTGTCTCATTGCCACAATAATAAGTAACCCTCCAAGAGCAAGCTTATTCACAGGGGAGAAAAAGCGGCCACGCTGAACAACATCGGTAACGGCTAGGGGGCAGGGCTGATCaa                                                                                                                                                                            >  2:396871/1‑250 (MQ=255)
                                                                               gTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCTTTATCTGGCTCTGGAGAAATCTTATAACAGAGGATACCCGCGCATGTTTCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACTCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                                                                                                  <  1:47367/251‑1 (MQ=255)
                                                                                   ggggAAGCATATTTCGGCAATCCAGAGACAGCGTCGTTATCGGGCTATGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGTGTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACTAAGAGGGTa                                                                                                                                                              <  2:257901/250‑1 (MQ=255)
                                                                                        agcttatgtCGGCAACCCAGAGACAGCGGCGTTAGCTGGCTCGGGAGAAAGCTTAAAACAGAGGATAACCGCTCATGGTTCTTGTCAAACCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGaa                                                                                                                                                         <  2:397750/243‑1 (MQ=255)
                                                                                             ttttCGTCAATCCAGAGACAGAGGCGTTAGCTGGCTTTGGAAAAAGCTTATAACAGAGGATACCCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTTCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGa                                                                                                                                                     <  2:148695/249‑1 (MQ=255)
                                                                                                              acatcggcgTTCTCGGGCTCTGGAGAAAGTGTAAAATAGGGGATAACCGGGTATGGTGCTTGGCAAACCGCAAACAGCGCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaagc                                                                                                                                   <  2:104163/247‑3 (MQ=255)
                                                                                                                                                   agagGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTTTTGGCGCACTGGGTCTGTTTGa                                                                                                >  2:88102/1‑248 (MQ=255)
                                                                                                                                                                    ggTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGCTCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGATTGAAGAGGGACAGGAACGAAgcg                                                                            >  2:455427/1‑251 (MQ=255)
                                                                                                                                                                      tGCTTGGCAAACCGCAACCAGACCCGACTCTCGAATGGTTATTGTCTCATTGCCACATTCATAAGTACCAATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAATCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTACTATCTGACTCAGGGTGATTTTATTGGCGAAATGGGCCTGTTTGCAGAGGGcagggacagtggcgcc                                                                          >  1:423873/1‑236 (MQ=255)
                                                                                                                                                                            gCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGt                                                                     >  1:447083/1‑250 (MQ=255)
                                                                                                                                                                                  cGCAACCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCg                                                              <  1:233183/251‑1 (MQ=255)
                                                                                                                                                                                     aaacagcaccaaCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGGAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAGTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGaaa                                                           <  2:355029/240‑1 (MQ=255)
                                                                                                                                                                                                              cTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAAttt                                  >  1:82226/1‑251 (MQ=255)
                                                                                                                                                                                                               ttGTCTAATTGCCACATTCATAAGTCCCCACCAAAGAGCAACTTTATTCACCAGGGTGAAAAATCGGAAACGCTGTATAACAGTGTTAAAGGCTTTGTGGCAGTGCTAATAAAAGACAAAGAGGGAAAAAAAATGTTCCCCCCCTACCTGAATCAGGTTGATTTTATTGCCGAACTGGGCCTGTTTGAAGAGGCCAAGGAACGTGGCCCTTGGGTACTTGCGAAAACCGCCTGTGAAGTGGCTGAAAttt                                  <  1:332670/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAAGAGCAAGTTTATCCACCAGGGTAAAAAGGGGGAACCGCTGTACTACTTCTTTAAAGGCTGTGTGGAATTGCTGATCAAAGAGGAAGAGGTTAAAAAAGTGTCCCTCTCTTACCTGAAGCGGGGTTATTTAATTGGGAATCGGGTCCTTTTATAAGGGGGCCAGGTGCGTGGGGCAGGGTTACGGCCGAAAACCTCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTa  <  1:184241/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
TACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTA  >  NC_000913/3485958‑3486446

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