Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,502 C→T T128I (ACT→ATT)  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,5020CT82.8% 52.0 / 7.3 30T128I (ACT→ATT) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (16/8);  total (19/11)
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.55e-01

CTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTA  >  NC_000913/3486259‑3486751
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
gtgtggaagtgatgaaaaaacaCGAAGAGGATAAAGAAATGTTCTGCCCCTATCTGAAGCAGGTTGATTTTATTGGCTACCTGGGCCTTGTTAAAGGGGGCCAGGAACTTAGCGAGGGGTTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCag                                                                                                                                                                                                                                                      <  1:47092/228‑1 (MQ=255)
    ggctgtgctaaacAAAGAAGAAGAGGGAAAAGAAAGAACCCTCCCCAATCTGAATCGGGGTGATTTTATGGGCAAACTGGTCCTGTTTAAAGGGGGCCAGAAACGTAGCGATTGGTTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAg                                                                                                                                                                                                                                                <  1:72132/239‑1 (MQ=255)
         tGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGAGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAACACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGc                                                                                                                                                                                                                                           >  1:384309/1‑251 (MQ=255)
                                                 cTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGCGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGcac                                                                                                                                                                                                      >  1:99116/1‑248 (MQ=255)
                                                  tATCTGAATCAGGGTGATTTTATTGACGAACTGGGCCTGTTTGCAGCGAGCCAGGAACGTAGCGCATGGGTACGTGCGAGAACCGCCTGTGAAGTGGCTGACATTTAGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGACTGCACGACTTTTCAGGGAAAGAGGGCAACCTGGCGTTCCTCGACGTGGCGGGCCGCATTGCCCAGAc                                                                                                                                                                                                  >  2:74499/1‑251 (MQ=255)
                                                  tATCTGAATCAGGGGTATTTTATTGACGAACTAGGCCTGTTTGAAGAGGGCCAGGAACATAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCACTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCTTCTGCAAGTCACTTAAGAGAAAGTGTGCAACCTGGCGTTCCACACCTTGACGGGCCGCTTTGCACAGAc                                                                                                                                                                                                  >  2:185089/1‑251 (MQ=255)
                                                     cTGAATCAGGGAGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGACAATTCGTACAAAAAAGTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCATCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACCTGACGGGCCGCATTGCACAGActct                                                                                                                                                                                               >  2:15138/1‑251 (MQ=255)
                                                       gAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTctgc                                                                                                                                                                                             >  1:155376/1‑251 (MQ=255)
                                                           cAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCCCAGACTCTGCTGa                                                                                                                                                                                          >  2:166405/1‑250 (MQ=255)
                                                                   ttttATTGGCGAACAGGGCCTGTTTGAAGAGGGCCAGGAAAGTAGCGCACGGGTACGTGCGCAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTCAACCCGGCCATTATGATGCGTTTGTCTGCACAGATGTCGCGTCGTCTGCAAGTCATTTCAGAGAAAGCGGGCAACCTGGCGTTCCTCGACGTGACGTGCCGCATTGCACCGACTCTGCTGAATCTGGCa                                                                                                                                                                                  >  2:180248/1‑250 (MQ=255)
                                                                             gcaCTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGAGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGACTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCCCAGCCTCTGCTGAATCTGGCAAAACAAcc                                                                                                                                                                          >  2:219255/3‑248 (MQ=255)
                                                                                                      aGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATg                                                                                                                                              >  1:359880/1‑251 (MQ=255)
                                                                                                         aaccgtggctaagggatacgggcaaaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCa                                                                                                                                            <  1:52618/227‑1 (MQ=255)
                                                                                                               ttagttgtgtaaCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTCCAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCAGTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGTCCGCATTGCACAGACTCTGCTGAATCTGGCAAATCAACCAGACGCTATGAATCCCCCGGACGGTATGCAAATCa                                                                                                                                       <  1:42559/239‑1 (MQ=255)
                                                                                                                   aTGGGGACGTGCGAAAACCGCCTGAGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTCGTCTGCACAGCTGACGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAAACTGGCGTTTCTCGACGTGACGTGCCGCATGGCACAGCCTCTTCTGAATCTGGCAAAACAACCAGACGCTATGACTAACCCAGATCGTATGCAAATCAAAAtt                                                                                                                                  >  2:304491/1‑250 (MQ=255)
                                                                                                                                  aaaccccccggtaaggtgcgtgaaTTTCCTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAAtt                                                                                                                   <  1:32497/229‑1 (MQ=255)
                                                                                                                                      cGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTc                                                                                                               >  1:127756/1‑251 (MQ=255)
                                                                                                                                        ccTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGGAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCCGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGCATCTGGCAAAACAACCAGAAGCTATGACTCACCCGGCAGGTATGCAACTCAAAATTACCCGTCAGGAAATTGGTc                                                                                                               >  2:368909/1‑249 (MQ=255)
                                                                                                                                         cTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActct                                                                                                                                                                                               >  1:128160/1‑168 (MQ=255)
                                                                                                                                         cTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActct                                                                                                                                                                                               <  2:128160/168‑1 (MQ=255)
                                                                                                                                             gAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGt                                                                                                        >  1:335908/1‑251 (MQ=255)
                                                                                                                                               aGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGCTGGCGCGTCGTCTGCAAGTCATTTCAGATAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGTCGCTGCTGAATCTGGCAAATCAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCg                                                                                                      <  1:178594/251‑1 (MQ=255)
                                                                                                                                                   gCTGAAATTTCGTACAAAATATTTCGACAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGATCCTCGACGTGACGGGCCGCATTGCACAGACACTGCTGACTCTGGCCAAACAACCAGACGCTATGACTCACCCGGCAGGGAACCAAAACACAATTACCCGTCTGGAACTTGGGCAGAGTGTCggttg                                                                                                  >  2:249015/1‑248 (MQ=255)
                                                                                                                                                                    aaaaTTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGCCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTgg                                                                                 <  1:178388/251‑1 (MQ=255)
                                                                                                                                                                                                              gCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGTTCAGAGTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCg                                       <  1:364460/251‑1 (MQ=255)
                                                                                                                                                                                                                             gATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAAcc                         >  1:152366/1‑250 (MQ=255)
                                                                                                                                                                                                                              tgggCGCGTCGTCTTCAAGTCATGTCAGAGAAAGTGGTCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTaaaa                           <  2:257575/245‑1 (MQ=255)
                                                                                                                                                                                                                              aTGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTaaaa                           >  1:257575/1‑247 (MQ=255)
                                                                                                                                                                                                                               tGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtc                      >  2:104509/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgtcgtCTGCAAGTCATTTCAGAGAAAGTGGGCAACCGGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGTTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATcgtcgt                 <  2:155376/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGCATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGAGCGGAACCTGATCTCCGCCCACGGTGAAACCATCGTCGTTTACGGCACTCGTTa  >  2:123131/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
CTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTA  >  NC_000913/3486259‑3486751

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