Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,822 A→T 13.3% intergenic (+123/‑53) ecnB → / → gdx bacteriolytic entericidin B lipoprotein/guanidinium exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,8220AT13.3% 61.9 / 2.3 30intergenic (+123/‑53)ecnB/gdxbacteriolytic entericidin B lipoprotein/guanidinium exporter
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base T (2/2);  total (15/15)
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.80e-01

AACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTT  >  NC_000913/4376681‑4376960
                                                                                                                                             |                                                                                                                                          
acccaaaGCGCAGCAATATGCAAGCACGGTACGACAGCGGTGTCGTGCCGTTTGGTTTTTCTGCGATAGTCACAAAGGTATTAGTTGAAATGCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                                                                                                                           <  2:210255/139‑1 (MQ=255)
   gAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAATTATCCGTTCAACCATCAGCATTGCCGGACGACATg                                                                                                                                        >  2:119887/1‑143 (MQ=255)
    acatcgcAGCAATAATCTAGAACGGTACGACAGCGGTTTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAACGTAATAGTCGAAATTCCCCTGCCACGTGGCAAAATATCCCTTCAACCATCAGCTTTGCAGGACGACCTGc                                                                                                                                       <  2:114796/139‑1 (MQ=255)
       gcgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaac                                                                                                                                    >  1:206832/1‑142 (MQ=255)
        agcagcaATAGGCAATAACGGTACGACAGCTGTGTCGTGCCGTTGGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAATCATCAGCTTTGCAGGTCGTCCTGCaaag                                                                                                                                   <  2:369239/142‑2 (MQ=255)
           agcaataagcaataACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAACAGTTGGAATTCCCCTGCCACCTGTCATAATCTCCGTTCATCCATCAGCTATGCCGTACCGCATGCAAACGcg                                                                                                                                >  2:577175/1‑142 (MQ=255)
                  agcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTc                                                                                                                         >  1:565824/1‑143 (MQ=255)
                      ataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCg                                                                                                                     <  1:197600/143‑1 (MQ=255)
                      agaaCGGTACGAGAGTGCGGTCGTGTCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTGCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCg                                                                                                                     <  2:773305/141‑1 (MQ=255)
                            gTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCGCTTTTCACCGGGGGCg                                                                                                               >  2:599882/1‑143 (MQ=255)
                                gACAGCTGTGTCGTGCCGTTTGTTTTTTCGTCGATAGACAAAAAGGTAATAGTGGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGccc                                                                                                           <  1:700439/143‑1 (MQ=255)
                                      tgtgTCGTGGCGTTTGTTTTGTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtc                                                                                                     <  2:9648/143‑1 (MQ=255)
                                             tcccGTGTGTTTTGCCTGCGATAGTCACAAAGGTAATAGTCGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCACCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAg                                                                                              <  2:336882/141‑1 (MQ=255)
                                              gCCGTTTGTTTTTTCTTCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGTCCTGCAAACGCCTCTTTTCTCCGGGGACGGCCCCAATTCTCCGGAGc                                                                                             <  2:279662/143‑1 (MQ=255)
                                                 gtttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTg                                                                                          >  1:181889/1‑143 (MQ=255)
                                                   ttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                        <  1:59366/143‑1 (MQ=255)
                                                   ttgtttTTTCCGCGATAGTCACACAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                        <  1:59364/143‑1 (MQ=255)
                                                             tGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGa                                                                              <  1:730244/143‑1 (MQ=255)
                                                               cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAGCGGGGACGGCCCCGATTCTCCGGAGCCGGATATGTCCGGGCtt                                                                            >  2:209022/1‑143 (MQ=255)
                                                                     tCACAAAGGTAATAGTTGACATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                      <  2:569705/143‑1 (MQ=255)
                                                                       acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGt                                                                    >  2:447817/1‑143 (MQ=255)
                                                                          aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAGCGCCTCGTTTCGCCGGGGGCGGCAATAATGCTCCGGAGCAAGATAGCTCCTGGGTTATGTGAGTTAt                                                                 <  2:17723/143‑1 (MQ=255)
                                                                                      gAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATGTGTCCTGGATTATCTTAGTTATTGCACGTctgcc                                                     >  2:244460/1‑142 (MQ=255)
                                                                                                        ggCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGt                                   >  1:100212/1‑143 (MQ=255)
                                                                                                           aaaaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCCATTCTCCGGAGCCTGAAATGTCATGAATTATCTTAGTTACTGCCCGTCTGCAGGCAATGGTATTGGCCGTTgg                                >  2:159567/1‑143 (MQ=255)
                                                                                                                atCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGa                           <  2:100212/143‑1 (MQ=255)
                                                                                                                       cAACCACCAGCTGTGCAGGCCGACCTGCAACCGCCGCCTTTCACCGGCGACGGCCCCAATTCTCCGGAGCCCGATATGTCCTGGATTAGCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATAc                    <  1:36654/143‑1 (MQ=255)
                                                                                                                           cATGAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCAc                >  1:772125/1‑143 (MQ=255)
                                                                                                                             tCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACgg              >  1:734313/1‑143 (MQ=255)
                                                                                                                                cgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGCGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCtt           >  2:389923/3‑143 (MQ=255)
                                                                                                                                  tttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGTTATGTACTGGAGTGGCCTAGTTATTGCTGGTCTGCTGGGAGTGGTATGGGCCGGTGGCCTGAAATATACCCCCggggtta         >  2:685361/1‑138 (MQ=255)
                                                                                                                                         gacgacCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGttt  <  1:70436/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
AACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTT  >  NC_000913/4376681‑4376960

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: