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% 138.9 / NA 48E13A (GAA→GCA) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (30/17);  total (31/17)
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.99e-01

AGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAG  >  NC_000913/3486017‑3486294
                                                                                                                                             |                                                                                                                                         
aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCACACCGCAACCAGACCCGACTCGCGCAt                                                                                                                                         >  2:1009041/1‑143 (MQ=255)
 gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATg                                                                                                                                        <  1:543583/143‑1 (MQ=255)
          cccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTCATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCACGGTTCTTGtc                                                                                                                               >  2:475668/1‑143 (MQ=255)
               ccccGGTAGCGGGAGGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCAtt                                                                                                                          <  2:295919/143‑1 (MQ=255)
               cccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCAtt                                                                                                                          >  2:1089215/1‑143 (MQ=255)
                      aGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTCATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACAt                                                                                                                   <  1:28468/143‑1 (MQ=255)
                               cATATTTCGGCAAGCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTa                                                                                                          <  1:132539/143‑1 (MQ=255)
                                atatTTCGGCAATCCAGAGACAGCGGCGTTATATGGCTCTGGAGAAAGCTTATAGCAGAGGATAACCACGCATGTTGCTGGGCAAACTGCGGACAGACTCGACTCACGCGGGGTTCTGGTCGCCGTGCCACAGTCATATGtcc                                                                                                         >  2:193378/1‑141 (MQ=255)
                                  atTTCGGCAACCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTAccc                                                                                                       <  1:710368/143‑1 (MQ=255)
                                   tttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCa                                                                                                      >  1:65141/1‑143 (MQ=255)
                                   tttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCa                                                                                                      >  1:973330/1‑143 (MQ=255)
                                      cGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCCCAAAGAGACCCGACTCTCGCATGATTCTTGTCTCAGTGCCACATTAAAAATTACCCATca                                                                                                   >  2:799711/1‑142 (MQ=255)
                                        gCAATCCAGACACAGCGGCGTTATCTGGCCCTGGAGAAAGCTTATAACAGAGGATCACCGCGCATGGTGCTTACCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCaa                                                                                                 <  1:488795/143‑1 (MQ=255)
                                            tCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGc                                                                                             <  1:170518/143‑1 (MQ=255)
                                                gagaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGc                                                                                         <  2:451022/143‑1 (MQ=255)
                                                   ccAGCGGGGTTATCTGGCTCTGGAGTAAGTGTATAACAGAGGATAACCGCTCATGGTGCTTGGCAAACCTCAAACAGACCCGACTCTCGCCTGGTTCTTGTCTCATTGCCACATTCATAAGTTCCCATCCAAGAGCACGCTTa                                                                                      <  2:1055794/142‑1 (MQ=255)
                                                     cgcgtcgTTATCTGGCACTTGCGAAAGCCTGTACCACAGCATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGGGTCTCGCACGGTTCTTGTCTCATTGCCACTTTCATAAGTACCCATCCAAGAGCACGCTTAtt                                                                                    <  2:491775/142‑1 (MQ=255)
                                                       cggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCa                                                                                  >  1:729648/1‑143 (MQ=255)
                                                       cggcgTGATCTGGCTGTGGAGAAAGCTCATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCa                                                                                  <  2:789817/143‑1 (MQ=255)
                                                        ggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCAGCACCACTCCCTTTTCGTGTCTCCATGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCAc                                                                                 >  1:261256/1‑143 (MQ=255)
                                                        ggcgTTATCTGGCTCTGGAGAAAGCATATAACAGAGGATAACCGCGCATGGTAAATTGCAAACCGCAAAATGACCGGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTAGCCATCCAAGAGCACGCTTATTCAc                                                                                 >  2:814189/1‑143 (MQ=255)
                                                         gcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCAcc                                                                                >  1:760447/1‑143 (MQ=255)
                                                          cggtATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCa                                                                               <  2:460874/140‑1 (MQ=255)
                                                          cgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCACCAGTGCCCATCCCAGAGCACGCTTCTTCAccc                                                                               >  2:232755/1‑142 (MQ=255)
                                                             tATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAggg                                                                            <  1:865112/143‑1 (MQ=255)
                                                                 tGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGAATGGTTATTGTCACATTGCCACATATATAAATAACAACACAAGAGCACTATTATTCACCAGAGTGaa                                                                        >  2:768350/1‑143 (MQ=255)
                                                                   gCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCGAAGAGCACGCTTATTGACAAGGGTGaaaa                                                                      >  1:938012/1‑143 (MQ=255)
                                                                     tctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAg                                                                    >  1:447709/1‑143 (MQ=255)
                                                                      ctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                   >  1:841559/1‑143 (MQ=255)
                                                                       tGGAGCAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCCAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACGTTCATCACTACCCATCCAAGCGCACGCTTCTCCACCAGGGTGACACAGCg                                                                  >  2:197663/1‑143 (MQ=255)
                                                                         gggaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGa                                                                <  2:726070/141‑1 (MQ=255)
                                                                              aGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACAGAAAACAGAGACGACTCTCGCATGGTTCTTGTGTCATTGCCAAATTCATAAGTACACATACAAGAGCACGCTTATTCACCAGGGTGAAAAAGAGGAAACGc                                                           >  2:162636/1‑143 (MQ=255)
                                                                                 ttATAACAGAGGATAACCGCGCATGGTGCTTGGCTAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAACGCGGAAACGCCGt                                                        >  2:759554/1‑143 (MQ=255)
                                                                                 ttATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGCGCACGCTTATTCACCCGGGCGACACAGCCGCAACGCTGt                                                        >  2:542573/1‑143 (MQ=255)
                                                                                   ataACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtac                                                      <  2:841559/143‑1 (MQ=255)
                                                                                            gATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCCTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGACACGCTGTACTACATcacg                                             >  1:511365/1‑140 (MQ=255)
                                                                                               aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGATCTTGTCTCATTGCAACATTTATAAGTACCCCACCAAGAGCACCCTTATTCACCAGCGTGAAAAACCGGAAACGCTTTCCTAAATCgcgcaa                                          >  2:660304/1‑138 (MQ=255)
                                                                                                    cgcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGctc                                     >  1:823198/1‑143 (MQ=255)
                                                                                                             cTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                            >  1:533798/1‑143 (MQ=255)
                                                                                                                  cAAACCGCAAACAGACCCGACCTCAGCTTGGTTCTTGTCTCATTGCCACATTCATGAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGTAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGa                       <  1:207248/143‑1 (MQ=255)
                                                                                                                       cGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCCTTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTGCTACATCGTTAAGGGCTCAGCAGCCGAGGTGATCaaa                  >  2:571724/1‑143 (MQ=255)
                                                                                                                       cGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaaa                  <  2:1028493/143‑1 (MQ=255)
                                                                                                                         cAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGa                >  2:1147346/1‑143 (MQ=255)
                                                                                                                            aCAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGa             >  2:300192/1‑143 (MQ=255)
                                                                                                                                aCCCGACTCTCGCGTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCCCGCTTATTCACCAGGGTGACAAAGCGGAAACGCTGTACCACCTCGTTAAAGGCTCTGGGGCAGCGCTGATCAAAGACGAagag         >  2:302079/1‑143 (MQ=255)
                                                                                                                                aCCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCACGAGCACGCTTGTTCACCAGCGTGACAAACCGGAAAAGCTACACTACGTCGTTAACGACACTTAGGGAGCGCTGATCAAAGCCAAagcg         >  2:287778/1‑141 (MQ=255)
                                                                                                                                aCCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACACATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAACGCTCTGTGGCAGTGCTGATCAAACACGAagag         >  1:1031897/1‑143 (MQ=255)
                                                                                                                                       tctcGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTCCTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAg  >  1:935997/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                         
AGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAG  >  NC_000913/3486017‑3486294

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: