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.T98.2% 183.9 / ‑3.5 54coding (1292/1617 nt)yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (29/24);  total (30/24)
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.33e-01

TAGATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGG  >  NC_000913/2763732‑2763995
                                                                                                                                   |                                                                                                                                      
tagATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGAtttttttt                                                                                                                               <  1:973401/139‑2 (MQ=255)
   aTTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGAttttttttgcg                                                                                                                            <  2:904848/139‑5 (MQ=255)
     tAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcg                                                                                                                           <  2:1048055/139‑1 (MQ=255)
     tAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcg                                                                                                                           >  1:539517/1‑139 (MQ=255)
      agatTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGCGGTTAGCATCGAAAGGTTTGCTGCTATCAAGAACTTCGTCAATGTCTTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcgc                                                                                                                          <  2:807945/139‑1 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCgg                                                                                                                        >  1:880594/1‑139 (MQ=255)
          tCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTa                                                                                                                      >  1:299701/1‑139 (MQ=255)
          tCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTa                                                                                                                      >  1:50132/1‑139 (MQ=255)
             tAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCg                                                                                                                   >  1:918762/1‑139 (MQ=255)
                  cACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCa                                                                                                              >  1:630676/1‑139 (MQ=255)
                    cGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAAt                                                                                                            >  1:798943/1‑139 (MQ=255)
                             caAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                                   <  1:853165/139‑1 (MQ=255)
                              aaaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTc                                                                                                  <  2:50132/139‑1 (MQ=255)
                               aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                                 <  2:798943/139‑1 (MQ=255)
                                aCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCt                                                                                                >  1:745366/1‑139 (MQ=255)
                                     aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcgc                                                                                                                          <  1:824277/108‑1 (MQ=255)
                                     aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcgc                                                                                                                          >  2:824277/1‑108 (MQ=255)
                                     aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                                 >  2:33582/1‑133 (MQ=255)
                                     aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                                 <  1:33582/133‑1 (MQ=255)
                                     aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTAc                                                                                           >  1:483373/1‑139 (MQ=255)
                                         gttcaaGGTGAGCATCGAAATGTTTGCTTCCATCAAGACTTTCGTCACTCTCGTGGCGAACCAGGTAGTTATCGTCTTTAATGACGACAATTTTTTTTGCGCGTCGTACCGACGCACTCCCAAATGTCCTCGTACAACt                                                                                       <  2:686835/134‑1 (MQ=255)
                                         gATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                       >  2:527188/1‑139 (MQ=255)
                                          aTCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                                 >  2:1006425/1‑128 (MQ=255)
                                          aTCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                                 <  1:1006425/128‑1 (MQ=255)
                                            cGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGAt                                                                                    >  1:1105466/1‑139 (MQ=255)
                                                   aGCATCGAAAGGTTGACTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                       <  2:588198/129‑1 (MQ=255)
                                                   aGCATCGAAAGGTTGACTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTATTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                       >  1:588198/1‑129 (MQ=255)
                                                     cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcg                                                                                                                           <  1:747127/91‑1 (MQ=255)
                                                     cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTgcgcg                                                                                                                           >  2:747127/1‑91 (MQ=255)
                                                                  gCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGAt                                                                                                                 <  1:247706/88‑1 (MQ=255)
                                                                  gCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGAt                                                                                                                 >  2:247706/1‑88 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           <  1:228096/139‑1 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           <  2:1105466/139‑1 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           <  2:73302/139‑1 (MQ=255)
                                                                      ccATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          <  2:918762/139‑1 (MQ=255)
                                                                      ccATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          <  2:299701/139‑1 (MQ=255)
                                                                                   cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGc                                                                                                               >  2:618779/1‑73 (MQ=255)
                                                                                   cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGc                                                                                                               <  1:618779/73‑1 (MQ=255)
                                                                                        aTGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCAcc                                        >  1:339126/1‑139 (MQ=255)
                                                                                           tCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAAc                                     <  2:745366/139‑1 (MQ=255)
                                                                                                   aaGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGaa                             >  1:629064/1‑139 (MQ=255)
                                                                                                    aGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAg                            >  1:984698/1‑139 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           >  2:779815/1‑107 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                           <  1:779815/107‑1 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                           >  1:533434/1‑139 (MQ=255)
                                                                                                     gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                           <  2:984698/139‑1 (MQ=255)
                                                                                                           aGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          >  2:205229/1‑102 (MQ=255)
                                                                                                           aGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                          <  1:205229/102‑1 (MQ=255)
                                                                                                              tATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCaa                   >  1:249769/1‑138 (MQ=255)
                                                                                                              tATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCaa                   <  2:249769/138‑1 (MQ=255)
                                                                                                                  gTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTg              >  2:146126/1‑139 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa             >  1:259466/1‑139 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa             >  2:99652/1‑139 (MQ=255)
                                                                                                                    cTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt            <  2:630676/139‑1 (MQ=255)
                                                                                                                    cTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt            >  2:853919/1‑139 (MQ=255)
                                                                                                                              gAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCACGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCGAAGCTGATCACCTTCTgg  >  2:783475/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
TAGATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGG  >  NC_000913/2763732‑2763995

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