Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,157 A→C E13A (GAA→GCA)  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,1570AC97.9% 135.4 / ‑4.0 47E13A (GAA→GCA) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (24/22);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

TTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAA  >  NC_000913/3486019‑3486296
                                                                                                                                          |                                                                                                                                           
ttGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGCAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGCCAAACCGCAAACAGACCCGACTCTCTCATgg                                                                                                                                         >  1:1003743/1‑143 (MQ=255)
   aTAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGTAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTc                                                                                                                                      <  1:738374/143‑1 (MQ=255)
   aTAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCCTGGtac                                                                                                                                      >  2:505174/1‑141 (MQ=255)
     aGCCCCTTCCCAGGGAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCtt                                                                                                                                    <  2:130584/143‑1 (MQ=255)
        cccTTGCAAGGTGGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAATCCGCAAACAGACCCGACTCTCGCATGGTTCTTGtc                                                                                                                                 <  2:779892/143‑1 (MQ=255)
        cccTTCCCAGGGGGCGGGAAGCATATTTCGGCAATACAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGtc                                                                                                                                 <  2:333014/143‑1 (MQ=255)
             cccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCAtt                                                                                                                            >  2:295356/1‑143 (MQ=255)
             cccAGGTAGCAGGAAGCATATTTCCGCAATCCAGAGACAGCGGCGTTATCTGCCTCTGGAGAAAGCTTATAACAGAGGATAACCCCCCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTTTCTCAtt                                                                                                                            <  2:662537/143‑1 (MQ=255)
               cAGGTAGCGGGAAGCCTATTTCGCGAATCCAGGGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGc                                                                                                                          <  2:807825/143‑1 (MQ=255)
                      cGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTc                                                                                                                   <  1:746897/143‑1 (MQ=255)
                       gggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTTTTTGccccctcta                                                                                                                  >  1:172075/1‑136 (MQ=255)
                          aatgatatTTCGGTAATCCAGAGACAGCGGCGTTATTTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATaa                                                                                                               <  2:375634/139‑1 (MQ=255)
                             cATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCACAAGTa                                                                                                            >  2:794505/1‑143 (MQ=255)
                              atatTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTAc                                                                                                           <  1:405791/143‑1 (MQ=255)
                                  ttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGTTTGTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCAt                                                                                                       <  2:366360/143‑1 (MQ=255)
                                      gCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCaa                                                                                                   >  1:543345/1‑143 (MQ=255)
                                            cAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATCGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGCTCCTGTCCCAGTGCAACACTCAGAAGTCCTCAGTCACGAGCAc                                                                                             >  2:932883/1‑143 (MQ=255)
                                             agagACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACACGACTCTCGCATGGTTCTTGTCTCATTGCCACATCAATAAGGACCCATACAAGACCACg                                                                                            >  2:818854/1‑143 (MQ=255)
                                              gagaCAGCGGCGTTATCTGGCTCAGGATAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGc                                                                                           <  2:948866/143‑1 (MQ=255)
                                                              cTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGa                                                                           >  1:269835/1‑143 (MQ=255)
                                                               tGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACGCATCCAAGAGCACGCTTATTCACCAGGGTGaa                                                                          <  2:542959/143‑1 (MQ=255)
                                                                    gtgtcgaAAGCTTATCACAGAGGATACCCGCGCAGGGCGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                     <  2:653418/138‑1 (MQ=255)
                                                                       gagaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCgca                                                                  >  2:342826/1‑141 (MQ=255)
                                                                        agaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCCCGCTTATTCACCCAGGTGAAAAAGCggga                                                                 >  1:148735/1‑141 (MQ=255)
                                                                                   aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtacta                                                      >  1:416687/1‑143 (MQ=255)
                                                                                    aCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtactac                                                     >  1:166626/1‑143 (MQ=255)
                                                                                     ccgagGACAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACa                                                    <  1:607592/141‑1 (MQ=255)
                                                                                          gATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCCTCCAAGAGCCCGCTTATTCACCCGCGCGACCACGCGGAAACGCTGTCCTCCATCGtt                                               >  1:69780/1‑143 (MQ=255)
                                                                                          gATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGGACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAGAAGCGGAAACGCTGTACTACATCGtg                                               >  1:127654/1‑142 (MQ=255)
                                                                                           aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCGGTACTACATCGTTa                                              >  2:439432/1‑143 (MQ=255)
                                                                                            tAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACACTCACAATTACCATTCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTaa                                             <  2:736174/143‑1 (MQ=255)
                                                                                                 gcgcGTGGTGCTTGGCAAACCGCAAACAGACCCCACTCTCGCATGGGTCTTGTCTCATTGCCACAGTCAGAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGct                                        <  2:50300/143‑1 (MQ=255)
                                                                                                           cTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCCTTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGT‑AAAAGGCGCAAACGCTGTCCTACATCGCTAAAGGCTCTGTGCCCGTg                             >  2:80444/1‑143 (MQ=255)
                                                                                                            ttGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACAGTCATAAGTACCCATCCAAGTGCACGCCGATTCACCACGGTGAACAAGCGGAAACGCGGTACTAGATCGTTAGAGGCTCTGTGGCAggg                             >  1:38761/1‑141 (MQ=255)
                                                                                                             tGGCAAACCGCAAACAGACCCGACTCTCGCATGCTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                            <  1:541697/143‑1 (MQ=255)
                                                                                                                    ccGCAAGCAGACCCGACTCTCGCATGGTTCTTGTCTTATTGCCACATTACTAAGTAACCTTCCCAGAGCGCGCTTATTCACCAGGGTGAAAAAGCGGAACAGTTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaa                     <  2:686008/143‑1 (MQ=255)
                                                                                                                     ggCAACCAGGCCCGCCTCTCGCATGGTTCGTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaaa                    <  2:34609/142‑1 (MQ=255)
                                                                                                                     cGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaaa                    <  2:543345/143‑1 (MQ=255)
                                                                                                                      gCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAg                   >  1:271710/1‑143 (MQ=255)
                                                                                                                          aCAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGa               >  1:509786/1‑143 (MQ=255)
                                                                                                                            aGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAag             >  1:826667/1‑143 (MQ=255)
                                                                                                                               cccGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAgg          <  2:416687/143‑1 (MQ=255)
                                                                                                                                ccGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAggg         >  1:556912/1‑143 (MQ=255)
                                                                                                                                ccGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATACAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAggg         >  2:732579/1‑143 (MQ=255)
                                                                                                                                ccGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATACAAGAGCACGCCTATTCACCAGGGAGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGACTCCGTGCCAGTGCTGATCAAAGCCGAAGAggg         >  2:247491/1‑143 (MQ=255)
                                                                                                                                  gACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGCGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGCTCAAAGACGAAGAGGGtg       >  2:347046/1‑142 (MQ=255)
                                                                                                                                       tcGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGaa  <  1:918/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAA  >  NC_000913/3486019‑3486296

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