Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,393 T→A 10.6% intergenic (‑54/‑30) ybjJ ← / → rcdA inner membrane protein YbjJ/DNA‑binding transcriptional regulator RcdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913887,3930TA10.6% 98.6 / 3.3 45intergenic (‑54/‑30)ybjJ/rcdAinner membrane protein YbjJ/DNA‑binding transcriptional regulator RcdA
Reads supporting (aligned to +/- strand):  ref base T (24/13);  new base A (3/2);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTG  >  NC_000913/887256‑887534
                                                                                                                                         |                                                                                                                                             
cATTAACAGGCCTGGCAAAAAGAAGAACAAAAACAGGGCCCGGCTTTTGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTac                                                                                                                                          >  1:777588/1‑143 (MQ=255)
 aTTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTaca                                                                                                                                         >  2:392211/1‑143 (MQ=255)
  ttAACAGGCCGGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                                                                        <  2:445490/143‑1 (MQ=255)
    aaCAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGCAGACAACAATAGCAAGGTTGTGTACAAATGAACACaa                                                                                                                                      >  2:339724/1‑143 (MQ=255)
    aaCACGACTGGAAAAAAGAAGAACATAAACAGCGCCCACGTTCGGAGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaa                                                                                                                                      <  1:464642/143‑1 (MQ=255)
      cgtgCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATGACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                                                                    <  1:192564/140‑3 (MQ=255)
              gCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATACGTTCACTATTGCggagag                                                                                                                            >  1:723858/1‑143 (MQ=255)
                aaaaaGCAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                                                                          <  2:312852/143‑1 (MQ=255)
                 aaaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATAACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACGCAATTGCAGAGAGAgg                                                                                                                         >  1:205725/1‑143 (MQ=255)
                     gaagcaCACAAACAGCGCGCCGGGTCGGCCTGTCAACGCATTACGAGAAGAATTGACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTGTCTCGGGTCACTCCACTTCCACACAATTGCAGAGAGAGGAAAt                                                                                                                     <  2:669417/143‑1 (MQ=255)
                      aagaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGCGAGCGGAGATg                                                                                                                    >  1:670085/1‑143 (MQ=255)
                         aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTGCACAATTGCAGAGAGAGGAAATGaaa                                                                                                                 >  1:159521/1‑141 (MQ=255)
                           cATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGCGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAAACAATTGCAGAGAGAGGAAATGACATg                                                                                                               >  1:664974/1‑143 (MQ=255)
                           cATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTGGTACACAATTTCATAGATATGAAATGACATg                                                                                                               >  1:563134/1‑143 (MQ=255)
                                  aGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCTCGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgcg                                                                                                        <  1:686599/143‑1 (MQ=255)
                                       ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGCGCGGAACTGACATGCGTCGCGATAAc                                                                                                   >  2:586498/1‑143 (MQ=255)
                                                    ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCACAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCGAACGACCCGCAACGGc                                                                                      >  2:236164/1‑143 (MQ=255)
                                                    ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                                                                      >  2:388440/1‑143 (MQ=255)
                                                     tcaaaccATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                                                                     <  1:444681/137‑1 (MQ=255)
                                                        aTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCCGAGAGGCGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                                                                  >  2:407264/1‑143 (MQ=255)
                                                           cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGCAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaa                                                                               >  2:160592/1‑143 (MQ=255)
                                                            aTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACGCAAGTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCGGCAACGGCGAGaaaaa                                                                              >  2:612229/1‑143 (MQ=255)
                                                            aTTACGTAAAGAATTTACGGTCATAGATCACGTAAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaa                                                                              <  1:16857/143‑1 (MQ=255)
                                                                 gTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTAt                                                                         >  2:646769/1‑143 (MQ=255)
                                                                    aagagTTTGCGGTCATAGAGCACGTCAGAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCa                                                                      <  2:191590/143‑1 (MQ=255)
                                                                    aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCa                                                                      >  2:502595/1‑143 (MQ=255)
                                                                              ggTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                                                            >  2:646156/1‑143 (MQ=255)
                                                                                     ggtCACGTCAAAATACGTAGCGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg                                                     <  2:246071/141‑1 (MQ=255)
                                                                                     gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATGTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg                                                     <  1:343711/143‑1 (MQ=255)
                                                                                           gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGACCACGCTGGAGGCGGTGAAAc                                               >  2:177680/1‑143 (MQ=255)
                                                                                               aaaTAAGAAGGGGGAAGACAACAGTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                           <  1:707283/143‑1 (MQ=255)
                                                                                               aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                           <  1:365281/143‑1 (MQ=255)
                                                                                               aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                           >  1:616355/1‑143 (MQ=255)
                                                                                                 aTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                                         >  2:205015/1‑143 (MQ=255)
                                                                                                      aagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCCACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATa                                    >  2:9238/1‑143 (MQ=255)
                                                                                                      aacagGGAAGACAACATTACCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATa                                    <  1:367745/140‑1 (MQ=255)
                                                                                                         acggAAGAGAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACCATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACAt                                 <  1:711440/141‑1 (MQ=255)
                                                                                                          tggAAGACAACATTAGCAAGGTTGTGTACGTGTTTACACAATTGCAGAGATAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATg                                <  1:111140/142‑1 (MQ=255)
                                                                                                          gggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATg                                >  2:639135/1‑143 (MQ=255)
                                                                                                               gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGtt                           >  1:626454/1‑143 (MQ=255)
                                                                                                               gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACCATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGACGCGGTGAAACTCTCCGGAATACATGCtcct                           >  2:369338/1‑140 (MQ=255)
                                                                                                                    gaTTGGCAAGGTTGTGTACATGTGTACGCGATTGGAGGGACAGGAAATGACATGCGTCGCGCTAACGATCCGCAGCGGCGAGAAAGAATTATCCAGGCCACGCTGTAGGCGGTGAAACTTTACGTAATACATGCTGTTACGcc                      <  2:6041/142‑2 (MQ=255)
                                                                                                                      ttAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAcc                    <  1:385218/143‑1 (MQ=255)
                                                                                                                            aGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTGACCCAGCGCaaaa              >  2:460236/1‑143 (MQ=255)
                                                                                                                             aattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAAt             >  1:275901/3‑143 (MQ=255)
                                                                                                                                        tttGTACACAATTGCAGAGAGAGGAAATGACATGTGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTAGGAACCGCAAAATTGATAACATTg  >  2:16969/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                             
CATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTG  >  NC_000913/887256‑887534

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