Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,778,123 G→T 9.3% intergenic (+341/+23) ypjF → / ← ypjA CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter YpjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,778,1230GT 9.3% 86.4 / 2.0 39intergenic (+341/+23)ypjF/ypjACP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter YpjA
Reads supporting (aligned to +/- strand):  ref base G (15/20);  new base T (2/2);  total (17/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.40e-01

ACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAA  >  NC_000913/2777982‑2778258
                                                                                                                                             |                                                                                                                                       
accaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGa                                                                                                                                        >  1:232163/1‑143 (MQ=255)
   aaTCATGATTGGACGGTGTAAGGACAACACCCACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACGGAtgt                                                                                                                                     >  1:216838/1‑143 (MQ=255)
         ggttGGACGGTGGAAGGACAACACCAACAAAAACAGGAAGTGAGAAGCCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                               <  2:689061/141‑1 (MQ=255)
           ttGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCtttt                                                                                                                             >  1:573404/1‑143 (MQ=255)
                  gtgtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGa                                                                                                                      <  1:309257/143‑1 (MQ=255)
                    gtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                    <  2:252047/143‑1 (MQ=255)
                      aaGGACAAGACAAACAAAAGCACGAAGTGAGAGGGCGCAGCAAAGCGCCGAGCGGACGGTGAGGGGAGATAAAAGTATACGCAAGGGAGCCGGGGCTCGTGGTAACATGAAGGCCCACAGATGTGGGCTTTTTCGTTGATGGt                                                                                                                  <  2:776843/143‑1 (MQ=255)
                      aaGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGt                                                                                                                  >  1:450612/1‑143 (MQ=255)
                          acaacaCCAACAAAAACAGGAAGTTAGAAGTCTAGGTGAAACACCCCCCAGACGGTGAGGAGACATATCAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGACAGCCCACCGATGTGCGCCTTCTCGCTCACGCTCAGa                                                                                                              >  2:897193/1‑143 (MQ=255)
                            aacaCCAACAAAAACAGGGAGGTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAAc                                                                                                            <  2:865989/143‑1 (MQ=255)
                                   acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGt                                                                                                     >  2:767721/1‑143 (MQ=255)
                                              aaGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCCGca                                                                                          >  1:684340/1‑142 (MQ=255)
                                              aaGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCt                                                                                          <  1:496528/143‑1 (MQ=255)
                                                     aaGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGt                                                                                   <  1:454306/143‑1 (MQ=255)
                                                      aGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGtt                                                                                  <  2:216838/143‑1 (MQ=255)
                                                            aGCAAAACACCGACCAGACGGTGAGGAGACATAGAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACgg                                                                            <  2:173945/143‑1 (MQ=255)
                                                            aGCAAAACACCGAACAGAGGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACgg                                                                            <  1:586642/143‑1 (MQ=255)
                                                                    acCGACCAGACGGTGAGGAGACATAGAGCGATACGCAAAGGAGCCCCTACTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCca                                                                    <  1:339671/143‑1 (MQ=37)
                                                                    acCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGCGGATCCca                                                                    >  2:470533/1‑143 (MQ=255)
                                                                            gACGGTGAGGAGACATAAAAGGCTACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCa                                                            <  2:728128/143‑1 (MQ=255)
                                                                                  gaggagACATCAAAGGCTACGCACAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                      <  1:457132/143‑1 (MQ=255)
                                                                                   gggcgaaaTAAAAAGATACGCGCAGGAGCCGCGGCTCCTGGTAACATGAACGCCCACATATGTGGGCTGTTTCGTTGAGGGTCAGAACGACCTGTTCACACCATCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGc                                                     <  2:864288/136‑1 (MQ=21)
                                                                                        aCATAGAAGTGTACGCAAAGGGGCCGCGGTTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGGCATTACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATac                                                <  2:847776/143‑1 (MQ=25)
                                                                                          aTAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccc                                              >  2:846408/1‑143 (MQ=255)
                                                                                                 gATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATAccccacccca                                       >  1:562047/1‑143 (MQ=255)
                                                                                                  agaCGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCaa                                      <  1:659176/141‑1 (MQ=255)
                                                                                                     cGCAAAGGAGCCGCGGCTCCTGGACTGGTGGAGTCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGAt                                   <  1:571910/143‑1 (MQ=18)
                                                                                                      gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATg                                  <  1:573639/143‑1 (MQ=255)
                                                                                                        aaaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAGGACACACCGGCACTATAGCTATACCCCACCCCCCGATGcc                                >  2:362675/1‑143 (MQ=25)
                                                                                                              gCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTCTACCCCACCCCAACATGCCCGcgca                          >  2:874401/1‑140 (MQ=255)
                                                                                                                cgcgGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                        <  1:767721/143‑1 (MQ=255)
                                                                                                                      tCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActg                  >  1:244384/1‑143 (MQ=255)
                                                                                                                           gTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCCCCCCAAGCTGCCCGCTTAACGTACTGCTGaa             >  2:798004/1‑143 (MQ=255)
                                                                                                                              aCATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCACAAGATGCCCGCTTCACGTACTGCTGAATGa          >  1:348676/1‑143 (MQ=255)
                                                                                                                               cgaaaaaGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAg         >  2:661362/5‑143 (MQ=255)
                                                                                                                                aTGAAAGCCCACCGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCATTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAgg        <  1:756051/143‑1 (MQ=255)
                                                                                                                                   aaaGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACCATGGCTCTACCCCCCCCCCAGATGCCCGCTTAACGTACTGCTGAATGAGGCtt     <  1:677179/143‑1 (MQ=255)
                                                                                                                                    aaGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCCTACCGTACTGCTGACTGAGGCctt    >  2:689262/1‑143 (MQ=17)
                                                                                                                                      gCCCACAGATGTGGGCTTTTTCTTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTaa  <  1:645704/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                       
ACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAA  >  NC_000913/2777982‑2778258

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: