Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,764,096 A→T I353N (ATC→AAC)  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,0960AT98.0% 191.0 / NA 49I353N (ATC→AAC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (23/25);  minor base C (1/0);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCAT  >  NC_000913/2763958‑2764240
                                                                                                                                          |                                                                                                                                                
aaCTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGc                                                                                                                                       >  1:189662/1‑150 (MQ=255)
 aCTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCaa                                                                                                                                     <  2:620874/151‑1 (MQ=255)
 aCTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCa                                                                                                                                      >  2:364333/1‑150 (MQ=255)
    cGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAAc                                                                                                                                    <  1:116801/149‑1 (MQ=255)
                tCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCa                                                                                                                       >  2:417046/1‑150 (MQ=255)
                   aaaGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCaaaa                                                                                                                    >  1:122323/1‑150 (MQ=255)
                      gCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTc                                                                                                                <  1:720506/151‑1 (MQ=255)
                                 tCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTc                                                                                                       <  2:562313/149‑1 (MQ=255)
                                 tCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCa                                                                                                      <  2:242615/150‑1 (MQ=255)
                                      aTCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGc                                                                                                >  1:603253/1‑151 (MQ=255)
                                        cgcgGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCg                                                                                               >  2:539644/1‑150 (MQ=255)
                                         gcggcgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGAt                                                                                             >  2:437014/1‑151 (MQ=255)
                                           ggcgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGAttt                                                                                           >  2:7592/1‑151 (MQ=255)
                                              gAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGa                                                                                         >  1:144690/1‑150 (MQ=255)
                                                     atGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGa                                                                                  <  1:615647/150‑1 (MQ=255)
                                                       gCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAAt                                                                                <  2:744674/150‑1 (MQ=255)
                                                             tttCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTa                                                                          <  2:511783/150‑1 (MQ=255)
                                                             tttCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTa                                                                          <  2:175021/150‑1 (MQ=255)
                                                               tCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGTTGCCGGCTCGATGAGTGCAACAACACTCCACCCAAAAATTCGTTGACGGCATGGAGCGAATTGAAGACGGAATTCTTTaat                                                                        >  1:553914/1‑148 (MQ=255)
                                                                 aaGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCg                                                                      >  1:546646/1‑150 (MQ=255)
                                                                 aaGAATTTTTGCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGt                                                                     <  2:500046/151‑1 (MQ=255)
                                                                  aGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCCGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGtt                                                                    >  2:584568/1‑151 (MQ=255)
                                                                   gAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtc                                                                   <  2:617094/151‑1 (MQ=255)
                                                                    aaTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtct                                                                  <  2:419178/151‑1 (MQ=255)
                                                                       tttttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGc                                                               >  1:651333/1‑151 (MQ=255)
                                                                        ttttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGc                                                               <  1:372014/150‑1 (MQ=255)
                                                                          ttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCttt                                                            >  1:655597/1‑151 (MQ=255)
                                                                            cGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCtttt                                                           <  2:490601/150‑1 (MQ=255)
                                                                            cGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTc                                                          <  2:614525/151‑1 (MQ=255)
                                                                             ggCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCa                                                         <  1:576408/151‑1 (MQ=255)
                                                                               ccccTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCaat                                                       >  2:147861/1‑151 (MQ=255)
                                                                                         ttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAg                                              <  2:460066/150‑1 (MQ=255)
                                                                                          tCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGtt                                            >  1:487076/1‑151 (MQ=255)
                                                                                           ccACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGtt                                            <  2:607578/150‑1 (MQ=255)
                                                                                              cGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATg                                        >  2:517518/1‑151 (MQ=255)
                                                                                              cGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGGGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATg                                        >  2:679838/1‑151 (MQ=255)
                                                                                               gTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGa                                       <  2:34328/151‑1 (MQ=255)
                                                                                               gTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGa                                       <  2:355213/151‑1 (MQ=255)
                                                                                                      aTGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCa                                 <  2:299488/150‑1 (MQ=255)
                                                                                                        gCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGa                               <  1:393313/150‑1 (MQ=255)
                                                                                                          aaGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACAt                            >  1:9440/1‑151 (MQ=255)
                                                                                                                tAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCttt                      >  2:440420/1‑151 (MQ=255)
                                                                                                                      ggctTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTc                 <  1:239761/150‑1 (MQ=255)
                                                                                                                          ttCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGga              <  1:122305/149‑1 (MQ=255)
                                                                                                                                gtgCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAAcc      >  2:510108/1‑151 (MQ=255)
                                                                                                                                   ctgcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCt     >  2:665346/1‑149 (MQ=255)
                                                                                                                                    tgcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCAt  >  2:516700/1‑151 (MQ=255)
                                                                                                                                    tgcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCAt  >  1:460252/1‑151 (MQ=255)
                                                                                                                                     gcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCAt  <  2:163141/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
AACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCAT  >  NC_000913/2763958‑2764240

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