Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,763,869 (T)7→8 coding (1285/1617 nt) yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,763,8621.T94.7% 122.9 / ‑0.1 38coding (1292/1617 nt)yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (16/20);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

TAGATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGG  >  NC_000913/2763732‑2763995
                                                                                                                                   |                                                                                                                                     
tagATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGAtttttttt                                                                                                                               >  2:104740/1‑138 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCgg                                                                                                                        >  2:825624/1‑139 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCgg                                                                                                                        >  2:509913/1‑139 (MQ=255)
                 ccACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGc                                                                                                               <  1:929648/139‑1 (MQ=255)
                           aacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAAt                                                                                                     >  2:342293/1‑139 (MQ=255)
                                         gATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                       >  1:551863/1‑139 (MQ=255)
                                         gATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                       <  1:843611/139‑1 (MQ=255)
                                             gAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGaa                                                                                                       <  1:2664/119‑1 (MQ=255)
                                             gAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGaa                                                                                                       >  2:2664/1‑119 (MQ=255)
                                               ggTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATccc                                                                                 <  2:93326/139‑1 (MQ=255)
                                                  gAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTa                                                                              >  1:84883/1‑139 (MQ=255)
                                                     cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTc                                                                                                  <  2:402892/116‑1 (MQ=255)
                                                     cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTc                                                                                                  >  1:402892/1‑116 (MQ=255)
                                                                  gCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGAc                                                              >  1:847154/1‑139 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           <  1:104740/139‑1 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           <  1:206262/139‑1 (MQ=255)
                                                                      ccATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          <  2:8786/139‑1 (MQ=255)
                                                                      ccATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          <  1:509913/139‑1 (MQ=255)
                                                                                   cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGaa                                                                                                       >  2:934394/1‑81 (MQ=255)
                                                                                   cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGaa                                                                                                       <  1:934394/81‑1 (MQ=255)
                                                                                   cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAg                                             >  2:784856/1‑139 (MQ=255)
                                                                                           tCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAAc                                     <  1:181474/139‑1 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGa                                                               <  1:72311/103‑1 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGa                                                               >  2:72311/1‑103 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                           <  1:1037930/139‑1 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                           <  1:22269/139‑1 (MQ=255)
                                                                                                             ttATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCaa                   <  1:848769/139‑1 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa             >  2:249642/1‑139 (MQ=255)
                                                                                                                    cTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt            >  1:806740/1‑139 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                                                                                        >  1:734201/1‑62 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                                                                                        <  2:734201/62‑1 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc           <  2:847154/139‑1 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc           >  2:487381/1‑139 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc           <  2:523705/139‑1 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc           <  1:249642/139‑1 (MQ=255)
                                                                                                                     tttAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc           <  1:784856/139‑1 (MQ=255)
                                                                                                                             cGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          <  2:291649/83‑1 (MQ=255)
                                                                                                                             cGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          >  1:291649/1‑83 (MQ=255)
                                                                                                                              gAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTgg  >  2:910913/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
TAGATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGG  >  NC_000913/2763732‑2763995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: