Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,764,549 T→C D202G (GAC→GGC)  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,764,5490TC100.0% 175.9 / NA 50D202G (GAC→GGC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/26);  total (24/26)

TGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTA  >  NC_000913/2764407‑2764693
                                                                                                                                              |                                                                                                                                                
tGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCg                                                                                                                                          >  1:351072/1‑151 (MQ=255)
        aTTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGttt                                                                                                                                   <  1:182447/150‑1 (MQ=255)
          tGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCg                                                                                                                                 <  2:407489/150‑1 (MQ=255)
               aTCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTg                                                                                                                           >  2:378615/1‑151 (MQ=255)
                     gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                                     >  1:134719/1‑151 (MQ=255)
                     gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                                     <  1:429686/151‑1 (MQ=255)
                             gctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCAt                                                                                                              <  2:371581/150‑1 (MQ=255)
                                  aaCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGAt                                                                                                        <  1:345187/151‑1 (MQ=255)
                                   aacACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGGCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGt                                                                                                            >  1:602387/3‑146 (MQ=255)
                                   aCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGt                                                                                                            >  2:330782/1‑146 (MQ=255)
                                       cTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTg                                                                                                    <  1:117058/150‑1 (MQ=255)
                                          aaGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAg                                                                                                <  2:491077/151‑1 (MQ=255)
                                                   tGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCg                                                                                        >  1:295368/1‑150 (MQ=255)
                                                     gTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTAt                                                                                     <  1:616518/151‑1 (MQ=255)
                                                        gTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGtt                                                                                  <  2:645696/151‑1 (MQ=255)
                                                             tAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTg                                                                              >  1:48371/1‑150 (MQ=255)
                                                                     gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAg                                                                      <  2:503079/150‑1 (MQ=255)
                                                                        atcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAg                                                                   >  1:156297/1‑150 (MQ=255)
                                                                         tcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGtt                                                                 >  2:476200/1‑151 (MQ=255)
                                                                           aGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTcc                                                               <  2:197611/151‑1 (MQ=255)
                                                                            gCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTcc                                                               >  1:162061/1‑150 (MQ=255)
                                                                             cACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTcc                                                               >  2:239385/1‑149 (MQ=255)
                                                                             cACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTcc                                                               >  2:161058/1‑149 (MQ=255)
                                                                             cACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTcc                                                               >  1:26068/1‑149 (MQ=255)
                                                                             cACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAg                                                             <  1:256554/151‑1 (MQ=255)
                                                                               cTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTgaga                                                                                                                        <  1:300034/90‑1 (MQ=255)
                                                                               cTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTgaga                                                                                                                        >  2:300027/1‑90 (MQ=255)
                                                                                tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTc                                                          >  1:317642/1‑151 (MQ=255)
                                                                                tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTc                                                          >  2:304645/1‑151 (MQ=255)
                                                                                 gtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTc                                                          <  1:605876/150‑1 (MQ=255)
                                                                                 gtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTc                                                          <  1:365080/150‑1 (MQ=255)
                                                                                   gAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTcgc                                                        <  2:329727/150‑1 (MQ=255)
                                                                                    aaTGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCCGTTCgcgc                                                      >  2:244161/1‑151 (MQ=255)
                                                                                         gTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGa                                                 <  1:299214/151‑1 (MQ=255)
                                                                                           cGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATg                                               <  1:378625/151‑1 (MQ=255)
                                                                                            gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATg                                               >  1:416011/1‑150 (MQ=255)
                                                                                            gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                              <  1:45112/151‑1 (MQ=255)
                                                                                                     cGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACtt                                      <  1:611845/150‑1 (MQ=255)
                                                                                                          tttcgcgTCACAGTATCAAGGATGGGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGAAGGCCAGCCGGATGTTCTTGATCCCTAGAAGTTCCAGATCGCGCAGGGATGAAAACACTTCCCATg                                >  2:448598/4‑151 (MQ=255)
                                                                                                          ttGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATg                                >  1:594752/1‑151 (MQ=255)
                                                                                                           tGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGt                               >  1:540554/1‑151 (MQ=255)
                                                                                                            gcgcGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGt                               <  1:116178/150‑1 (MQ=255)
                                                                                                             cgcgTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGt                               >  1:118202/1‑149 (MQ=255)
                                                                                                             cgcgTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGt                               >  1:641481/1‑149 (MQ=255)
                                                                                                             cgcgTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAg                             <  1:278521/151‑1 (MQ=255)
                                                                                                                                     cACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGc     <  1:239390/151‑1 (MQ=255)
                                                                                                                                      aCCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCt    >  1:171954/1‑151 (MQ=255)
                                                                                                                                       ccAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCtt   <  1:337323/151‑1 (MQ=255)
                                                                                                                                       ccAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCt    <  2:145012/150‑1 (MQ=255)
                                                                                                                                        cAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTta  <  2:416000/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTA  >  NC_000913/2764407‑2764693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑